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ABSTRACT 


A  physical  understanding  and  consequent  mathematical 
modeling  of  RF  energy  in  naval  indoor  environments  is  of  vital 
importance  to  the  usability  and  effectiveness  of  communication 
systems  used  by  Navy.  Over  the  last  few  years,  there  is  a 
growing  interest  in  placing  Wireless  Local  Area  Networks  (WIANs) 
in  ships  and  submarines.  Especially  large  ships  yet  to  be 
constructed,  are  designed  with  increased  electronic  systems  but 
limited  personnel.  Reliable  electronic  systems  will  be  crucial 
for  efficient  ship  operation  and  survivability. 

This  thesis  investigates  the  feasibility  of  deploying  a 
physical  model  called  Numerical  Electromagnetic  Code-Basic 
Scattering  Code  (NEC-BSC)  to  simulate  confined  naval  compartments 
in  the  2.4  GHz  Industrial  Scientific  Medical  (ISM)  band.  More 
specifically,  the  coverage  area,  the  number  and  positions  of 
transmitters  and  observation  points  and  the  statistics  of  Radio 
Frequency  (RF)  signal  distribution  were  described.  The  area 
specifically  targeted  for  this  research  was  a  typical  two-story 
missile  room.  Additionally,  some  important  conclusions  regarding 
the  validity  of  NEC-BSC  for  indoor  applications  are  presented  and 
some  recommendations  for  future  research  are  provided. 
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EXECUTIVE  SUMMARY 


Indoor  wireless  systems  are  already  used  in  a  huge 
variety  of  office,  residential  and  even  military 
environments.  Two  types  of  models,  the  empirical  and  the 
physical  ones,  are  powerful  tools,  capable  of  simulating  a 
given  confined  area  and  collecting  significant  databases. 

The  Numerical  Electromagnetic  Code-Basic  Scattering 
Code  (NEC-BSC) ,  developed  at  Ohio  State  University  under 
government  contract,  is  a  well-known  physical  model  in 
electromagnetic  field  theory.  Over  the  last  ten  years  it 
has  been  of  great  use  in  analyzing  scattering  and  radiation 
problems  and  visualizing  the  flow  of  energy  from  antennas. 
Given  that  the  code's  applicability  has  recently  been 
proved  reliable  for  complex  buildings  [37] ,  we  tested  its 
performance  to  a  typical  missile  room  of  a  military  ship  in 
the  2.4  GHz  ISM  band. 

At  first,  the  case  of  the  electric  field  returned  by  a 
smooth  perfectly  conducting  convex  surface  was 
investigated.  Then,  the  results • were  compared  to  the  ones 
obtained  by  the  reliable  software  called  "PATCH",  within 
the  framework  and  limitations  of  the  Method  of  Moments 
(MOM)  .  The  comparison  offered  encouraging  results.  In 
addition,  other  simple  geometries  were  simulated  in  order 
to  gain  familiarity  with  the  code  and  its  limitations.  The 
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goals,  however,  were  to  construct  a  typical  missile  room 
with  material  (metallic)  uniformity,  examine  the  electric 
field  strength  coverage  maps  and  perform  a  statistical 
analysis  of  the  electric  field  strength  variation. 

To  this  end,  the  antenna  was  placed  in  four  different 
positions  and  data  were  collected.  Then  the  excess  "room 
gain"  over  free  space  of  the  same  volume  was  calculated 
within  certain  distances  from  the  antenna,  called  "distance 
groups."  Next,  the  probability  density  function  (pdf)  and 
cumulative  distribution  function  (cdf)  of  the  variation  of 
the  "room  gain"  variable  was  introduced.  It  appeared  that 
the  key  parameter  in  defining  the  mean  value  of  the  excess 
gain  was  the  relative  number  of  observation  points  that 
were  shadowed  by  the  obstacles  and  thus  experience  non- 
line-of-sight  propagation  attenuation.  In  any  case, 
however,  the  distribution  followed  the  lognormal  one  and 
its  pdf  and  cdf  obeyed  the  Gaussian  ones.  Finally,  the 
mean  value  of  the  "room  gain"  ranged  from  -0.1  dB  to  +  2.91 
dB. 
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I. 


INTRODUCTION 


A.  WIRED  VS  WIRELESS  CC»4MCINICATION 

The  last  century  has  seen  the  rise  of  telephony,  as 
well  as  the  successful  deployment  of  wireless 
communications.  •  The  introduction  of  these  inventions 
underscored,  without  doubt,  a  fundamental  change  in  the  way 
people  communicated  and  was  a  critical  step  toward 
shattering  the  barriers  of  time  and  space.  To  date, 
however,  the  hindrance  of  location  has  not  been  conquered, 
since  most  individuals  are  still  confined  to  wired 
equipment  for  communications  [1] . 

Rapid  changes  in  word  communications  have  been 

prompted  by  demand  for  more  reliable  and  convenient  ways  to 

communicate.  Especially  during  the  past  decade,  wireless 

services  proved  to  be  cost  effective  solutions  for 

providing  data  due  to  their  inherent  advantages  over  the 

conventional,  wired  services.  These  advantages  include 

simplicity  and  ease  of  use,  integrity  and  reliability, 

compatibility  with  existing  networks  without  the  need  for 

expensive  and  time-consuming  rewiring.  In  addition,  there 

is  more  immunity  to  natural  and  man  made  disasters, 

flexibility  of  shifting  terminals  and  equipment  around, 

relatively  easy  maintenance,  providing  services  for  a 
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larger  number  of  users,  roaming  and,  finally,  elimination 
of  wiring  around  and  inside  buildings.  The  introduction  of 
successful  wireless  systems  is  comparable  in  significance 
with  that  of  any  major  invention  made  during  the  20th 
century.  The  deployment  of  wireless  systems  is  increasing, 
while  the  wired  communications  is  holding  steady  if  not 
decreasing. 

All  the  aforementioned  advantages  of  using  wireless 
communication  mainly  in  small  or  large  offices,  public 
services  buildings  and  even  military  environments,  have 
built  a  potential  market  of  a  very  large  number  of  users. 
In  order  to  meet  these  demands,  new  systems  have  been  and 
are  being  designed  to  provide  services  of  the  best  possible 
quality  in  indoor  and  outdoor  environment  [2] . 

B.  WIKELESS  PROPAGATION  BACKGROUND 

Radio  communication,  to  date,  is  enjoying  the  fastest 
growth  period  in  history,  since  enabling  technologies 
permit  wide  spread  deployment.  During  the  past  decade, 
radio  communications  business  has  grown  by  orders  of 
magnitude,  fueled  by  large-scale  circuit  integration, 
digital  and  radio  frequency  (RF)  circuit  fabrication 
improvements  and  other  miniaturization  technologies  which 
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make  radio  equipment  smaller,  less  heavy  and  more 
affordable  and  reliable  [3] . 

The  Ultra  High  Frequency  (UHF)  electromagnetic 
spectrum  ranges  from  300  MHz  to  3  GHz.  There  is  much  more 
frequency  space  available  for  communications  in  this  band 
than  in  the  entire  spectrum  below.  Moreover,  high 
directivity  obtainable  at  microwave  frequencies  makes  reuse 
of  these  frequencies  possible, in  the  same  area. 

The  most  basic  model  of  radio  wave  propagation 
involves  the  ideal  "free  space"  model  which  treats  the 
region  between  the  transmit  and  receive  antennas  as  being 
free  of  all  objects  that  might  absorb,  refract,  reflect, 
diffract  or  scatter  radio  frequency  energy.  It  thus 
follows  that,  within  this  region,  the  medium  behaves  as  a 
perfectly  homogeneous  and  nonabsorbent  medium. 
Furthermore,  earth  is  treated  as  being  infinitely  far  away 
from  the  propagating  signal  (or,  equivalently,  as  having  a 
reflection  coefficient  that  is  negligible) . 

In  free  space  the  spreading  of  RF  energy  creates 
apparent  attenuation  between  the  transmitter  and  receiver 
according  to  an  inverse-square  law.  In  typical 
broadcasting  applications  radio  waves  emanate  from  a  point 
source  of  radio  energy,  traveling  in  all  directions  along 
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straight  lines,  filling  the  entire  spherical  volume  of 
space  with  radio  energy  that  varies  in  strength  with  a 
l/(range)^  rule  (or  20  dB  per  decade  increase  in  range). 
The  received  power,  expressed  in  terms  of  transmitted 
power,  is  attenuated  by  a  factor  called  (free  space)  path 
loss . 

Calculating  free  space  loss  is  quite  simple.  Consider  a 
transmitter  with  power  Pt  coupled  to  an  antenna,  which 
radiates  equally  in  all  directions  (an  isotropic  antenna) . 
At  a  distance  d  from  the  transmitter,  the  radiated  power  is 
distributed  uniformly  over  an  area  of  4d^  (i.e.-the  surface 
area  of  a  sphere  of  radius  d) ,  so  that  the  power  flux 
density  is: 

s=  ^ ,  Watts/m^  .  (1.1) 

The  transmission  loss  then  depends  on  how  much  of  this 
power  is  captured  by  the  receiving  antenna.  If  the  capture 
area,  or  effective  aperture  of  this  antenna  is  Ag  then  the 
power,  which  can  be  delivered  to  the  receiver  (assuming  no 
mismatch  or  feed  line  losses),  is  simply: 
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P^  =  sA^  Watts. 


(1.2) 


For  the  hypothetical  isotropic  receiving  antenna  and  a 
wavelength  X  we  have: 


A 


e 


(1.3) 


Combining  all  equations,  we  obtain: 


ATid 


Watts. 


(1.4) 


The  free  space  path  loss  between  isotropic  antennas  is 
then  Pt/Pr.  For  this  case  of  idealized  propagation,  received 
signal  power  is  very  predictable. 

In  reality,  however,  propagation  rarely  follows  the 
free  space  model.  Three  major  mechanisms,  which  are  going 
to  be  examined  in  the  next  pages,  dominate  in  the  radio 
propagation:  (1)  reflection,  (2)  diffraction  and  (3) 
scattering.  All  three  of  these  phenomena  result  in  radio 
signal  distortions,  additional  losses,  and  cause  signals  to 
fade. 
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C.  THE  PHYSICS  OF  PROPAGATION 

The  phenomena,  which  govern  radio  propagation,  can  be 
very  complex  and  diverse.  However,  we  can  categorize  them 
according  to  the  three  major  mechanisms  of  object 
interaction  (reflection,  diffraction  and  scattering) . 

Reflection  occurs  when  a  propagating  electromagnetic 
wave  impinges  upon  an  object,  which  has  very  large 
dimensions  compared  to  the  wavelength  of  the  propagating 
wave.  Reflections  occur  from  the  surface  of  the  earth  and 
from  buildings  and  walls,  and  reflected  waves  can  interfere 
constructively  or  destructively  at  the  receiver  [3] . 

Diffraction  occurs  when  the  radio  path  between  the 
transmitter  and  receiver  is  obstructed  by  an  impenetrable 
surface.  The  secondary  waves  resulting  from  the 
obstructing  surface  are  formed  throughout  the  space  and 
even  behind  the  obstacle,  giving  rise  to  an  apparent 
bending  of  waves  around  the  obstacle,  enabling  reception 
even  when  a  line-of-sight  (LOS)  path  does  not  exist  between 
transmitter  and  receiver.  Diffraction  is  assumed  to  be  a 
strictly  local  phenomenon;  this  means  that  the  diffracted 
field  depends  only  on  the  power  of  the  incident  field  at 
the  exact  diffraction  point,  and  on  the  local  geometry  of 
the  diffracting  object  [3] . 
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At  high  frequencies,  diffraction,  like  reflection, 
depends  on  the  geometry  of  the  object,  as  well  as  the 
amplitude,  phase,  and  polarization  of  the  incident  wave  at 
the  point  of  diffraction.  Because  of  diffraction  RF  energy 
can  be  received,  eventhough  the  source  is  "shadowed"  by  an 
obstacle  [3]  , 

Scattering  occurs  when  the  medium  through  which  the 
wave  travels  contains  objects  ("obstacles")  with  dimensions 
that  are  small  compared  to  the  wavelength,  and  where  the 
number  of  obstacles  per  unit  volume  is  large.  Scattered 
waves  are  produced  by  rough  surfaces,  atmospherics  (rain 
and  snow),  or  by  other  irregularities  in  the  channel. 
Foliage,  street  signs,  and  lampposts  induce  scattering  in  a 
mobile  communications  system.  Figure  1  shows  the  different 
propagation  mechanisms  in  the  presence  of  an  arbitrary, 
impenetrable  obstacle. 
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multiple  signal  propagation  paths  exist,  caused  by  whatever 
phenomenon,  the  actual  received  signal  level  is  a  vector 
sum  of  the  all  signals  incident  from  all  directions.  Some 
signals  will  aid  the  direct  path  (assuming  that  the  direct 
path  exists),  while  other  signals  will  subtract  (or  tend  to 
vector  cancel)  from  the  direct  signal  path.  The  total 
composite  phenomenon  is  called  multipath. 

As  a  mobile  receiver  moves  through  a  specific  indoor 
space,  the  three  phenomena  affect  the  instantaneous 
received  signal  in  different  ways.  For  a  clear  LOS  path  to 
the  transmitter,  diffraction  and  scattering  will  not 
dominate.  On  the  other  hand,  in  an  environment  where  there 
are  many  obstacles  of  different  sizes  and  constituent 
materials  and  the  direct  channel  is  bloclced  between  the 
receiver  and  the  transmitter,  diffraction  as  well  as 
scattering  are  most  likely  to  dominate  [5] . 

Two  kinds  of  multipath  exist:  (1)  specular  multipath 
arising  from  discrete,  coherent  reflections  from  smooth 
metal  surfaces,  and  (2)  diffuse  multipath  arising  from 
diffuse  scatterers  and  sources  of  diffraction  (the  visible 
glint  of  sunlight  off  a  choppy  sea  is  such  an  example)  . 
Both  forms  of  multipath  are  undesired  for  radio 
communications.  Diffuse  multipath  provides  a  sort  of 
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background  "noise"  level  of  interference,  while  specular 
multipath  can  actually  cause  complete  signal  outages  and 
radio  "dead  spots"  within  a  building. 

Another  phenomenon,  called  scintillation,  describes 
the  rapid  fluctuation  of  the  phase  and  intensity  of  a  radio 
signal  that  has  passed  through  the  earth's  ionosphere, 
typically  on  a  satellite-to-ground  propagation  channel. 
This  phenomenon  is  similar  to  the  twinkling  of  stars  in  the 
night  sky,  except  in  this  case  the  fluctuations  are  caused 
by  small-scale  variations  (irregularities)  in  the 
ionospheric  plasma  density  along  the  propagation  path 
followed  by  the  signal . 

D.  FUNDAMENTALS  OF  ANTENNAS 

1 .  Radiation  Pattern 

The  antenna' s  radiation  pattern  is  the  graphical 
representation  of  the  radiation  characteristics  of  the 
antenna  as  a  function  of  space  coordinates.  The  most  usual 
representation  is  the  three-dimensional  spatial 
distribution  of  radiated  energy  along  an  observer's 
position  at  a  constant  distance  from  the  antenna.  The 
graph  of  the  spatial  variation  of  the  electric  field  along 
that  path  is  called  a  field  pattern. 
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An  isotropic  antenna  is  defined  as  a  hypothetical 
lossless  antenna  that  radiates  or  receives  equally  in  all 
directions.  Clearly,  an  isotropic  antenna  is  a  fictitious 
entity,  since  even  the  simplest  antenna  has  some  degree  of 
directivity.  However,  an  isotropic  radiator  represents  a 
convenient  reference  for  comparing  the  directional 
properties  of  physical  antennas.  On  the  other  hand,  a 
directional  antenna  radiates  or  receives  electromagnetic 
energy  more  effectively  in  some  directions  than  in  others. 
As  an  example,  the  radiation  pattern  of  a  half-wave  linear 
dipole  is  shown  in  Figure  2. 


Figure  2.  Power  Pattern  of  a  Half-Wave  Dipole  {from  [6]). 
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2. 


Near  Field  and  Far  Field  Regions 


The  antenna's  electromagnetic  (EM)  field  can  be 
classified  according  to  the  relative  (in  terms  of  operating 
wavelength)  distance  to  the  antenna  as: 

•  The  reactive  near  field  region  adjacent  to  the 
antenna.  In  this  case,  the  field  patterns  change 
rapidly  and  reactive  energy  is  stored  around  the 
antenna 

•  The  radiating  near  field  region  (for  large 
antennas  known  as  the  Fresnel  region) .  Again, 
the  field  patterns  change  rapidly,  but  radiating 
energy  oscillates  towards  the  antenna. 

•  The  furthermost,  the  far  field  region,  commonly 
referred  to  as  the  Fraunhofer  region,  has  no 
reactive  energy  and  is  represented  by  a  spherical 
wave  (for  any  realistic  antenna  of  finite  size) 
propagating  away  from  the  antenna. 

Certain  behavior  characteristics  of  EM  fields  dominate 
at  one  range  from  the  radiating  antenna,  while  others, 
completely  different,  can  dominate  at  another  distance.  In 
order  to  categorize  characteristics  of  electromagnetic 
fields  as  a  function  of  distance  from  the  radiating  source, 
we  define  the  following  boundaries:  the  "Near-Field 
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boundary, "  the  "Far-Field  boundary"  and  the  "Transition 
Zone  boundary, "  which  has  a  combination  of  the 
characteristics  found  in  both  fields.  The  regional 
boundaries  are  usually  categorized  as  a  function  of  the 
operating  wavelength  {X)  . 

We  must  state,  however,  that  these  regions  categorize 
behaviors,  which  may  vary  even  within  each  region.  In 
addition,  the  boundaries  for  these  regions  are  approximate 
"rules  of  thumb"  (more  precise  boundaries  can  be  defined 
based  primarily  on  antenna  type  and  antenna  size,  and  even 
then  the  experts  differ) .  Figure  3  shows  these  regions  and 
boundaries . 
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Figure  3.  Antenna  Field  Regions  for  Typical  Antennas  (from 
Ref .  [  7  ]  )  . 

When  placing  an  antenna,  it  is  of  vital  importance  to 
remove  all  objects  from  at  least  the  near  field,  since 
their  presence  can  alter  the  design  radiation  and  impedance 
characteristics  unpredictably . 

E.  OBJECTIVE  AND  APPROACH 

The  objective  of  this  thesis  is  to  evaluate  the 

Numerical  Electromagnetic  Code  -  Basic  Scattering  Code 

(NEC-BSC) ,  for  shipboard  (indoors)  wireless  propagation 

prediction.  Using  the  code,  we  will  predict  the  signal 

strength  distribution  for  wireless  communication  in  the  2.4 
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GKz-  Instrumentation-Scientific-Medical  (ISM)  band  for  a 

given  antenna  location  in  a  shipboard  compartment. 
Simulation  results  will  be  presented  graphically  and  the 
signal  strength  statistics  will  be  extracted  from  the 
resulcs. 

F.  THESIS  OUTLINE 

This  thesis  consists  of  six  chapters  and  is  organized 
as  follows.  Chapter  I  compares  the  wired  and  wireless 
communication  and  provides  the  introduction  of  the  RF 
propagation.  Chapter  II  describes  the  theoretical 
principles  of  indoor  propagation.  Chapter  III  presents  the 
NEC-BSC  software  overview.  Chapter  IV  describes  some  basic 

examples  and  the  associated  results  obtained  by  using  NEC- 
BSC  code.  Chapter  V  provides  the  development  of  a 
realistic  model  of  a  shipboard  compartment,  the  simulation 
results  and  the  statistics  for  signal  propagation. 
Finally,  chapter  VI  concludes  the  thesis  and  present 
recommendations  for  future  work. 
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II.  INDOOR  PROPAGATION  MODELS  AND  THESIS  OBJECTIVES 


A.  COMPARISON  BETWEEN  INDOOR  AND  OUTDOOR  MODELS 

Outdoor  environments  are  characterized  by  significant 
signal  strength  fluctuations  caused  by  reflection, 
diffraction,  and  scattering.  These  signal  strength 
variations  may  have'  a  range  of  30  to  40  dB  at  Ultra  High 
Frequencies  (UHF)  and  account  for  some  of  the  difficulties 
that  a  designer  of  reliable  radio  communications  systems 
has  to  cope  with.  For  example,  if  we  shift  the  receiver  by 
a  distance  of  a  half -wavelength  in  a  random  direction,  the 
received  signal  magnitude  may  change  drastically.  The 
basic  signal  attenuation  with  range  noticed  in  the  real 
world  describes  the  so-called  "large  scale"  effects,  while 
the  signal  strength  fluctuations  with  motion  are  reported 
"small  scale"  effects  [8]  . 

Outdoor  propagation  models  presented  by  Ibrahim  [Ref. 
10],  Okumura  [11]  and  Hata  [12],  have  been  developed  for 
large  cells,  whereas  for  indoor  communications  the  suitable 
cell  diameter  is  substantially  smaller  and  such  cells  are 
usually  termed  as  nanocells  and  picocells.  Hence,  the 
outdoor  models  are  inadequate  for  accurate  indoor 
propagation  prediction. 
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Indoor  propagation,  although  dominated  by  the  same 
mechanisms,  can  not  be  treated  as  an  extension  of  outdoor 
propagation.  For  example,  signal  levels  depend  on  whether 
interior  doors  and  windows  are  open  or  closed.  Also,  the 
location  of  antennas  affects  the  large-scale  phenomena, 
since  an  antenna  mounted  at  a  floor  level  in  a  partitioned 
room  receives  different  signals  than  that  same  antenna 
mounted  above  the  floor.  Moreover,  in  a  room  crowded  with 
obstacles,  all  propagation  mechanisms  become  more  and  more 
complex  and  the  .smaller  propagation  range  makes  it  even 
more  difficult  to  ensure  far-field  radiation  for  all 
plausible  types  of  antennas  and  receiver  positions. 

B.  MODELS  FOR  INDOOR  PROPAGATION 

To  make  things  more  challenging,  it  is  impossible  to 
design  an  "RF  friendly"  building,  free  of  multipath 
reflections  from  numerous  smooth  surfaces,  diffraction 
around  sharp  corners  or  scattering  from  walls  and  ceilings. 
In  addition,  the  RF  signal  is  subject  to  temporal 
variations  due  to  the  motion  of  people,  doors,  etc.  [13], 
[14],  [15],  [16]. 

Radio  wave  propagation  inside  metal-walled 
compartments  such  as  a  missile  room  of  a  military  ship  is 
dominated  by  strong  reflections  off  the  compartment  walls. 
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strength  at  2.4  GHz  can  vary 


In  such  an  environment ,  signal 
dramatically  over  a  distance  of  a  fev.*  centimeters.  Despite 
these  problems,  vv-e  are  still  capable  of  predicting  radio 
propagation  in  a  typical  room^  or  a  shipboard  compartment. 
To  this  effect  a  number  of  models  have  been  developed  and 
sem.i-empirical  formulas  for  indoor  signal  propagation 
prediction  developed  from;  propagation  studies.  Figure  4 
illustrates  the  results  of  a  study  of  a  signal  at  1.8  GHz, 
penetrating  a  vj-indov.?  of  a  topical  room;. 


Figure  4.  The  Penetration  of  a  Mobile  Phone  Signal  at  1.8 
GKz  Through  a  Windov;  in  a  Reinforced  Concrete  Wall  is 
Studied  by  Means  of  a  Tv/O-Dim-ensional  Model.  Right: 

Geom^etry .  Left:  Field  Levels  in  dB  Due  to  the  Line  Source 
in  the  Left  Upper  Corner.  (From;  [17]). 

Generally,  models  can  be  classified  as  either 
statistical,  based  on  em.pirical  data,  or  deterministic, 
valid  for  the  soecific  environmient 


under  consideration. 


Some  models  do  not  belong  in  either  category,  since  they 
are  a  combination  of  those  two.  These  are  referred  to  as 
hybrid  models. 

Statistical  models  [18]  and  [3],  are  based  on  fitting 
curves  or  analytical  expressions  to  observed  empirical 
data,  such  that  the  propagation  mechanisms  are  accounted 
for  in  an  average  sense.  Since  reality  seldom  corresponds 
to  relatively  simple  mathematical  models,  a  certain  range 
of  error  cannot  be  avoided.  The  accuracy  of  these  models 
depends  on  the  accuracy  of  the  measurements  taken,  as  well 
as  che  similarities  of  the  room  to  be  examined  with  the 
(usually  simpler)  room  in  which  the  results  were  collected. 
However,  the  statistical  models  are  simple  to  implement  and 
computationally  efficient. 

With  Deterministic  models,  (also  called  site  specific 
(SISP)  models  [19],  [20],  [21],  [22],  [23],  [24],  [25] 
specific  environments  are  characterized  individually 
(rooms,  buildings,  tunnels,  shipboard  compartments  etc. ) . 
In  this  way  the  model  accuracy  is  improved,  as  long  as 
there  is  detailed  database  available  for  the  geometry, 
materials,  location/motion  of  people  etc  for  the  indoor 
environment  under  consideration.  The  obvious  and  very 
important  drawback  in  this  case  is  that  even  if  the 
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geometry /media  database  is  available,  the  algorithms  for 
indoor  propagation  calculation  are  typically  very  complex 
and,  depending  on  the  accuracy  desired,  model  calculations 
can  be  extremely  time  consuming. 

In  general,  various  "ray  tracing"  algorithms  are 

employed  for  indoor  propagation  prediction.  Such 

algorithms  typically  implement  an  exhaustive  search  of  a 
"ray  tree"  as  follows:  first  and  foremost,  a  ray  tracing 
code  detects  if  there  is  a  LOS  path  between  the  receiver 
and  the  transmitter  and  determines  the  RF  energy  received 

via  the  direct  path.  Then,  depending  on  the  source  ray 

algorithm,  the  program  detects  if  an  object  intersection 
occurs.  If  such  is  the  case,  code  distinguishes  the 
phenomenon  as  scattering  or  diffraction  or  reflection  and 
computes  the  received  signal,  otherwise  process  stops  and  a 
new  ray  is  "launched."  Typically,  as  the  number  of  rays 

encountered  in  the  problem  increases  linearly,  the  number 
of  required  computation  increases  exponentially. 
Therefore,  for  given  computational  resources,  there  is  a 
limit  on  the  relative  size  of  the  indoor  geometry  'that  a 
particular  code  can  be  applied  to. 

Hybrid  models  [19],  [26]  can  be  a  successful 

combination  of  the  statistical  and  site-specific  models  by 
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trading  accuracy  for  computational  efficiency.  A  survey  of 
field  strength  prediction  models  is  presented  in  [27] . 

NEC-BSC  utilizes  ray-tracing  technique  based  on  the 
Uniform  Theory  of  Diffraction  (UTD)  and  the  geometry  and 
materials  data  for  the  specific  environment  to  be  analyzed. 
Therefore  it  can  be  classified  as  a  physical  model. 

C.  POTENTIAL  BENEFITS  FOR  THE  NAVY 

As  mentioned  previously,  the  majority  of  materials 
used  in  a  shipboard  compartment  are  metals  (good 
conductors)  that  can  be  approximated  as  Perfect  Electric 
Conductors  (PECs).  However,  the  geometry  of  a  shipboard 
compartment  (on  the  scale  of  the  operating  wavelength)  can 
be  very  complex  and  the  behavior  of  the  transmitted  signals 
can  be  difficult  to  predict.  Furthermore,  very  limited  and 
often  classified  research  and  documentation  is  available 
for  shipboard  wireless  propagation. 

In  1998,  a  bistatic  measurement  system  was  developed 
for  the  0.8  to  3  GHz  frequency  range  [28].  Experiments 
took  place  onboard  the  former  USS  SHADWELL  for  a  variety  of 
types  of  confined  rooms,  locations  of  transmitters  and 
receivers,  as  well  as  open  or  .closed  door  configurations. 
The  results  led  to  the  conclusion  that  the  expected 
elimination  of  RF  energy  due  to  the  ship's  metallic 
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structure  (bulkheads,  tanks,  overheads,  watertight  doors, 
etc.)  did  not  apply.  The  major  phenomenon  observed,  was  an 
approximate  10-dB  signal  attenuation  for  each  closed 
passageway  door  and/or  bulkhead  [28] . 

Let  us  now  examine  the  motives  that  urged  the  Navy  to 
sponsor  these  and  other  similar  measurements.  Over  the  last 
few  years,  there  is  a  growing  interest  in  placing  Wireless 
Local  Area  Networks  (WLANs)  in  ships  and  submarines.  Future 
generations  of  large  ships  are  being  designed  with  more  and 
more  electronic  systems  but  with  fewer  personnel. 
Currently,  it  is  generally  agreed  that  WLANs  are  more 
efficient  than  conventional  LANs.  First,  WLANs  can  provide 
services  almost  everywhere  on  the  ship,  without  the 
physical  limitations  that  cables  have.  WLANs  also  mitigate 
problems  due  to  physical  damage  to  wires  or  fiber  optic 
cables.  Finally,  a  sophisticated  WLAN  can  also  offer 
installation  flexibility,  reduced  cost  and  ability  to 
transfer  vast  amounts  of  sensor  data  related  to  condition, 
operation  and  damage  control  of  the  majority  of  ship's 
departments . 

D.  OPERATION,  IMPLEMENTATION  AND  SAFETY  ISSUES 

So  far,  the  most  common  and  convenient  application  has 
been  the  hand-held  walkie-talkies  that  crewmembers  use. 
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During  state  of  readiness,  communication  can  be  established 
from  various  locations  even  if  all  intermediate  watertight 
doors  are  closed.  However,  this  application  could  lead  to 
mild  to  severe  Electromagnetic  Interference  (EMI)  problems 
[28]  . 

For  all  potential  WLAN  applications,  EMI  is  the 
system  "bottleneck."  Therefore,  WLANs  have  to  be  designed 
using  the  best  combination  of  transmitter  power,  receiver 
performance  and  operating  frequency.  In  1985,  the  Federal 
Communication  Commission  (FCC)  authorized  license-free 
spread  spectrum  RF  communications  for  low  power  (less  than 
1  watt)  transmitters  in  the  910  MHz,  2.44  GHz,  and  5.8  GHz 
ISM  bands  [29]  .  This  enabled  a  rapid  growth  of  spread 
spectrum  WLAN  industry,  primarily  for  commercial 
applications.  Lately,  the  technology  is  being  applied  to 
military  use  as  well.  Finally,  since  the  transmitter  power 
levels  of  WLAN  systems  are  typically  less  than  1  Watt  and 
radiated  power  density  is  attenuated  at  least  20  dB  per 
decade  increase  in  range,  crewmembers  present  in  a  WLAN 
coverage  area  are  not  at  risk  to  high  RF  exposure. 
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III.NEC-BSC  SOFTWARE  OVERVIEW 


A.  INTRODUCTION  TO  NEC-BSC 

NEC-BCS,  developed  at  Ohio  State  University  under 
government  contract,  has  been  used  over  the  last  ten  years 
for  solving  EM  scattering  and  radiation  problems  [34] . 
More  specifically,  this  user-oriented  code  is  capable  of 
analyzing  near  and  far  field  radiation  and  scattering  in  an 
arbitrarily  shaped  room/ compartment  at  UHF  (300  MHz  to  3 
GHz)  and  above.  Rooms  and  compartments  can  be  modeled  by 
combining  a  number  of  basic  structural  shapes:  multi-sided 
flat  plates  and  finite  elliptic  curved  surfaces  like 
cylinders,  cone  frustums,  and  finite  composite  ellipsoids. 

Over  the  past  years,  NEC-BSC  has  been  successfully 
used  to  predict  scattering  from  superstructures  of  ships, 
wings  of  aircraft,  large  tanks  and  indoor  compartments  of 
space  stations,  simulated  using  perfectly  conducting 
materials.  The  newest  version  of  NEC-BSC  (4.2) 

incorporates  a  thin  dielectric  slab  option,  hence  the  code 
can  be  used  to  model  propagation  within  buildings  with  non- 
metallic  walls  [38] . 

B.  NEC-BSC  MODELS  AND  LIMITATIONS 

NEC-BSC  utilizes  advanced  formulations  based  on  the 
Geometric  Theory  of  Diffraction  (GTD)  and  Uniform  Theory  of 


Diffraction  (UTD)  [30],  [31]  to  solve  problems  involving 
electromagnetic  coupling,  scattering,  diffraction  and 
computation  of  near-  and  far-field  antenna  patterns. 

Geometrical  Theory  of  Diffraction  is  a  high-frequency 
technique  suitable  for  the  analysis  of  antennas  that  are 
large  in  terms  of  wavelength  (□)  and  located  in  a  complex 
environment.  Keller  [32],  [33]  originally  formulated  GTD 
as  an  extension  of  Geometrical  Optics  (GO)  .  GTD  proved  to 
be  very  reliable  in  calculating  the  non-zero  fields  in  a 
shadow  region.  This  was  accomplished  by  introducing  a  set 
of  diffracted  rays  analogous  to  the  reflected  and 
transmitted  rays  of  GO. 

In  NEC-BSC,  ray-tracing  algorithms  define  the 
parameters  of  the  field  incident  onto  various  obstacles. 
The  diffracted  rays  are  then  defined  by  UTD  for  all 
incident  rays,  and  added  to  the  geometrical  optics  terms 
for  the  observation  points  of  interest.  Since  most  of  the 
scattered  and  diffracted  rays  will  encounter  obstacles 
themselves,  a  rapidly  increasing  number  of  higher  order 
scattered  and  diffracted  rays  will  be  created.  In  order  to 
limit  the  number  of  computations  for  a  solution  with 
acceptable  accuracy,  a  sophisticated  code  is  required  to 
"trace  out"  the  various  possible  combinations  that  will  not 


26 


contribute  significantly  to  the  total  field  at  the  selected 
observation  points.  In  this  manner  only  the  components 
with  the  strongest  contributions  need  to  be  calculated, 
thus  increasing  the  code's  computational  efficiency. 

There  are,  however,  some  limitations  of  this  user- 
oriented  code,  which  basically  occur  due  to  the  very  nature 
of  the  algorithms  employed  [34]  .  Since  all  diffracted 
components  are  calculated  using  UTD,  each  plate  should  have 
edges  of  at  least  a  wavelength  long.  In  terms  of  curved 
surfaces  their  major  and  minor  radii  and  length  should  be 
at  least  a  wavelength  in  extent.  Furthermore,  for 
dielectric  (non-metallic)  obstacles,  the  source  has  to  be 
located  at  least  a  wavelength  away  from  the  surface  and  the 
incident  field  angle  of  incidence  should  not  be  too  close 
to  the  grazing  angle.  Also,  the  dielectric  layer  must  be 
thin,  at  most  half  of  the  operating  wavelength.  For 
example,  at  2.4  GHz  (12.5cm  wavelength)  dielectric 
thickness  must  be  less  than  6.25  cm.  Moreover,  antenna 
elements  have  to  be  at  least  half  a  wavelength  from  all 
edges.  In  any  case,  the  antenna  can  not  be  mounted  on 
curved  surfaces . 

Limitations  on  the  NEC-BSC  accuracy  also  stem  from  the 
high  order  term  elimination.  Eliminating  higher  order 
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diffracted  rays  improves  computational  efficiency,  but  it 
also  makes  the  code  unsuitable  for  certain  types  of  EM 
calculations.  For  example,  higher  order  diffracted  rays 
are  important  for  accurately  predicting  the  radar  cross 
section  of  a  complex  structure.  ,  Furthermore,  flat  plates 
and  simple  curved  surface  models  are  not  always  adequate 
for  modeling  complex  environments  found  within  buildings  or 
shipboard  compartments.  Due  to  all  these  limitations,  NEC- 
BSC  program  proved  to  give  accurate  results  for 
approximately  the  first  25  dB  of  the  results  dynamic  range. 

.As  a  whole,  the  present  NEC-BSC  execute  module  will 
accommodate  models  with  the  following  parameter  limits: 


Pattern  points 

1801 

Defined  pattern  points 

= 

181 

Defined  frequencies 

= 

181 

Plates 

= 

200 

Plate  edges 

= 

30 

Curved  entities 

= 

200 

Curved  sections 

= 

10 

Material  types 

= 

25 

Material  layers 

= 

20 

Sources 

— 

300 
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•  Receivers 


300 


Antenna  data  points 

=  91 

Wires 

=  200 

C.  MODELING  THE  STRUCTURE 

Most  of  the  structures  of  practical  interest  can  be 
successfully  approximated  by  a  combination  of  multi-sided 
flat  plates  and  some  finite  elliptic  surfaces.  These  can 
consist  of  either  good  conductors  (even  PECs)  or 
dielectrics.  A  group  of  flat  plates,  each  having  a  certain 
number  of  edges,  can  model  most  of  the  objects  with  slowly 
varying  radii  of  curvature.  In  either  case  there  is  a 
certain  procedure  that  one  has  to  follow  in  order  to 
construct  an  arbitrary  object. 

If  a  plate  has  to  be  modeled,  one  first  has  to  specify 
the  corners  in  a  three  dimensional  Cartesian  coordinate 
system.  The  only  condition  is  that  all  these  corners 
should  form  a  flat  plate.  If  the  points  are  difficult  to 
specify,  the  user  can  simply  construct  the  plate  in  two 
dimensions  and  then  rotate  or  translate  it  into  its  desired 
location  in  three  dimensions  using  the  appropriate  commands 
provided  by  the  software. 

A  plate  can  be  attached  to  another  plate  in  NEC-BSC  in 
two  different  ways.  First,  two  plates  can  share  a  pair  of 
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corners,  forming  a  wedge.  The  angle  of  the  wedge  is 
defined  by  the  angle  of  the  two  plates,  based  on  the 

direction  of  antenna's  illumination.  Alternately,  one 
plate  can  be  mounted  slightly  within  the  volume  of  the 
other.  By  piecing  together  a  number  of  such  isolated 
plates,  a  larger,  more  complex  general  structure  can  be 
approximated  with  sufficient  accuracy. 

NEC-BSC  also  has  the  capability  to  implement 
electrically  thin  dielectric  layers  in  order  to  model 
structures  like  non-metallic  doors,  walls,  glass  windows, 
furniture,  etc.  For  flat  dielectric  plates,  the  user  has 
to  specify  the  thickness  (not  exceeding  one  half  of  the 

operating  wavelength) ,  relative  permitivity,  relative 
permeability,  and  the  loss  tangents  for  the  material. 

If  a  surface  with  a  non-trivial  curvature  has  to  be 
modeled,  one  can  use  elliptic  cylinders,  cone  frustums  and 
composite  ellipsoids.  In  NEC-BSC  all  curved  surfaces,  are 
defined  by  their  center  coordinates  and  their  major  and 
minor  axes.  In  case  of  an  elliptic  cylinder,  the  position 
of  the  first  and  the  last  cap  is  given  with  respect  to  the 

cylinder  axis.  The  cone  frustum  is  characterized  by  the 

position  of  its  rim  junctions.  Finally,  the  ellipsoid  is 
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constructed  by  setting  different  values  of  radii  curvature 


along  its  maximum  and  minimum  z-axis  positions. 

D.  MAJOR  COMPONENTS  OF  NEC-BSC  MAIN  PROGRAM 

The  NEC-BSC  code  is  constructed  in  a  systematic 
manner,  meaning  that  all  procedures  are  categorized  and 
executed  in  modular  sections.  Figure  5  shows  the  major 
components  of  the  NEC-BSC  program. 


riiiitialize  Data] 


Figure  5.  Block  Diagram  of  NEC-BSC  Main  Program  (from 
[34]  )  . 
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In  the  input  part  the  problem  geometry  is  specified. 
The  second  section  prepares  all  necessary  parameters  for 
later  use.  In  the  following  section  the  fields  resulting 
from  individual  propagation  mechanisms  are  calculated  and 
added  to  form  the  solution.  Finally,  all  the  outputs  are 
provided  to  the  user,  either  in  the  form  of  numerical 
results  or  as  ray  path  visualization. 

E.  TABLE  OF  COMMANDS 

A  listing  of  NEC-BSC  commands  is  provided  in  the 
preliminary  draft  of  the  NEC-BSC  user's  manual.  (The 
listing  of  commands  is  still  incomplete  at  the  time  of  this 
writing)  .  The  commands  can  be  grouped  into  six  broad 

categories:  the  geometry  control  commands,  the  structure 

definition  commands,  the  antenna  definition  commands,  the 

pattern  control  commands,  the  output  control  commands,  and 

the  program  control  commands .  In  creating  a  NEC— BSC  input 
file  the  order  of  groups  should  be  as  given  above,  but 
there  is  no  required  order  for  the  commands  within  the  same 
group.  The  default  unit  is  meters  for  all  lengths  and  the 
default  frequency  is  299.9  MHz.  All  angles  are  specified 
in  degrees.  The  following  examples  illustrate  most  of 
these  points  and  validate  a  sample  NEC-BSC  solution. 
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IV. 


DEVELOPMENT  OF  BASIC  EXAMPLES 


In  this  chapter,  NEC-BSC  results  for  several  sample 
problems  are  presented.  All  pertinent  Matlab,  NEC-BSC,  and 
Mathcad  files  are  included  in  Appendices  A  through  C, 
respectively . 

The  sample  problems  have  been  selected  to  illustrate 
various  features  of  NEC-BSC.  Each  problem  has  been 

selected  in  order  to  show  specific  code  capabilities  and 
convenient  options,  as  well  as  some  unresolved  problems 
still  present  in  the  most  recent  version  of  the  code.  The 
collection  of  results  provides  an  insight  into  the 
tradeoffs  of  the  NEC-BSC  use  for  indoor  propagation 
prediction . 

A.  COMPARISON  OF  NEC-BSC  AND  "PATCH"  CODE  RCS  PREDICTIONS 

The  problem  of  EM  scattering  by  a  smooth,  perfectly 

conducting  convex  surface,  is  investigated  in  this  section. 

The  scattered  field  results  obtained  from  the  "PATCH"  code 

[9]  will  be  compared  to  the  scattered  field  calculated  by 

the  NEC-BSC  for  the  same  geometry.  The  first  example 

involves  an  infinitely  thin  electric  dipole,  shown  as  a  red 

dot  in  Figure  6,  located  parallel  and  relatively  close  to  a 

metallic  circular  cylinder.  The  cylinder  axis  coincides 

with  the  z-axis  of  the  coordinate  system.  The  operating 
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frequency  is  2.54  GHz,  which  corresponds  to  the  operating 
wavelength  of  11.81  cm.  The  obstacle  is  a  metallic 

cylinder  1  meter  in  length  and  having  a  17.5  cm  radius. 
The  antenna  is  1  meter  from  the  surface  of  the  cylinder. 
The  observation  points  of  interest  form  a  horizontal  circle 
of  360  points,  5  meters  away  from  the  cylinder  axis,  as 
shown  in  the  Figure  6. 


Figure  6.  Geometry  of  the  PEC  Convex  Surface  and  the 
Dipole. 

The  validity  of  NEC-BSC  results  is  shown  by  comparing 
the  NEC-BSC  results  with  the  solution  obtained  for  the  same 
problem  by  using  the  "PATCH"  code.  The  results  from  the 
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two  codes  are  shown  combined  in  a  single  graph  in  Figure  7 . 


Electric  field  measured  by  patch  (blue)  and  NEC-BSC 


Figure  7.  Comparison  of  Electric  Fields  Calculated  by  NEC- 
BSC  and  PATCH  Codes. 

It  can  be  clearly  stated  that  the  match  is  not 
perfect.  Although  the  number  of  minima  and  maxima  are  the 
same  and  the  "shape"  of  the  field  variation  with  angle  is 
almost  identical,  there  is  approximately  1.7  dB  difference 
between  the  values  at  each  point  of  the  graph.  Moreover, 
in  the  very  deep  shadov/  area  exactly  behind  the  cylinder, 
there  seems  to  be  som^e  differences  as  well,  although  at  a 
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relatively  low  level.  Despite  the  small  differences,  there 
is  good  enough  agreement  between  the  two  to  validate  the 
NEC-BSC  code  results. 

NEC-BSC  has  been  also  extensively  validated  by  its 
author.  User's  manual  [34]  contains  several  examples 
involving  metallic  cylinders,  where  patterns  were  taken  in 
many  different  cuts,  for  many  source  positions  and  for 
different  polarization.  For  these  examples,  comparison  of 
NEC-BSC  predictions  with  the  measured  results  shows  a  very 
good  agreement. 

B.  SIMPLE  BOOM  WITH  METALLIC  WALLS 

This  example  considers  a  simple  room  with  PEC  walls 
and  the  following  dimensions: 

•  Length  of  six  meters 

•  Width  of  twelve  meters 

•  Height  of  twelve  meters. 

The  room  walls  are  modeled  using  six  PEC  plates. 
Inside,  there  is  a  half-wavelength  antenna  located  at 
(3,1,11)  with  respect  to  the  origin,  as  shown  in  Figure  8 
(the  antenna  is  depicted  as  a  red  dot)  .  The  operating 
frequency  is  2.4  GHz,  which  corresponds  to  a  wavelength  of 
12.5  cm.  A  large  metallic  cylinder,  representative  of  a 
missile  tube,  is  positioned  in  the  room  such  that  the 
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antenna  illuminates  only  about  half  of  the  surface  of  the 
cylinder.  Our  objective,  in  this  specific  example,  is  to 
visualize  the  propagation  paths,  especially  the  non-LOS 
ones.  The  observation  points  are  thus  positioned  in  a 
semicircle  and  appear  as  gray  dots  in  the  figure  below. 


Figure  8.  Room  With  PEC  Walls  and  a  Single  PEC  Cylinder. 

Running  the  NEC-BSC  program  for  even  this  relatively 
simple  geometry  results  in  a  plethora  of  rays.  Hence,  even 
if  we  rotate  or  zoom  into  the  figure,  we  will  not  get  a 
clear  picture  of  the  individual  propagation  mechanisms 
involved. 
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NEC-BSC,  however,,  provides  the  user  with  two  very 
convenient  options.  The  first  is  that  the  user  can  adjust 
the  dynamic  range  of  the  results  (in  the  ray  menu)  such 
that  only  the  strongest  rays  are  displayed  (e.g.,  for  LOS 
paths  and  single  reflections  only) .  For  instance,  in 
Figure  9  only  the  rays  within  1  dB  of  the  strongest  ray  are 
shown.  Second,  we  can  select  a  specific  ray  and  transfer 
it  to  a  figure,  which  only  shows  the  geometry  of  the  room 
(the  "drag  and  drop"  method) .  Therefore  a  "hand-picked" 
number  of  rays  can  be  superimposed  to  the  room  geometry. 
To  illustrate  this  for  the  room  with  the  single  cylinder 
example,  only  the  diffracted  rays  to  some  certain 
observation  points  are  shown  in  Figure  10. 


Figure  9.  Propagation  Paths  for  Signals  Within  1  dB  of  the 
Strongest  Signal  Path. 
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Figure  10.  Sample  Diffracted  Rays,  "Isolated"  from  Figure 
7. 

C.  COMPLEX  ROCM  WITH  METALLIC  WALLS 

Let  us  now  consider  a  more  complex  room  with  PEC  walls 
and  the  same  dimensions  as  the  room  from  the  preceding 
section: 

♦  Length  of  six  meters 

♦  Width  of  twelve  meters 

♦  Height  of  twelve  meters. 

Six  PEC  plates  form  the  room.  Inside,  there  is  a 
half-wavelength  antenna  located  at  (3,3,1)  with  respect  to 
the  origin,  depicted  as  the  red  dot.  The  operating 
frequency  is  2.45  GHz,  which  corresponds  to  the  wavelength 
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of  12.245  cm.  Two  large,  vertical  metallic  cylinders, 
representing  missile  tubes,  are  placed  in  the  center  of  the 
room.  In  addition,  two  large  metallic  tanks  are  located 
along  the  front  and  the  rear  wall  of  the  room,  as  shown  in 
Figure  11.  The  objective,  in  this  specific  example,  is  to 
obtain  relative  signal  strength  "map"  of  the  room, 
revealing  the  areas  where  the  received  signal  is  strong, 
weak  or  virtually  non-existent. 


Figure  11.  Room  With  PEC  Wall,  Two  PEC  Cylinders  and  Two 
PEC  Storage  Tanks 


NSC-BSC  does  not  have  the  3D  "relief  map"  plot 
capability  to  represent  the  calculated  fields.  A  program 
has  been  v/ritten  in  Mathcad  to  create  the  relative  signal 
strength  "relief"  (3D)  map.  More  specifically,  using  the 
Notepad  application,  we  convert  the  data  file  that  NEC-3SC 
creates  for  the  eleccric  field  at  the  observation  points  to 
an  ASCII  text  file  ( * .  txt  extension).  The  .ASCII  data  file 
is  next  imported  into  a  Mathcad  script,  and  the  electric 
field  m;agnitude  is  plotted  in  dB,  as  shown  in  the  ne.xt 
Figure . 

The  observation  points  are  all  in  the  same  plane  so 
that  the  tv/o  horizontal  axes  are  the  coordinates  in  the 
observation  plane  and  the  color-coded  vertical  axis  (the 
"height")  indicates  the  electric  field  strength. 


Figure  12.  Signal  Coverage  of  the  Room,  in  Which  Red 
Corresponds  to  0  dB  and  Blue  to  -25  d3. 
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Figure  13.  Top  View  of  the  Signal  Coverage  of  the  Room,  in 
Which  Red  Corresponds  to  0  dB  and  Blue  to  -25  dB. 

The  red  color  identifies  the  areas  of  the  highest 
signal  strength  v^hile  blue  identifies  the  areas  of  weakest 
signal.  It  is  important  to  remember  that  although  the  NEC- 
BSC  will  report  signal  strengths  with  very  high  dynamic 
ranges,  the  NEC-BSC  user's  manual  [34]  cautions  that  only 
about  30  dB  of  the  dynamic  range  (30  dB  down  from  the 
highest  signal  level)  should  be  considered  reliable. 
Therefore,  one  can  opt  to  place  a  threshold  at  30,  dB  below 
the  largest  signal  strength  and  thus  establish  a  scale  with 
red  color  indicating  the  largest  signal  strength  and  blue 
color  indicating  signal  strengths  at  30  dB  lower  (than  the 
maximum) . 
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As  expected,  the  signal  strength  for  obseirvation 
points  v;ithin  PEC,  objects  (the  two  cylinders  and  the  tanks) 
are  represented  by  blue  color,  while  the  area  close  to  the 
antenna  has  a  large  number  of  red  "spikes"  corresponding  to 
the  standing  wave  peaks.  The  pronounced  standing  waves  are 
a  consequence  of  strong  reflection  off  the  PEC  v/alls. 
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D.  OPENINGS 

In  this  section  the  performance  of  NEC-BSC  is 
evaluated  for  openings  such  as  doors,  hatches,  and  windo'ws. 
The  opening  will  be  constructed  according  to  the  method 
recommended  in  the  NEC-BSC  user's  manual.  An  electrically 
large  opening,  showing  in  Figure  14,  'will  be  constructed 
using  a  number  PEC  plates  attached  together  (by  the  virtue 
of  shared  corner  pairs) .  These  individual  plates  form  a 
larger  flat  plate  with  the  opening  in  the  middle.  A  half- 
wavelength  antenna  is  positioned  on  one  side  of  the  opening 
and  a  set  of  observation  points  is  defined  on  the  opposite 
side  as  shov.mi  in  Figure  15.  In  addition  to  the  direct  rays 
passing  through  the  opening  we  expect  to  see  the  rays 
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Observation 


Figure  15-  The  Flat  PEC  Plate  With  the  Opening  in  the 
Middle,  as  Constructed  by  NEC-BSC. 

The  NEC-BSC  results  for  the  PEC  plate  with  the  opening 
show  that,  in  addition  to  the  rays  diffracted  from  the 
edges  of  the  opening,  there  are  non-physical  rays 
diffracted  by  the  seams  of  the  individual  plates,  as  shown 
in  Figures  16  and  17. 
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Figure  16.  Sample  Propagation  Paths  for  the  PEC  Plate  With 
the  Opening. 


Figure  17.  Non-Physical  Rays  "Diffracted"  by  the  Plate 
"Seams 

Unfortunately,  this  problem  was  not  resolved  either 

by  constructing  the  large  plate  using  overlapping  (by  0-1 

of  the  wavelength)  individual  plates,  or  by  using  the 

multiple  date  (PM)  command-  We  have  therefore  concluded 
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that  the  most  likely  cause  is  an  unreported  {and  perhaps  an 
undececced)  "bug''  in  uhe  code  that:  could  afreet  the 
accuracy  of  the  code  in  certain  cases  (for  example  for 


room;- to- room; 

propagation 

sim.ulation)  .  A  v;ay  to 

circumvent 

the  "leakage 

through  the 

sea.m;S"  problem,  has  been 

discovered 

by  trial  and  error. 

In  this  approach  one  starts  with  a  single,  solid  PEC 
plate.  Next,  a  second  plate  is  defined  such  that  it 
coincides  with  the  desired  opening  in  the  PEC  plate.  This 
second  plate  is  assigned  electrical  properties 
(permittivity  and  permeability)  of  free  space.  The  NEC-BSC 
apparently  treats  the  second  plate  as  an  "opening,"  v;ith 
direct  and  diffracted  rays,  but  'without  the  non-physical 
rays  diffracted  from;  the  "seam;S"  (because  there  are  no 
sea.m.s).  Figure  18  shov.'S  the  LOS  and  the  diffracted  paths 
for  the  "seamdess"  PEC  plate  with  the  opening.  The 
di fractions  from;  the  edges  of  the  opening  are  evident. 
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igure  18.  Propagation  Paths  for  the  "Seamless"  Plate  With 
he  "Opening"  Formed  by  a  "Free  Space  Plate". 
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V. 


TYPICAL  MISSILE  ROOM 


A.  GENERAL  DERCRIPTION 

As  mentioned,  providing  good  signal  coverage  for 
wireless  communications  on  a  military  ship  can  be  a  very 
challenging  task.  One  of  the  objectives  is  to  provide  the 
best  coverage  with  the  smallest  number  of  "access  point" 
stations  ("base"  stations).  To  accomplish  this,  computer 
simulations  may  be  used  for  the  particular  shipboard 
compartments  to  determine  the  signal  coverage  for  various 
combinations  of  the  numbers  and  locations  of  antennas  (with 
multiple  antennas  it  is  assumed  that  they  transmit 
identical  signals;  that  is,  the  "multicast"  operation  is 
assumed) .  The  requirement  to  minimize  the  computational 
time  must  be  balanced  against  the  desired  accuracy  of  the 
solution.  Therefore,  the  optimum  level  of  geometric  detail 
in  the  compartment  models  results  from  a  tradeoff  between 
the  desired  accuracy  and  the  computational  resources 
available  (CPU  time  and  memory) . 

Simulation  of  a  "typical"  shipboard  missile  room  is 
presented  in  this  chapter.  Although  the  geometry  of  missile 
room  does  not  represent  the  missile  compartment  of  any 
particular  ship,  it  is  typical  of  missile  compartments  for 
submarines . 


51 


The  following  software  versions  were  used  for  the 
simulations : 

•  NEC-BSC  graphical  user  interface  (GUI)  called 
Workbench,  version  4.1.37  and 

•  NEC-BSC  executable  version  4.2-06. 

These  were  the  most  recent  releases  by  Prof.  Ron  Marhefka, 
Ohio  State  University,  at  the  time  of  completion  of  this 
thesis. 

In  shipboard  environments  the  majority  of  materials 
involved  are  metallic.  Therefore,  although  lossy 

dielectric  materials  -  are  also  present  in  shipboard 
compartments  (people  for  example) ,  propagation  within 
shipboard  compartments  will  be  dominated  by  reflections  and 
diffraction  of  metallic  surfaces.  Furthermore,  the  metals 
involved  may  be  approximated  as  PEC. 

The  large,  parallelogram- shaped  compartment  shown  in 
Figure  9  has  the  following  dimensions: 

•  Length  of  twenty  meters 

•  Width  of  ten  meters 

•  Height  of  ten  meters. 

Six  PEC  plates  form  the  compartment.  Located  at  the  first 
floor  is  a  half -wave length  antenna,  depicted  as  the  red  dot 
in  the  figures.  The  operating  frequency  is  2.4  GHz,  which 
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corresponds  to  the  wavelength  of  12.5  cm.  Six  large, 
vertical  metallic  cylinders,  representing  missile  tubes, 
are  placed  in  the  room.  Each  tube  is  2  meters  in  diameter 
and  10  meters  high. 

Furthermore,  there  are  three  large  tanks,  which  cover 
the  walls  parallel  to  the  x-axis.  The  front  and  the  rear 
tanks'  height  are  2.5  and  0.5  meters  respectively,  while 
the  middle  one's  width  is  1  meter  and  the  height  is  5 
meters  (same  as  the  height  of  the  first  floor's  ceiling). 
In  addition  to  these  six  tanks,  there  is  one  large  tank  of 

2.5  meters  height  along  the  wall  parallel  to  y-axis  on  the 
rear  side  of  the  room. 

One  hatch  is  located  in  the  middle  of  the  front 
bulkhead.  Since,  hatches  do  not  start  at  the  floor  level 
for  the  safety  reasons  (water  containment) ,  the  hatch  is 
positioned  30  centimeters  from  the  floor  and  its  height  is 

1.5  meters.  Two  similar  hatches  are  located  on  the  first 
and  the  second  overheads . 

The  hatches  can  be  modeled  as  either  closed  or  opened. 
Therefore,  propagation  from  one  compartment  section  to  the 
others  through  open  hatches  can  be  simulated.  Finally,  a 
nine-scep  stairway,  complete  with  metallic  handrails, 
connects  the  first  and  the  second  floor.  The  geometry  of 
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the  NEC-BSC  missile  compartment  model  is  shown  in  Figure 
19.  The  top  view  of  the  geometry  of  the  compartment  is 
provided  in  Figure  20. 


Figure  19.  General  View  of  a  Typical  Missile  Shipboard 
Compartment . 
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B.  NEC-BSC  MODEL  SPECIFICS 

The  NEC-BSC  model  described  in  the  preceding  section 
met  all  the  requirements  with  respect  to  the  limitations  of 
the  code,  stated  in  Chapter  III.  In  addition,  all  plates 
were  modeled  in  a  way  that  the  code  could  infer  which  side 
of  the  plate  was  illuminated  by  the  antenna.  This  was 
accomplished  by  defining  all  plates  according  to  the  right- 
hand  rule  described  in  Chapter  V  of  the  NEC-BSC  manual .  To 
eliminate  NEC-BSC  run-time  "warnings"  concerning  the  wall- 
attached  plates  of  the  rear  tank,  we  had  to  define  them  as 
having  corners  not  exactly  identical,  but  positioned  a 
small  amount  (0.5A)  through  the  surfaces  to  which  '  there 
were  being  attached. 

In  the  cases  where  we  had  to  place  identical  objects 

in  different  locations,  the  NEC-BSC  "Relative  Rotate- 

Translate"  (RR)  command  was  used.  This  command  enables  the 

user  to  translate  and/or  rotate  the  coordinate  system  used 

to  define  the  input  data  in  order  to  simplify  the 

specification  of  the  set  of  missile  tubes  and  the  numerous 

« 

stairs.  Unlike  the  RT  command,  which  always  referred  back 
to  the  reference  coordinate  system,  the  RR  command 
accumulated  the  translations  and/or  rotations  necessary  for 
the  desired  positions.  Therefore,  at  the  end  of  the 
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procedure  it  was  necessary  to  refer  to  the  global 
coordinates  again,  by  placing  and  then  commenting  out  a 
test  object  at  the  origin  that  verified  the  accuracy  of  the 
procedure . 

Finally,  this  NEC-BSC  model  did  not  account  for  the 
presence  of  the  people.  Although  presence  and/or  motion  of 
people  can  cause  up  to  10  dB  attenuation  of  the  average 
signal  strength  [35] ,  the  trend  towards  reduced  ship 
manning  [36]  through  increased  automation  and  sensor 
monitoring  via  ship-wide  networks,  means  that  the  presence 
of  a  large  number  of  personnel  in  a  missile  compartment  is 
unlikely . 

C.  SIMULATION  RESULTS 

1 .  Antenna  Located  in  the  Middle  of  the  Room 

In  this  case,  both  antenna  and  observation  area  are 
located  on  the  first  floor.  The  source  is  located  in  the 
middle  of  the  compartment,  2  meters  above  the  floor  level, 
while  the  field  observation  points  form  a  horizontal 
rectangle  of  400  by  200  points  (for  a  total  of  80,000 
observation  points)  at  the  same  height  as  the  antenna. 
Recall  that  the  operating  wavelength  is  12.5  cm,  while  the 
X  and  y  dimensions  of  the  compartment  are  2  0  by  10  meters 
respectively.  Therefore,  the  minimum  spatial  sampling  that 
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allows  us  to  capture  wave  minima  and  maxima  requires  a 
sampling  step  of  one  half  of  the  operating  wavelength, 
which  corresponds  to  at  least  320*160  =  51,200  observation 
points.  Increasing  the  density  of  the  observation  points 
results  in  a  more  accurate  "map"  of  the  field,  but  at  the 
expense  of  increased  computational  time.  It  was  decided 
that  80,000  points  represent  a  good  compromise  between  the 
detail  level  and  the  computational  time  required  for  the 
solution . 
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Antenna’s 

location 


Figure  20.  Top  View  of  the  Compartment. 

Despite  the  reasonable  computing  resources  used  (a  PC 
with  a  single  Pentium  III  600  MHz  processor  and  768 
Megabytes  of  RAM) ,  the  simulation  run  time  was  still 
approximately  seven  hours  for  the  80,000  observation 
points . 
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By  extracting  the  numerical  results  from.  NEC-BSC  and 
im.porting  them  to  Mathcad,  a  "miap"  of  the  signal  strength 
within  the  compartment  was  obtained,  shown  in  Figure  21. 
It  should  also  be  noted  that  the  threshold  for  the  signal 
strength  at  the  observation  points  was  set  at  40  dB  below 
the  peak  observed  signal  strength.  As  stated  before,  the 
dynamic  range  for  which  the  NEC-BSC  has  good  accuracy  is 
approximately  25  d3.  Nevertheless,  a  larger  dynamic  range 
has  been  selected,  in  order  for  the  shadow  areas  and  the 
obstacles  to  be  clearly  visible  in  the  signal  strength 
"map." 


Figure  21.  RF  Energy  Propagation  Inside  the  Room  (2  Meters 
Above  the  Floor) . 
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2 .  Results  for  the  Antenna  in  an  Arbitrary  Location 

Again,  both  antenna  and  observation  points  are  located 
on  the  first  floor.  The  antenna,  however,  is  now  placed  at 
an  "arbitrary"  location  (no  longer  in  the  middle  of  the 
compartment)  1  meter  above  the  floor,  while  the  observation 
points  form  a  horizontal  rectangle  of  400  by  200  points 
(for  a  total  of  80,000  observation  points)  and  1.5  meters 
above  the  floor.  The  top  view  of  the  compartment  is 
provided  in  Figure  22.  The  arrow  indicates  the  antenna 
location . 


60 


Antenna’s 

location 


I - 1 

Figure  22.  Top  View  of  the  Missile  Compartment  With  the 
Antenna  in  an  Arbitrary  Location. 

The  "maps"  of  the  signal  strength  are  shown  in  Figures 
23  and  24,  for  the  dynamic  ranges  of  33  and  22  dB 
respectively.  In  the  33-dB  dynamic  range  "map"  we  notice 
the  apparent  signal  "leakage"  into  the  missile  tubes  and 
the  tanks.  Since  the  tubes  and  the  tanks  are  non- 
penetrable  (modeled  as  PEC)  these  are  clearly  consequence 
of  the  code's  limited  dynamic  range.  Reducing  the  dynamic 
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range  to  22  dB  eliminates  most  of  the  "non-physical"  rays, 
as  evident  in  Figure  24. 


Figure  23.  Signal  Strength  "Map"  for  the  Missile  Room  Using 
a  Dynamic  Range  of  33  dB. 


Figure  24 .  Signal  Strength  "Map"  for  the  Missile  Room  Using 
a  Dynamic  Range  of  22  dB. 
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3. 


Results  for  two  Antennas  in  the  Missile  Room 


Next,  the  results  are  next  presented  for  a  "typical" 
wireless  LAN  installation  where  a  number  of  antennas  may  be 
used  to  increase  the  uniformity  of  the  signal  coverage 
within  the  user  space.  In  particular,  it  is  assumed  that 
the  two  antennas  transmit  identical  signals,  a  situation 
commonly  referred  to  as  "multicasting."  The  two  antennas 
are  centrally  located,  close  to  the  ceiling,  at  the  height 
of  4  m  above  the  floor,  and  8  m  apart  from  each  other,  as 
shown  in  Figure  25.  The  80,000  observation  points  again 
form  a  horizontal  rectangle  (400  by  200  points)  at  1  meter 
above  the  floor  (approximately  waist  level) . 
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Position  of  the  two  antennas 


I - 1 

Figure  25.  Top  View  of  the  Compartment  With  the  Two 
Antennas . 

By  extracting  the  NEC-BSC  solution  (numerical  results 
for  the  field  strength)  and  importing  it  into  Mathcad,  we 
obtained  the  "map"  of  the  signal  strength  within  the 
compartment,  shown  in  Figure  26.  The  dynamic  range  for 
this  "map"  was  limited  to  25  dB.  As  expected,  the  "map" 
shows  that  a  much  more  uniform  coverage  is  obtained  with 
two  antennas.  There  is  even  coverage  behind  the 

cylinders,  due  to  reflected  and  "creeping"  wave 
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contributions . 


The  computer  running  time  for  this 


simulation  was  approximately  12  hours. 


Figure  26-  Signal  Propagation  Inside  the  Missile  Room,  4 
Meters  Above  the  Floor  and  Using  a  Dynamic  Range  of  25  dB. 

4 .  Floor-to-Floor  Propagation 

The  simulations  have  confirmed  the  20  to  30  dB  dynamic 
range  limitation  of  for  NEC-BSC  stated  in  the  NEC-BSC 
User''s  Manual.  The  limited  dynamic  range  suggests  that 
NEC-BSC  will  not  be  well  suited  for  applications  where  the 
desired  signal  levels  could  be  very  low,  such  as  in  the 
case  of  floor-to-floor  propagation.  Nevertheless, 

simulations  were  perfomed  to  test  the  NEC-BSC  performance 
for  room-to-room.  and  floor-to-floor  propagation.  Presented 
here  are  the  NEC-BSC  results  for  floor-to-floor  propagation 
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signal  strength  for  the  floor-to-floor  propagation,  shown 
in  Figure  27.  Since  floor-to-floor  propagation  incurs  a 
considerable  loss  (no  LOS  path)  we  had  to  use  75-dB  dynamic 
range  for  the  "map."  This  is  clearly  beyond  the  stated 
reliable  dynamic  range  of  the  NEC-BSC  and  is  manifested  by 
what  appears  to  be  the  signal  "leakage"  at  the  "seams" 
between  the  plates  forming  the  floor/ceiling  and  the 
"seams"  between  the  tubes  and  the  floor/ceiling,  evident  in 
Figure  28.  Nevertheless,  the  cods  properly  predicts  the 
highest  signal  strength  for  the  vicinity  of  the  stairwell 
opening.  The  computer  run  time  for  this  simulation  was 
approximately  12  hours - 


Figure  28.  Signal  Propagation  Inside  the  Second  Floor, 
Caused  by  an  Antenna  Located  in  the  First  one  and  Using  a 
Dynamic  Range  of  75  dB. 
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However,  this  simulation  does  give  us  some  insight 
about  the  diffraction  and  reflection  mechanisms  that  take 
place  in  communication  between  floors.  These  phenomena, 
which  dominate  in  room-to-room  cases,  can  be  generally  seen 
in  Figure  28.  The  run  time  for  this  program  was 
approximately  12  hours. 

D.  PROPAGATION  STATISTICS  FOR  THE  MISSILE  COMPARTMENT 

1 .  Theoretical  Background 

Thus  far,  we  have  presented  the  signal  strength  "maps" 
for  particular  compartments  using  the  numerical  output 
results  of  the  NEC-BSC  and  Mathcad.  In  this  section  the 
objective  is  to  extract  the  statistics  of  the  excess  path 
gain/loss  (relative  to  free  space  path  loss)  for  the  entire 
floor  of  a  particular  compartment  and  for  different  antenna 
placements.  The  result  will  be  a  random  variable  with  a 
certain  probability  density  function  (pdf)  expressing  the 
likelihood  of  a  particular  excess  path  gain/loss. 

The  path  loss  for  free  space  can  be  easily  calculated 
for  any  distance  between  the  transmitter  and  an  observation 
point.  Adding  a  random  variable  with  a  known  pdf  to  the 
deterministic  free  space  loss  allows  us  to  calculate  the 
mean  value,  dispersion,  etc.  for  the  compartment  path  loss 
random  variable.  This  applies,  in  statistical  sense,  to 
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any  two  points  (the  transmitter  and  the  observation  point) 
inside  the  particular  compartment,  regardless  of  their 
exact  location  as  long  as  the  distance  between  them  is  the 
same. 

For  each  transmit  antenna  location,  two  NEC-BSC  input 
files  were  created:  one  for  the  compartment  and  the  other 
for  the  free  space  of  exactly  the  same  shape  and  volume  as 
the  compartment.  All  other  parameters  (antenna  type  and 
position,  radiated  frequency,  observation  points  placement) 
were  kept  identical.  The  two  output  ASCII  files,  resulting 
from  the  execution  of  the  input  files,  were  formatted  using 
the  Windows  WordPad  editor  and  imported  into  Mathcad.  Two 
column  matrices  were  created  in  Mathcad  for  each  input 
file,  with  the  first  column  containing  the  relative 
distance  between  the  transmit  antenna  and  the  observation 
point  and  the  second  column  containing  the  corresponding  E- 
field  strength. 

The  matrices  are  then  rearranged  in  the  following 
manner.  First,  the  distances  (from  minimum  to  maximum)  are 
discretized  into  distance  "bins,"  each  bin  50  cm  (4 
wavelengths)  wide.  Next,  all  values  within  the  same 
distance  bins  are  placed  on  the  same  row  of  a  new  "matrix" 
(this  "matrix"  has  a  variable  number  of  columns) .  The 
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distance  "bins"  geometrically  correspond  to  concentric 
spherical  shells  of  50  cm  thickness,  all  centered  at  the 
transmit  antenna  location.  The  number  of  observation 
points  within  each  "shell"  varies  with  the  shell  distance 
from  the  "origin"  (the  transmit  antenna  location) . 

The  next  step  is  to  form  the  ratios  of  the  compartment 
and  "free  space"  fields  at  the  same  observation  points. 
Expressing  this  ratio  in  dB  creates  a  new  variable  that  is 
referred  to  as  the  "excess  compartment  gain"  R.  Positive 
values  of  the  excess  compartment  gain  indicate  that  the 
signal  in  the  compartment  exceeds  the  free  space  signal  at 
the  same  location,  while  the  negative  values  of  the  excess 
compartment  gain  indicate  that  the  signal  in  the 
compartment  is  weaker  than  the  free  space  signal  at  the 
same  location. 


R  =  20  log 


^free 


(5.1) 


In  Equation  5.1,  R  denotes  the  "excess  room  gain"  with 
respect  to  free  space.  Since  for  a  given  distance  bin 
there  is  a  large  nvimber  of  values  (because  of  a  typically 
veiry  large  nxamber  of  observation  points)  the  values  of  R 
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from  any  given  bin  form  a  random  variable.  The  random 
variables  for  the  individual  distance  bins  can  be  further 
combined  to  form  a  single  random  variable  that  now 
incorporates  the  excess  gain  statistics  for  the  entire 
compartment,  regardless  of  the  distance. 

The  probability  density  function  (pdf)  and  the 
cumulative  distribution  function  (cdf)  for  the  excess  gain 
random  variable  are  compared  to  the  Gaussian  pdf  and  cdf 
with  the  same  mean  and  variance.  The  pdf's  and  the  cdf's 
are  plotted  versus  the  excess  gain  in  dB.  The 

distributions  obtained  via  simulation  are  plotted  on  the 
same  graphs  with  the  Gaussian  distributions,  to  illustrate 
the  "match"  between  the  two.  In  addition,  the  histograms 
for  the  distance  groups  are  also  presented. 

In  order  to  obtain  the  aforementioned  plots,  we  first 
had  to  resolve  an  important  issue.  Invoking  the  volumetric 
plot  command  (VN)  results  in  a  number  of  observation 
points.  The  observation  points  are  uniformly  spaced  and 
thus  would  be  positioned  inside  as  well  as  outside  non- 
penetrable  objects  like  the  cylinders,  the  tanks,  even  the 
steps  of  the  stairways.  Since  the  interiors  of  non- 
penetrable  objects  are  not  considered  as  "legitimate" 
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transceiver  locations,  the  observation  points  within  non- 
penetrable  objects  had  to  be  eliminated. 

A  simple  attempt  to  circumvent  this  problem  by  the 
multiple  selective  use  of  the  same  command  failed  since  the 
code  executed  only  the  last  of  the  coinmands.  Therefore,  we 
had  to  apply  a  rather  time  consuming  procedure:  we 
truncated  the  missile  room  and  the  free  space  area  into  a 
number  of  obstacle-free  observation  regions.  We  ran  the 
program  once  for  each  of  the  observation  regions,  and  then 
"stacked"  (using  Mathcad)  the  data  resulting  from  each  run 
to  form  a  single  "output"  file. 

The  results  of  the  simulations  presented  in  this 
chapter  indicate  that  the  Gaussian  distribution  best 
describes  the  missile  compartment  excess  gain  (in  dB) 
random  variable.  It  has  been  reported  [1]  that  the 
lognormal  distribution  also  applies  to  indoor  propagation 
between  observation  points  that  have  different  distances 
from  the  RF  transmitter. 

Furthermore,  it  has  been  reported  recently  [38]  that 
the  lognormal  distribution  applies  to  the  "same  distance 
points"  that  form  concentric  cycles  around  an  arbitrarily 
placed  source  for  indoor  propagation  within  a  building. 
The  same  approach  has  been  applied  in  this,  thesis  to 


metallic  compartments/  with  similar  results,  that  is  the 
lognormal  distribution  for  the  excess  gain  random  variable. 
Recall  that  the  missile  compartment  is  composed  of  PEC 
materials,  hence  it  can  be  characterized  as  highly 
reflecting/diffracting  environment . 

Due  to  the  absence  of  absorption  and  refraction, 
reflected,  diffracted,  and  "creeping"  waves  dominate  in  the 
propagation  mechanisms.  Each  of  the  multiple  reflections 
that  a  single  "ra/'  experiences  within  the  compartment 
contributes  a  factor  (the  reflection  coefficient  to  the 
magnitude  of  the  ray  at  the  particular  observation  point). 
Therefore  the  observed  ray  magnitude  is  a  product  of  a 
number  of  factors.  The  logarithm  of  the  product  (that  is 
the  path  gain/loss  in  dB)  is  equal  to  the  sum  of  the 
logarithms  of  each  term  in  the  product.  Since  each 
reflection  is  by  and  large  independent  from  the  others 
(involving  different  locations,  objects,  and  incident 
angles)  the  sum  of  logarithms  tends  to  be  a  normally 

distributed  random  variable,  hence  the  lognomal 

distributions . 

2 .  Results 

Let  us  now  examine  the  results  for  several 
transmitting  antenna  locations.  The  observation  points  are 
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always  located  on  a  horizontal  plane,  but  the  observation 
plane  height  {above  the  floor  level)  may  differ  for  each 
individual  case. 

2 . 1  First  Location 

The  transmitting  antenna  is  located  at  (14,5,1)  and 
the  observation  points  lay  in  a  horizontal  plane  at  Z  = 
1.5m,  as  show  in  Figure  29.  The  various  cells  are  shown  in 
Figures  30  through  40. 


Location  of 
Antenna 


Figure  29.  Side  View  of  the  Compartment,  the  Observation 
Area,  Forming  a  Line,  and  the  antenna  at  location  1. 
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Figure  30.  Top  View  of  the  Compartment,  the  First 
Observation  cell,  and  the  antenna  at  location  1. 
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Figure  31.  Top  View  of  the  Compartment,  the  Second 
Observation  cell,  and  the  antenna  at  location  1. 
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Figure  32.  Top  View  of  the  Compartment,  the  Third 
Observation  cell,  and  the  antenna  at  location  1. 
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Figure  33.  Top  View  of  the  Compartment,  the  Fourth 
Observation  cell,  and  the  antenna  at  location  1. 
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Figure  34.  Top  View  of  the  Compartment,  the  Fifth 
Observation  cell,  and  the  antenna  at  location  1. 
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Figure  35.  Top  View  of  the  Compartment,  the  Sixth 
Observation  cell,  and  the  antenna  at  location  1. 
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Figure  36.  Top  View  of  the  Compartment,  the  Seventh 
Observation  cell,  and  the  antenna  at  location  1. 
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Figure  37.  Top  View  of  the  Compartment,  the  Eighth 
Observation  cell,  and  the  antenna  at  location  1. 
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Figure  38.  Top  View  of  the  Compartment,  the  Ninth 
Observation  cell,  and  the  antenna  at  location  1. 
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Figure  39.  Top  View  of  the  Compartment,  the  Tenth 
Observation  cell,  and  the  antenna  at  location  1. 
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Figure  40.  Top  View  of  the  Compartment,  the  Eleventh 
Observation  cell,  and  the  antenna  at  location  1. 

The  excess  gain  (in  dB)  pdf  and  cdf,  obtained  by 
averaging  over  all  distance  groups,  are  shown  in  Figures  41 
and  42,  respectively.  Plotted  on  the  same  graphs  are  the 
lognormal  (Gaussian)  pdf  and  cdf,  to  visualize  the 
agreement  between  the  lognormal  distribution  model  and  the 
data  obtained  by  simulation. 
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Antenna  at  location  1 
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Figure  41.  Pdf  of  the  Excess  Gain  and  the  Lognormal 
(Gaussian)  Pdf,  Transmitting  Antenna  at  Location  1 


Antenna  at  location 


Figure  42.  Cdf  of  the  Excess  Gain  and  Lognormal  (Gaussian) 
Cdf/  Antenna  at  Location  1 

As  shown,  there  is  a  good  match  between  the  excess 
gain  pdf  (obtained  by  simulation)  and  the  lognormal  pdf. 
The  oscillations  ("dispersion")  of  the  simulation  results 
about  lognormal  distribution  are  inherent  to  the 
statistical  nature  of  the  simulation  results.  The 
"dispersion"  can  be  reduced  by  increasing  the  number  of 
observation  points  but  at  the  expense  of  increasing  the 
computational  times.  The  agreement  is  even  better  for  the 
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cdf's,  V7hich  could  have  been  expected  as  the  integration 
"averages  out"  the  statistical  dispersion  of  the  pdf.  Note 
that  the  pdf  oscillations  are,  on  the  average,  equally 
likely  to  be  either  above  or  below  the  lognormal  pdf,  hence 
they  tend  to  cancel  out  when  the  pdf  is  integrated  to 
obtain  the  cdf.  From  the  cdf  we  note  that  the  compartment 
(for  this  particular  antenna  location)  has  a  0.5 
probability  of  behaving  exactly  like  free  space;  that  is, 
it  experiences  neither  gain  nor  loss  relative  to  free  space 
propagation. 

2 . 2  Second  Location 

In  this  case  the  transmitting  antenna  is  located  at 
(2,2,2)  and  the  observation  points  are  in  a  horizontal 
plane  at  2  =  2m.  The  results  are  illustrated  in  the  Figures 
4  3  and  4  4 . 
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this  "thresholding'^  causes  a  distinct  peak  to  appear  in  the 
pdf,  distorting  the  pdf's  shape  relative  to  the  expected 
lognormal  distribution.  To  illustrate  the  effects  of  the 
dynamic  range  limiting  on  the  pdf  obtained  via  simulation, 
the  pdf  for  the  40-d3  dynamic  range  is  shown  in  Figure  45. 
Compared  with  the  30-dB  dynamic  range  pdf,  the  40-dB 
dynamic  range  pdf  shows  a  more  promdnent  peak,  shifted 
further  to  the  left  (larger  values  of  excess  path  loss). 

In  order  to  mitigate  the  effects  of  the  dynamic  range 
limitation  inherent  to  the  NEC-BSC  (as  a  mater  of  fact  all 
simulation  codes  have  a  finite  dynamic  range,  due  to  the 
"numerical  noise")  a  "fix"  has  been  implemented  in  Mathcad. 
The  essence  of  the  "fix"  was  to  replace  the  value  of  the 
"offending"  peak  with  the  value  that  corresponds  to  the 
lognormal  distribution  and  then  scale  the  resulting 
distribution  such  that  its  integral  is  equal  to  1  (as  it 
has  to  be  for  any  pdf).  Shown  in  Figure  45  are  the  pdf's 
for  the  30  and  40  dB  dynamic  ranges  and  the  "corrected"  (in 
the  m^anner  described  above)  pdf  for  the  30-dB  dynamic 
range.  It  is  evident  that  the  outlined  procedure  results 
in  a  pdf  that  is  again  close  to  the  lognormal  pdf. 
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2 . 3  Third  Location 


In  this  case  the  transmitting  antenna  is  located  at 
(2, 2, 2.2)  and  the  observation  points  are  in  a  horizontal 
plane  at  Z  =  2.2m  (at  the  same  height  as  the  transmitting 
antenna)  .  Note  that  the  difference  between  this  and  the 
preceding  case  is  that  the  antenna  and  the  observation 
points  have  been  translated  upward  by  20  cm.  The  reason  to 
do  this  is  to  examine  the  consequences  of  a  translation  by  a 
relatively  small  distance,  coirmensurate  with  the  operating 
wavelength  (approximately  1.5  wavelengths  in  this  case).  In 
addition,  the  observation  points  do  not  fall  on  the  wide 
horizontal  steps  of  the  stairway.  The  results  are  shown  in 
Figures  46  and  47. 
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Third  location  of  antenna 
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Figure  46.  Pdf  of  the  Excess  Gain  and  the  Lognormal 
(Gaussian)  Pdf,  Transmitting  Antenna  at  Location  3. 


Room  gain  (db) 


Figure  47.  Cdf  of  the  Excess  Gain  and  the  Lognormal 
(Gaussian)  Cdf,  Antenna  at  Location  3. 


2.4  Fotarth  Location 


In  this  case  the  transmitting  antenna  is  located  at 
(4,3,1)  and  the  observation  points  are  in  a  horizontal  plane 
at  Z  =  1  m  (at  the  same  height  as  the  transmitting  antenna) . 
This  is  the  last  case  considered,  for  which  the  transmitting 
antenna  is  placed  at  an  "arbitrary^'  location  (not  in  a  plane 
of  syircmetry,  or  the  like) .  However,  the  observation  points 
are  still  in  a  horizontal  plane,  at  the  same  height  as  the 
transmitting  antenna.  Following  the  same  procedure  as  for 
the  preceding  three  locations,  we  obtain  the  results  shown 
in  Figures  50  and  51. 
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Figure  50.  Pdf  of  the  Excess  Gain  and  the  Lognormal 
(Gaussian)  Pdf,  Transmitting  Antenna  at  Location  4. 
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Figure  51.  Cdf  of  the  Excess  Gain  and  the  Lognorxaal 
{Gaussian)  Cdf,  Antenna  at  Location  4. 


The  data  for  the  four  antenna  locations  used  in  the 
simulations  are  presented  in  Table  5.1. 


Antenna 

Antenna 

[  Observation 

1 

location 

Coordinates  ; 

[  points 

Excess 


Excess  gain 


Second 


Third 


2,2,2 


1.5 


Z=2 


gain  St.  Deviation 


Mean  (dB) 


-0.10 


1.43 


2.91 


0.5 


(dB) 


5.31 


4.85 


3.95 


3.9 


’able  5.1.  Data  for  the  Four  Propagation  Simulations. 
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VT.  UNRESOLVED  ISSUES,  RECOMMENDATIONS,  AND 

CONCLUSIONS 

A.  UNRESOLVED  PROBLEMS 

1 .  Floor- to-Floor  Propagation 

The  case  in  which  the  antenna  was  located  inside  the 
first  floor  of  the  missile  room,  and  the  observation  points 
in  the  second  (and  vice  versa),  conld  not  be  thoroughly 
examined  due  to  the  problems  caused  by  the  "ray  leakage," 
as  reported  in  Chapter  IV.  Therefore  the  compartment-to- 
compartment  propagation  statistics  were  not  addressed. 
Note  that  the  limited  dynamic  range  is  also  an  issue  with 
regard  to  compartment-to-compartment  coupling  that  "by 
definition"  involves  large  attenuation. 

2.  Problematic  PN  Command 

The  Near  Zone  Pattern  (PN)  command  enabled  us  to 
define  the  near  zone  pattern  coordinate  system,  the 
volumetric  pattern  cut,  and  the  spatial  range  of  desired 
observation  points.  However,  as  reported  in  Chapter  V,  the 
multiple  and  selective  use  of  the  PN  command  to  exclude  the 
undesired  observation  points  (inside  non-penetrable 
objects)  was  unsuccessful,  because  the  code  executed  only 
the  last  of  the  multiple  PN  commands.  The  truncation  of 
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the  room  into  several  regions,  exclusive,  of  non-penetrable 
objects,  did  work  but  at  the  ej^ense  of  excluding  some 
points  located  near  the  non-penetrable  obstacles  and  it 
increased  time  requirements  for  processing  of  multiple  data 
files. 

3 .  Running  Time 

All  simulation  programs  executed  in  Chapter  V  were 
time  consuming  and  had  run  times,  ranging  from  seven  to 
twelve  hours.  As  expected,  time  reguirem.ents  depended  on 
the  complexity  of  the  coirpartment ,  the  number  of  the 
observation  points  and  the  nuit^Der  of  antennas  used. 

B.  RECOMMENDATIONS 

All  recommendations  for  further  study  are  made  in  the 
light  of  existing  information,  based  on  NEC-BSC  Workbench 
version  4.1.37  and  executable  NEC-BSC  FORTRAN  program, 
versions  4.2-06. 

Some  commands,  namely  BF,  FT,  LY,  NC,  NG,  NM,  NP,  NR, 
NS,  NX,  PT,  RI,  SI,  TA,  TD,  TW,  TY,  VM,  were  not  supported 
by  the  Coimnand  Wizard  window  of  the  code's  Graphical  User 
Interface  (GUI) .  Moreover,  the  function  of  the  Command  ZC, 
provided  by  Workbench,  remains  mysterious,  since  there  is  no 
description  for  it,  either  "on-line"  or  in  the  User's 
manual . 
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In  addition,  the  Zoom  In  command  in  the  main  toolbar 
was  inconvenient,  because  the  user  was  not  able  to  specify 
the  desired  scale,  which  by  default  was  inadequately  large. 
The  accuracy  of  the  Zoom  was  also  questionable,  since  it 
altered  the  location  of  certain  objects,  presented  the 
cylinder  caps  as  polygons,  and  showed  some  rays  diffracted 
from  cylinders  as  if  the  rays  penetrated  the  cylinders  and 
then  diffracted  from  the  inside  of  the  cylinders. 

It  should  also  be  stated  that  the  very  useful  and 
convenient  "drag  and  drop"  option  couldn't  be  used  in  all 
cases;  the  user  could  not  always  distinguish  and  select  an 
individual  ray  among  a  plethora  of  others.  Finally,  the 
coordinate  system,  of  the  code,  which  by  default  placed  the 
origin  in  the  left  upper  corner  of  the  display,  used 
unknown  units  and  was  difficult  to  work  with. 

A  recent,  extensive  study  [38]  includes  a  number  of 
recommendations  that  would  significantly  improve  the 
applicability  of  the  NEC-BSC  for  indoor  propagation 
simulation.  In  addition  to  these,  parallel  execution 
option  [37]  would  also  decrease  the  run-times,  especially 
in  complex  cases  that  approach  the  limitations  of  the  code, 
as  reported  in  Chapter  III. 
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Moreover,  a  considerable  improvement  in  the 
calculation  speed  and  memory  requirements  could  be  made  if 
the  near-field  computation  by  NEC-BSC  did  not  include  both 
the  Electric  and  the  Magnetic  field  components 
simultaneously.  In  situations  where  the  computation  of  the 
signal  strength  was  the  goal,  the  calculation  of  the  H- 
fieid  could  not  be  skipped.  The  user  should  be  able  to 
choose  the  particular  quantity  that  the  code  should 
calculate  in  order  to  reduce  the  conputational  time. 

Another  factor  that  would  improve  code  performance  is 
the  support  of  ail  commands  in  the  Workbench.  In  addition, 
the  Users  Manual  should  be  completed  and  should  include  all 
the  figures  omitted  in  the  previous  versions,  along  with 
more  examples  that  illustrate  the  workflow. 

Finally,  the  creation  of  the  problem  geometry  will  be 
greatly  simplified  if  future  versions  of  NEC-BSC  would 
support  data  import  in  a  standard  CAD  file  format  {.dxf, 
.dwg,  etc.). 

C.  CC»?CLUSIC»JS 

The  primary  objective  of  this  research  was  to  provide 
simple  statistical  models  for  indoor  propagation  in  complex 
environments,  such  as  shipboard  compartments.  The  emphasis 
in  this  research  was  on  the  basic  principles  and  a  simple 
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statistical  description  for  a  representative  shipboard 
environment  rather  than  on  encyclopedic  compilation  for 
indoor  propagation  data  for  vastly  diverse  complex 
environments .  Previous  work  using  NEC-BSC  for  indoor 
propagation  sim.ulation,  involving  environments  with  no 
material  uniforrrdty  [38]  showed  satisfactory  agreement  with 
both  the  results  reported  in  140] ,  and  with  the 
measurements  conducted  in  the  Bullard  Hall  at  the  Naval 
Postgraduate  School  in  Monterey,  California  [38] . 

In  our  case,  as  well  as  for  the  Bullard  Hall  case 
[38],  the  '^room  excess  gain"  random  variable  v/as  defined  to 
describe  hov;  the  propagation  in  the  indoor  environment 
differs  from  the  deterministic  (formula-predictable)  free 
space  propagation.  The  simulation  results  in  both  cases 
show  that  the  excess  gain  (in  dB)  random  variable  could  be 
approximated  as  lognormal  (Gaussian)  distributed  [39] . 
Note  that  this  random  variable  is  representative  for  all 
distance  combinations  and  locations  of  the  source  and 
observation  points  within  a  single  compartment.  Although 
the  shape  of  the  probability  density  function  for  the 
excess  gain  (in  dB)  random  variable  does  not  depend  on  the 
particular  choice  of  the  source  location  or  the  compartment 
details,  the  parameters  defining  the  pdf,  such  as  the  mean 
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value  and  the  variance,  do  depend  on  the  position  of  the 
source  and  the  observation  points  relative  to  the  obstacles 
present  in  the  compartment - 

It  appears  that  the  key  parameter  in  defining  the  mean 
value  of  the  excess  gain  is  the  relative  number  of 
observation  points  that  are  shadowed  by  the  obstacles  and 
thus  experience  non-line-of-sight  propagation  attenuation. 
If  there  is  a  large  number  of  observation  points  in  deep 
shadows,  which  depends  primarily  on  the  location  of  the 
source  relative  to  the  obstacles  in  the  compartment,  the 
mean  value  of  the  excess  gain  is  negative,  indicating  a 
loss  {on  the  average)  relative  to  the  free  space 
propagation.  If  on  the  other  case  most  of  the  observation 
points  are  either  direct  line-of-sight  or  are  obstructed  by 
low-loss  penetrable  obstacles,  the  mean  value  of  the  excess 
gain  (in  dB)  becomes  positive,  indicating  that  the  signal 
strength  in  that  compartment  exceeds,  on  average,  signal 
strength  in  free  space - 

The  later  can  be  physically  understood  by  observing 
that  the  signal  energy  in  free  space  propagates  outward, 
with  the  corresponding  decrease  in  the  power  density,  while 
in  an  indoor  environment  -  some  of  the  energy  is  reflected 
back  by  the  walls,  window  panes,  etc.,  thus  increasing  the 
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average  power  density  within  the  compartment.  The 
particular  compartment  analyzed  in  this  research  and  for 
the  several  transmitting  antenna  locations  is  by  and  large 
"free  space"  neutral,  with  the  excess  gain  mean  value 
between  -0.1  and  +2.9  dB,  depending  on  the  source  location. 
The  restrictions  of  tim.e  and  computer  resources  limited 
this  research  to  only  one  complex  compartment. 
Furthermore,  the  computer  resources  available  would  not 
allow  the  placement  of  observation  points  at  separations 
below  "two  per  operating  wavelength"  to  allow  very  fine 
sampling  of  the  fields  inside  the  compartments. 

These  facts,  as  well  as  the  code  Imitation  with 
respect  to  the  dynamic  range,  should  be  taken  into  account, 
and  offer  plenty  of  opportunities  for  further  research. 
Nonetheless,  the  results  definitely  confirmed  the 
applicability  of  the  excess  gain  random  variable  approach  in 
general  and  the  lognormial  distribution  model  in  particular. 
NEC-BSC  is  generally  accepted  as  a  very  useful  tool  for 
modeling  a  variety  of  radiation  and  scattering  phenomena. 
Over  the  years  it  was  reported  to  provide  reliable  results 
especially  when  used  for  outdoor  environments.  Recently, 
other  studies  (  [38]  and  the  present  one)  showed  quite 
promising  results  in  using  NEC-BSC  for  modeling  indoor 
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propagation  as  well.  However,  definite  conclusions  require 
extensive  further  investigation  and  modeling  of  a  multitude 
of  diverse  indoor  environments. 
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APPENDIX  A.  MATIAB  FILES 

This  is  the  Matlab  file  used  to  read  the  data,  resulted 
from  "PATCH"  software  and  create  the  plots  necessary  for 
the  comparison  between  that  and  NEC-BSC  (Case:  cylinder  and 
dipole) 

%  plot  near  field  data 
clear 

rad=pi/180; 
iff=l; 
if  iff==0 
load  amptffl.m 
load  amppffl.m 

amptff=amptf fl ;  amppf f=amppf f 1; 

?= [0:180]  ; 
figure (1) , clf 

plot(P,amptff  (1:181)  ,P,amppff  (1:181) ,  '  — ')  ,grid 
xlabel ( ’ Phi ,  deg ' ) 

ylabel (' Total  electric  field  (dB) ' ) 

title ('Far  field,  a=0.175m,  L=lm,  theta=90,  dip  @  phi=0') 

axis( [0,180,-30,0] ) 

end 


%  currents  (get  from  outnear  files) 

Js=0.17859E+01-^-j*0.20941E4-01;  iplt=0;  %  2.54  GHz  outnear 

2 

%  Js=0.0041847-fj*0.40483E^00;  iplt=0;  %  1  GHz  outnear  3 
%  Js=0.18605EH-01-j *0.83378;  iplt=l;  '  %  1  GHz  dipole 
outnear  6  (20  segs) 

%  Js=0.96774E+01-j*0.48566E+01;  iplt=l;  %  1  GHz  skinny 

dipole  20  segs  outnear  7 

%  Js=2.4268-j*0. 62429;  iplt=0;  %  outnear  9  (longer  dip 
for  1  GHz) 

I=Js*.005;  %  I=Js* (width  of  dipole) 

V=l; 

Z=V/I; 

%  D=20*logl0 (abs (I) * -19) ;  %  for  outnear  3  L  is  not  wavl/2 

but  0.19 

D=20*logl0 (abs (I) ) ;%  factor  of  nonunit  amplitudecurrent 
load  nearpts2.m. 
nearpts=nearpts2 ; 

[nrow,ncol] =size (nearpts) ; 
npts=nrow/2 ; 
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Ez ( 1 : npts ) =nearpts (1:2: nrow, 7 ) + j  *nearpts (1:2: nrow, 8 ) ; 

Ey  ( 1 :  npts )  =nearpts  (1:2: nrow,  5 )  +  j  *nea3rpts  (1:2  :nrow,  6 )  ; 

Ex ( 1 ; npts ) =nearpts (1:2: nrow, 3 ) + j  *nearpts (1:2: nrow, 4 ) ; 

if  iplt==0, 

load  nearfldptslSO.m;  nearfldpts=n€arfldptsl80 ; 
end 

if  iplt==l,  load  nearfldptsz.m;  nearfldpts=nearfldptsz; 
end 

x=nearf Idpts ( : , 1 ) ;  y=nearf Idpts ( : , 2 ) ;  z=near f Idpts ( : , 3 ) ;  % 
set  for  180  points 
R=sqrt(x.''2+y.^2+z.''2)  ; 

T=acos (z. /R) ; 

P=atan2  (y,x)  ; 

U=cos (T) . *cos (P) ; 

V=cos(T) .*sin{P) ; 

W=-sin(T} ; 

Y=cos (P) ; 

Z=sin(P) ; 

Et=U.  *■  transpose  (Ex)  +V.  *  transpose  (Ey)  +W.  *  transpose  (Ez)  ; 

Ep=-Z .  *transpose  (Ex)  -fY .  *  transpose  (Ey)  ; 

Emax=sqrt (Et . *conj (Et) ) ;  %  +Sp. *conj (Ep) ) ; 

%  Emax=sqrt  (Ex.  *conj  (Ex)  -rEy .  *conj  (Ey)  4-Sz  .  *conj  (Sz) )  ; 
Edb=20*logl0 (Emax) -D; 
load  gtddat.m 

E=gtddat(: ,6) ;  E (361) =E (1) ; 

if  iplt==0  %  plot  cylinder  and  dipole 

E2=[Edb  fliplr (Edb) ] ;  %  mirror  image  of  pattern 

?2= [P/rad  360-fliplr (P/rad) ] ; 
figure (2) , clf 

plot {P2,E2, 'h' , [0:360] ,E, 'r--') ,grid 
xlabel ( ' Phi ,  deg ' ) 
ylabel {' Total  electric  field  (dB)') 

title { 'Electric  field  by  patch  (blue)  and  NEC -BSC  (red)') 
end 


if  iplt==l  %  just  plot  dipole  data 

for  i=l: length (T) 

TT(i)=T(i) ; 

if  P(i)==pi,  TT(i)=2*pi-T(i) ; 

end 

end 
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plot  (TT/rad, Edb, '  - ' ) /ilold  on 
xlabel { ' Theta ,  deg ' ) 
ylabeK'  Ie]  ,  dB') 

title  ( 'half -tvave  dipole  field  at  R=5in  (f=l  GHz)') 
axis{[0,360,-20,30]) 
figure (3) , clf 

polardb { (TT/rad) . ' , Edb, ' - ' ) 

title {'short  dipole  field  at  R=5m  (f=l  GHz)') 

%  analytical  dipole  field  (current  is  I  from  above) 
Et=377/2/pi/5*abs (cos {pi/2 . *cos (TT) ) . /sin(TT+le-3) ) 
Etdb=20*logl0(Et) ; 
figure (2) 

plot (TT/rad. Etdb, ' — ' ) 
axis( [0.360,-20,30] ) 
hold  off 
end 

End  of  all  Matlab  files 
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APPENDIX  B.  NEC-BSC  INPUT  FILES 


CE:  This  is  the  input  file  used  in  the  con5)arison  between 
CE:  NEC-BSC  ar>ri  PATCH  software  in  the  case  of  the  cylinder 
CE:  and  the  dipole 

ZD: 

0 

DN: 

1 

FR:  2.54 
2.54000 


US :  v/avelengths 
0 

SG: 


1.00000, 

0.000000, 

0.000000 

0.000000, 

0.000000, 

90.0000, 

0.000000 

-2,  0.500000,  0.000000 
1-00000,  0.000000 

CG: 

0.000000, 

0.000000, 

0.000000 

0.000000, 

0.175000, 

0.000000, 

0.175000 

90.0000, 

0.000000 

-0.500000 

,  90.0000, 

0.500000, 

90.0000 

PN: 

0.000000, 

0.000000, 

0.000000 

0.000000, 

F 

0.000000, 

90.0000, 

0.000000 

5.00000, 

90.0000,  0 

.000000 

0.000000, 

0.000000, 

1.00000 

360 

PL: 

T 

LY 

T 

XQ:  Execute  the  program 

EN:  End 


107 


CEs  ISiis  is  the  input  file  used  in  the  simple  metallic  room 

CE:  Six  box  room  with  a  cylindrical  obstacle  inside  it 

ZD:  type  of  material:  PEC 
0 

?R:  we  will  use  a  frequency  of  2.4  GH 
2.40000 

US: 

0 


SG:  Source  geometry  (3,1,11). 

3.00000,  1.00000,  11-0000 
0.000000,  0.000000,  90.0000,  0.000000 
-2,  0.500000,  0.000000 
1.00000,  0.000000 


PG:  First  plate 
4,  0  . 

o.oooboo,  0.000000,  0.000000 
0.000000,  12.0000,  0.000000 
6.00000,  12.0000,  0.000000 
6.00000,  0.000000,  0.000000 

PG:  second  plate 
4,  0 

0.000000,  0.000000,  0.000000 
6.00000,  0.000000,  0.000000 
6.00000,  0.000000,  12.0000 
0.000000,  0.000000,  12.0000 

PG :  Third  plate 
4,  0 

6.00000,  0.000000,  12.0000 
6.00000,  0.000000,  0.000000 
6.00000,  12.0000,  0.000000 
6.00000,  12.0000,  12.0000 

PG:  Forth  plate  ==front 
4 ,  0 

6.00000,  12.0000,  12.0000 
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6.00000,  12.0000,  0.000000 
0.000000,  12.0000,  0.000000 
0.000000,  12.0000,  12.0000 

PG:  Fifth  plate 
4,  0 

6.00000,  0.000000,  12.0000 
6.00000,  12.0000,  12.0000 
0.000000,  12.0000,  12.0000 
0.000000,  0.000000,  12.0000 


PG:  Sixth  plate 
4,  0 

0.000000,  0.000000,  12.0000 
0.000000,  12.0000,  12.0000 
0.000000,  12.0000,  0.000000 
0.000000,  0.000000,  0.000000 

CG:  Put  a  cylinder  as  an  obstacle  inside  the  room 
3.00000,  6.00000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
1.00000,  1.00000 

0.000000,  90.0000,  6.00000,  90.0000 


PN:  Near  zone  cut 
3.00000,  6.00000,  3.00000 
0.000000,  0.000000,  90.0000,  0.000000 
F 

1.50000,  90.0000,  0.000000 
0.000000,  0.000000,  1.00000 
180 

PL: 

T 

LY 

T 

XQ:  Execute  the  program 

EN:  End 

CE:  This  is  the  input  file  used  in  the  con^lex  metallic 
room 

CM:  Six  box  room  (6,12,12)  with  two  cylindrical  obstacles 

CE:  and  two  tanks  inside 
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ZD:  type  of  material:  PEC 
0 

FR:  we  will  use  a  frequency  of  2.45  GH 
2.45000 

US: 

0 


SG :  Source  geometry  (3,3,1) . 

3.00000,  3.00000,  1.00000 
0.000000,  0.000000,  90.0000,  0.000000 
-2,  0.500000,  0.000000 
1.00000,  0.000000 

VN:  VOLUMETRIC  PLOT  at  Z=1 . 5  m. 

0.000000,  0.000000,  0.000000 

0.000000,  0.000000,  90.0000,  0.000000 

T 

0.000000,  0.000000,  1.50000 
0.100000,  0.100000,  0.000000 
1,  60,  120 
PG:  First  plate 
4,  0 

0.000000,  0.000000,  0.000000 
0.000000,  12.0000,  0.000000 
6.00000,  12.0000,  0.000000 
6.00000,  0.000000,  0.000000 


PG:  second  plate  . 

4,  0 

0.000000,  0.000000,  0.000000 
6.00000,  0.000000,  0.000000 
6.00000,  0.000000,  12.0000 
0.000000,  0.000000,  12.0000 


PG:  Third  plate 
4,  0 

6.00000,  0.000000,  12.0000 
6.00000,  0.000000,  0.000000 
6.00000,  12.0000,  0.000000 
6.00000,  12.0000,  12.0000 


PG:  Forth  plate  ==front 
4,  0 
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6.00000,  12.0000,  12.0000 

6.00000,  12.0000,  0.000000 

0.000000,  12.0000,  0.000000 

0.000000,  12.0000,  12.0000 


PG:  Fifth  plate 
4,  0 

6.00000,  0.000000,  12.0000 
6.00000,  12.0000,  12.0000 
0.000000,  12.0000,  12.0000 
0.000000,  0.000000,  12.0000 

PG:  Sixth  plate 
4,  0 

0.000000,  0.000000,  12.0000 
0.000000,  12.0000,  12.0000 
0.000000,  12.0000,  0.000000 
0.000000,  0.000000,  0.000000 

CG:  Put  a  cylinder  as  an  obstacle  inside  the  room 
1.50000,  6.00000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
O'.SOOOOO,  0.500000 

0.000000,  90.0000,  12.0000,  90.0000 


CG:  Put  a  second  cylinder  as  an  obstacle  inside  the  room 
4.50000,  6.00000,  0.000000 
0.000000,  0.000000,  90.0000,  0-000000 
0.500000,  0,500000 

0.000000,  90.0000,  12.0000,  90.0000 


PG:  first  tank  front 
4,  0 

2.00000,  1.00000,  0.000000 
2.00000,  1.00000,  6.00000 
4.00000,  1.00000,  6.00000 
4.00000,  1.00000,  0.000000 


PG:  first  tank  side 
4,  0 

4.00000,  0.000000,  6.00000 
4-00000,  0.000000,  0.000000 
4.00000,  1.00000,  0.000000 
4.00000,  1.00000,  6.00000 

PG:  first  tank  side 
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4,  0 

2.00000,  0.000000,  6.00000 
2.00000,  1.00000,  6.00000 
2.00000,  1.00000,  0.000000 
2.00000,  0.000000,  0.000000 


PG:  first  tank  up 
4,  0 

2.00000,  0.000000,  6.00000 
4.00000,  0.000000,  6.00000 
4.00000,  1.00000,  6.00000 
2.00000,  1.00000,  6.00000 


PG:  second  tank  front 
4,  0 

2.00000,  11.0000,  0.000000 
4.00000,  11.0000,  0.000000 
4.00000,  11.0000,  6.00000 
2.00000,  11.0000,  6.00000 

PG:  second  tank  side 
4,  0 

4.00000,  11.0000,  6.00000 
4.00000,  11.0000,  0.000000 
4.00000,  12.0000,  0.000000 
4.00000,  12.0000,  6.00000 

PG:  second  tank  side 
4,  0 

2.00000,  11.0000,  6.00000 
2.00000,  12.0000,  6.00000 
2.00000,  12.0000,  0.000000 
2.00000,  11.0000,  0.000000 


PG:  second  tank  uo 


4,  0 
2.00000, 
4.00000, 
4.00000, 
2.00000, 

PL: 

T 

XQ: 


11.0000, 

11.0000, 

12.0000, 

12.0000, 


6.00000 

6.00000 

6.00000 

6.00000 


112 


EN:  End 


CE:  This  is  the  input  file  used  for  openings 

CE:  right  hand  ok,  commands  ok  problem  with  penetration 

FR:  We  will  use  a  frequency  of  2.4  GHz  (ISM  band) 
2.40000 


US: 

0 


PG :  platel  up 
4 ,  0 

6.00000,  -6.00000,  12.0000 
6.00000,  -6.00000,  8.00000 
6.00000,  6.00000,  8.0000 
6.00000,  6.00000,  12.00000 


PG:  plate  2 
4,  0 

6.00000,  6.00000, 
6.00000,  3.00000, 
6.00000,  3.00000, 
6.00000,  6.00000, 


12.0000 

12.0000 

0.000000 

0.000000 


PG:  plate  3 
4,  0 

6.00000,  6.00000,  0.000000 
6.00000,  6.00000,  4.00000 
6.00000,  -6.00000,  4.00000 
6.00000,  -6.00000,  0.000000 


PG:  plate  4 
4,  0 

6.00000,  -6.00000,  0.000000 
6.00000,  -3.00000,  0.000000 
6.00000,  -3.00000,  12.0000 
6.00000,  -6.00000,  12.0000 

SG:  Source  in  height  6 
0.000000,  0.000000,  6.000000 
0.000000,  0.000000,  90.0000,  0.000000 
-1,  1.00000,  0.000000 
1.00000,  0.000000 
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PN:  Near  zone  cut 


9.00000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
T 

3.00000,  0.000000,  0.000000 
0.000000,  0.000000,  2.00000 
7 

PL: 

T 

LY 

T 

XQ:  Execute  the  program 

EN:  End 

CE:  ®his  is  the  input  file  used  for  pseudo-openings 

CE:  right  hand  ok 

PR:  We  will  use  a  frequency  of  2.4  GHz  (ISM  band) 
2.40000 

US: 

0 


PG:  wall 
4,  0 

6.00000,  -6.000000,  12.0000 
6.00000,  -6.000000,  0.000000 
6.00000,  6.0000,  0.000000 
6.00000,  6.0000,  12.0000 


ZD:  opening- thin  material  slab 
1 
1 

0.0500000,  1.01000,  0.000000,  1.00000,  0.000000 

PG:  door 
4,  0 

6.00000,  -3.000000,  8.0000 
6.00000,  -3.000000,  4.000000 
6.00000,  3.00000,  4.000000 
6.00000,  3.0000,  8.0000 


SG:  Source  in  height  6 

0.000000,  0-000000,  6.00000 
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0.000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0.000000 
1.00000,  0.000000 

PN:  Near  zone  cut 
9.00000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
T 

3.00000,  0.000000,  0.000000 
0.000000,  0.000000,  2.00000 
7 

PL: 

T 

LY 

T 

XQ:  Execute  the  program 

EN:  End 

CE:  This  is  the  input  file  used  for  one  Antenna  located  in 
CE:  the  middle  of  the  room 

CE:  Rear  tank  plates  overlap  to  avoid  diffraction 
CE:  inside  them 

CM:  Dimensions  of  the  massile  room:  Dim(20, 10, 10)  Origin  at 
CE:  the  most  FORWARD  and  PORT  point  of  the  default  view 

CM:  All  materials  would  be  Perfect  Electric  Conductors, 

CE:  unless  otherwise  defined. 

ZD:  P.E.C. 

0 

FR:  We  will  use  a  frequency  of  2.4  GHz  (ISM  band) 

2.40000 

US:  (Necessary)  in  Wavelengths 
0 


SG:  Put  the  source  at  (10,5,2)  position 
10.0000,  5.00000,  2.00000 
0.000000,  0.000000,  90.0000,  0.000000 
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-1,  0.500000.  0.000000 

1.00000.  0.000000 

VN:  VOLUMETRIC  PLOT  at  Z=2  m. 

0.000000.  0.000000,  0.000000 
0.000000,  0.000000.  90.0000,  0.000000 
T 

0.000000,  0.000000,  2.00000 
0.0500000,  0.0500000,  0.000000 
1,  400,  200 

CE:  Construct  the  room,  according  to  the  RIGHT-HMD  RULE 

CM:  Since  we  can  not  simulate  the  opening  with  an  adequate 
CM:  material,  we  define  five  different  plates  for  the 
CE:  Forward  Level. 

CM:  PG:  room: (20, 10, 10)  FORWARD  level 
CM:  4,  0 

CM:  0.000000,  10.0000,  0.000000 
CM:  0.000000,  10.0000,  10.0000 
CM:  0.000000,  0.000000,  10.0000 
CE:  0.000000,  0.000000,  0.000000 

CE:  Construct  all  the  plates,  but  the  opening  itself 


PM: 

12 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0 


0.000000 

10.0000, 

10.0000, 

0.000000 

5.50000, 

10.0000, 

10.0000, 

0.000000 

5.50000, 
0.000000 

4.50000, 

4.50000, 


0.000000 

0.000000 

0.300000 

0.300000 

1.80000 

1.80000 

10.0000 

10.0000 

0.300000 

1.80000 

0.300000 

1.80000 


1 

2 

3 

4 
0 
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9 

3 
6 

5 
0 

4 

11 

12 

10 
0 

10 

6 

7 

8 
0 
0 

CM:  In  case  we  want  the  wall  without  the  opening  (solid) , 
CE:  we  uncomment  the  following  command: 


CM:  PG:  This  is  the  opening  in  the  forward  bulkhead 
CM:  4,  0 

CM:  0.000000,  4.50000,  0.300000 
CM:  0.000000,  4.50000,  1.80000 
CM:  0.000000,  5.50000,  1.80000 
CE:  0.000000,  5.50000,  0.300000 

PG:  room: (20,10,10)  GROUND  level 
4,  0 

0.000000,  0.000000,  0.000000 
20.0000,  0.000000,  0.000000 
20.0000,  10.0000,  0.000000 
0.000000,  10.0000,  0.000000 


PG:  room:  (20,10,10)  UPPER  level 
4,  0 

0.000000,  10.0000,  10.0000 
20.0000,  10.0000,  10.0000 
20.0000,  0.000000,  10.0000 
0.000000,  0.000000,  10.0000 


PG:  room:  (20, 10, 10)  AFT  level 
4,  0 

20.0000,  0.000000,  0.000000 
20.0000,  0.000000,  10.0000 
20.0000,  10.0000,  10.0000 
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20-0000,  10.0000,  0.000000 


?G:  PORT  BULKHEAD 
4,  0 

0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  10.0000 
20.0000,  0.000000,  10.0000 
20.0000,  0.000000,  0.000000 

PG:  STARBO.ARD  BULKHEAD 
4,  0 

0.000000,  10.0000,  0.000000 
20.0000,  10.0000,  0.000000 
20.0000,  10.0000,  10.0000 
0.000000,  10.0000,  10.0000 

CE:  First  OVERHEAD:  This  is  an  intermediate  floor  at  z 

CE : level 


PG:  level: (0, 0, 5)  first  walking  floor 
4,  0 

0.000000,  0.000000,  5.00000 
0.000000,  10.0000,  5.00000 
20.0000,  10.0000,  5.00000 
20.0000,  0.000000,  5.00000 

CE:  Construct  all  the  plates,  but  the  opening  itself 


CM: 

PG:  RIGHT 

CM: 

4,  0 

CM: 

0.000000,  5.50000, 

5.00000 

CM: 

0.000000,  10.0000, 

5.00000 

CM: 

20.0000,  10.0000, 

5.00000 

CM: 

CM: 

20.0000,  5.50000, 

5.00000 

CM: 

PG:  MIDDLE 

CM: 

4,  0 

CM: 

15.7500,  5.51000, 

5.00000 

CM: 

20.0000,  5.51000, 

5.00000 

CM: 

20.0000,  4.49000, 

5.00000 

CM: 

CM: 

15.7500,  4.49000, 

5.00000 

CM: 

PG:  MIDDLE 

CM: 

4,  0 

CM: 

0.000000,  5.51000, 

5.00000 

CM: 

14.7500,  5.51000, 

5.00000 
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CM:  14.7500,  4.49000,  5.00000 
CM:  0.000000,  4.49000,  5.00000 
CM: 

CM:  ?G:  LEFT 
CM:  4,  0 

CM:  0.000000,  0.000000,  5.00000 
CM:  0.000000,  4.50000,  5.00000 
CM:  20.0000,  4.50000,  5.00000 
CE:  20.0000,  0.000000,  5.00000 

CM:  In  case  we  want  the  overhead  without  the  opening 
CM: (solid) , 

CE:  we  uncomment  the  following  command: 

CM:  PG:  This  is  an  opening  in  the  overhead 
CM:  4,  0 

CM:  15.7500,  4.50000,  5.00000 
CM:  14.7500,  4.50000,  5.00000 
CM:  14.7500,  5.50000,  5.00000 
CE:  15.7500,  5.50000,  5.00000 


CE:  These  are  the  6  missile  txibes 

CG:  First  (FORWARD  PORT)  missile  tube 
6.00000,  3.00000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
1.00000,  1.00000 

0.000000,  90.0000,  10.0000,  90.0000 

CG:  Second  (FORWARD  STARBOARD)  missile  txxbe 
6.00000,  7.00000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
1.00000,  1-00000 

0.000000,  90.0000,  10.0000,  90.0000 


RR:  Make  the  rem.aining  symmetric  tubes.  First  shift 
4.00000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 


CG:  Third  (PORT)  missile  txibe 
6.00000,  3.00000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
1.00000,  1.00000 

0.000000,  90.0000,  10.0000,  90.0000 
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CG:  Forth  (STARBOARD)  missile  tube 
6.00000,  7.00000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
1.00000,  1.00000 

0.000000,  90.0000,  10.0000,  90.0000 

RR:  Make  the  remaining  synmetric  tubes.  Second  shift 
4.00000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 

CG:  Fifth  (PORT)  missile  tube 
6.00000,  3.00000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
1.00000,  1.00000 

0.000000,  90.0000,  10.0000,  90.0000 


CG:  Sixth  (STARBOARD)  missile  tube 
6.00000,  7.00000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
1.00000,  1.00000 

0.000000,  90.0000,  10.0000,  90.0000 

CE:  We  have  to  revert  to  global  coordinates.  AGAIN 


RR:  Back  to  global  coordinates 
-8.00000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0-000000 

CM:  CG:  Test  tube  at  new (0,0,0).  Should  be  placed  at 
origin.  CM:  OK! 

CM:  0.000000,  0.000000,  0.000000 

CM:  0.000000,  0.000000,  90.0000,  0.000000 

CM:  1.00000,  1.00000 

CE:  0.000000,  90.0000,  10.0000,  90.0000 

CE:  FIRST  FLOOR  INTERIOR,  TANKS,  STAIRS,  OPENING 

CM:  Construct  the  TANKS  inside  the  first  floor.  There  will 
CM:  be  three  tanks  along  the  PORT  side  of  the  room.  The 
CM:  FORWARD  and 

CM:  AFT  one  will  have  half  the  height  of  the  room  (2.5m) 
and 

CM:  the  middle  one  will  go  up  to  the  OVERHEAD  (5m)  . 

CM:  Similarly 

CM:  and  syxttmetrically,  will  be  the  STARBOARD  side  of  the 
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CM:  room. 

CM:  In  addition,  one  big  tank  will  be  placed  in  the  AFT 
part 

CE:  with  a  height  of  2.5m. 

CM:  We  CAN  NOT  consider  PORT  tanks,  as  if  there  was  only 
one  CM:  large, 

CM:  along  the  PORT  side  and  an  additional  one  in  the  center 
CM:  because 

CM:  of  the  warnings  V7e  get.  For  the  first, 

CM:  length  20m,  width  0.5m  and  height  2.5m  (half  the  way 
CM:  up)  . 

CM;  We  only  need  to  put  the  STARBOARD  and  UPPER  plate  of 
it, 

CE:  since  the  rest  are  already  put  as  walls  in  the  room. 

PG:  PORT  large  tank,  STARBOARD^-  forward  plate  only 
4,  0 

0.000000,  0.500000,  0.000000 
0.000000,  0.500000,  2.50000 
5.00000,  0.500000,  2.50000 
5.00000,  0.500000,  0.000000 

PG:  PORT  large  tank,  UPPER  -fforward  plate  only 
4,  0 

0.000000,  0.500000,  2.50000 
0.000000,  0.000000,  2.50000 
5.00000,  0.000000,  2.50000 
5.00000,  0.500000,  2.50000 

PG:  PORT  large  tank,  STARBOARD^-  aft  plate  only 
4,  0 

15.0000,  0.500000,  0.000000 
15.0000,  0.500000,  2.50000 
20.0000,  0.500000,  2.50000 
20.0000,  0.500000,  0.000000 

PG:  PORT  large  tank,  UPPER  +aft  plate  only 
4,  0 

15.0000,  0.500000,  2.50000 
15.0000,  0.000000,  2.50000 
20.0000,  0.000000,  2.50000 
20.0000,  0.500000,  2.50000 

PG:  PORT  additional  tank,  FORWARD  plate  only 
4,  0 
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5.00000,  0.000000,  0.000000 
5.00000,  0.000000,  5.00000 
5.00000,  1.00000,  5.00000 
5.00000,  1.00000,  0.000000 

PG:  PORT  additional  tank,  AFT  plate  only 
4 ,  0 

15.0000,  0.000000,  0.000000 
15.0000,  1.00000,  0.000000 
15.0000,  1.00000,  5.00000 
15.0000,  0.000000,  5.00000 

PG:  PORT  additional  tank,  STARBOARD  plate  only 
4,  0 

5.00000,  1.00000,  0.000000 
5.00000,  1.00000,  5,00000 
15.0000,  1.00000,  5.00000 
15.0000,  1.00000,  0.000000 

PG:  PORT  additional  tank,  UPPER  plate  only 
4,  0 

5.00000,  0.000000,  5.00000 
15.0000,  0.000000,  5.00000 
15.0000,  1.00000,  5.00000 
5.00000,  1.00000,  5.00000 

PG:  STARBOARD  large  tank,  PORT+  forward  plate  only 
4,  0 

0.000000,  9.50000,  0.000000 
5.00000,  9.50000,  0.000000 
5.00000,  9.50000,  2.50000 
0.000000,  9.50000,  2.50000 

PG:  STARBOARD  large  tank,  UPPER  and  forward  plate  only 
4,  0 

0.000000,  9.50000,  2.50000 
5.00000,  9.50000,  2.50000 
5.00000,  10.0000,  2.50000 
0.000000,  10.0000,  2.50000 

PG:  STARBOARD  large  tank,  PORT  and  aft  plate  only 
4,  0 

15.0000,  9.50000,  0.000000 
20.0000,  9.50000,  0.000000 
20.0000,  9.50000,  2.50000 
15.0000,  9.50000,  2.50000 
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PG:  STAR30APD  large  tank,  UPPER+  aft  plate  only 
4,  0 

15.0000,  9.50000,  2.50000 
20.0000,  9.50000,  2.50000 
20.0000,  10.0000,  2.50000 
15.0000,  10.0000,  2.50000 

PG:  STARBOARD  additional  tank,  FORWARD  plate  only 
4,  0 

5.00000,  9.00000,  0.000000 
5.00000,  10.0000,  0.000000 
5.00000,  10.0000,  5.00000 
5.00000,  9.00000,  5.00000 

PG:  STARBOARD  additional  tank,  AFT  plate  only 
4,  0 

15.0000,  9.00000,  0.000000 
15.0000,  10.0000,  0.000000 
15.0000,  10.0000,  4.99000 
15.0000,  9.00000,  4.99000 

PG:  STARBOARD  additional  tank,  PORT  plate  only 
4,  0 

5.00000,  9.00000,  0.000000 
15.0000,  9.00000,  0.000000 
15.0000,  9.00000,  5.00000 
5.00000,  9.00000,  5.00000 

PG;  STARBOARD  additional  tank,  UPPER  plate  only 
4,  0 

5.00000,  9.00000,  5.00000 
15.0000,  9.00000,  5.00000 
15.0000,  10.0000,  5.00000 
5.00000,  10.0000,  5.00000 

CM:  We  now  construct  the  AFT  tank,  whose  height  is  2.5m.  We 
only 

CE:  need  the  FORWARD  and  the  UPPER  part 

PG:  AFT  tank  FORWARD  plate 
4,  0 

18.0000,  0.4500000,  0.000000 
18.0000,  0.4500000,  2.50000 
18.0000,  9.550000,  2.50000 
18.0000,  9.550000,  0.000000 
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PG:  AFT  tank  UPPER  plate 
4,  0 

18.0000,  0.4500000,  2.50000 
20.0000,  0.4500000,  2.50000 
20.0000,  9.550000,  2.50000 
18.0000,  9.550000,  2.50000 

CM:  The  Offset  Multi-level (up  and  down)  Ladderway  will  be 
put  in 

CM:  the  center  of  the  room.  First  we  put  a  vertical  plate 
at  x=17 

CE:  and  2=0.5.  This  does  not  coinat  as  a  stair 


PG:  Initial  vertical  plate,  before  first  stair 
4,  0 

17.0000,  4.50000,  0.000000 
17.0000,  4.50000,  0.500000 
17.0000,  5.50000,  0.500000 
17.0000,  5.50000,  0.000000 

PG:  First  Stair  (two  plates  vertical  each  other) . 

Horizontal 

4,  0 

17.0000,  4.50000,  0.500000 
17.0000,  5.50000,  0.500000 

16.7500,  5.50000,  0.500000 

16.7500,  4.50000,  0.500000 

PG:  First  Stair  (two  vertical  plates) .  Vertical 
4,  0 

16.7500,  5.50000,  0.500000 

16.7500,  4.50000,  0.500000 

16.7500,  4.50000,  1.00000 

16.7500,  5.50000,  1.00000 

RR:  Make  the  remaining  symmetric  stairs.  First  shift 
-0.250000,  0.000000,  0.500000 
0.000000,  0.000000,  90.0000,  0.000000 

PG:  Second  Stair  (two  vertical  plates) .  Horizontal 
4,  0 

17.0000,  4.50000,  0.500000 
17.0000,  5.50000,  0.500000 

16.7500,  5.50000,  0.500000 

16.7500,  4.50000,  0.500000 
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PG:  Second  Stair  (two  vertical  plates) .  Vertical 
4,  0 

16.7500,  5.50000,  0.500000 

16.7500,  4.50000,  0.500000 

16.7500,  4.50000,  1.00000 

16.7500,  5.50000,  1.00000 


RR:  Make  the  remaining  symmetric  stairs.  Second  shift 
-0,250000,  0.000000,  0.500000 
0.000000,  0.000000,  90.0000,  0.000000 


PG;  Third  Stair  (two  vertical  plates) .  Horizontal 
4,  0 

17.0000,  4.50000,  0.500000 
17.0000,  5.50000,  0.500000 

16.7500,  5.50000,  0.500000 

16.7500,  4.50000,  0.500000 

PG:  Third  Stair  (two  vertical  plates) .  Vertical 
4,  0 

16.7500,  5.50000,  0.500000 

16.7500,  4.50000,  0.500000 

16.7500,  4.50000,  1-00000 

16.7500,  5.50000,  1.00000 

RR:  Make  the  remaining  symmetric  stairs.  Third  shift 
-0.250000,  0.000000,  0.500000 
0.000000,  0-000000,  90.0000,  0-000000 


PG:  Forth  Stair  (two  vertical  plates) .  Horizontal 
4,  0 

17.0000,  4.50000,  0.500000 
17.0000,  5.50000,  0.500000 

16.7500,  5.50000,  0.500000 

16.7500,  4.50000,  0.500000 

PG:  Forth  Stair  (two  vertical  plates) .  Vertical 
4,  0 

16.7500,  5.50000,  0.500000 

16.7500,  4.50000,  0.500000 

16.7500,  4.50000,  1.00000 

16.7500,  5.50000,  1.00000 

RR;  Make  the  remaining  symmetric  stairs.  Forth  shift 
-0.250000,  0.000000,  0.500000 
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0.000000,  0.000000,  90.0000,  0.000000 


PG:  Fifth  Stair  {two  vertical  plates) .  Horizontal 
4,  0 

17.0000,  4.50000,  0.500000 
17.0000,  5.50000,  0.500000 

16.7500,  5.50000,  0.500000 

16.7500,  4.50000,  0.500000 

PG:  Fifth  Stair  (two  vertical  plates) .  Vertical 
4,  0 

16.7500,  5.50000,  0.500000 

16.7500,  4.50000,  0.500000 

16.7500,  4.50000,  1.00000 

16.7500,  5.50000,  1.00000 


KR:  Make  the  remaining  symmetric  stairs.  Fifth  shift 
-0.250000,  0.000000,  0.500000 
0.000000,  0.000000,  90.0000,  0.000000 

PG:  Sixth  Stair  (two  vertical  plates) .  Horizontal 
4,  0 

17.0000,  4.50000,  0.500000 
17.0000,  5.50000,  0.500000 

16.7500,  5.50000,  0.500000 

16.7500,  4.50000,  0.500000 

PG:  Sixth  Stair  (two  vertical  plates) .  Vertical 
4,  0  ‘ 

16.7500,  5.50000,  0.500000 

16.7500,  4.50000,  0.500000 

16.7500,  4.50000,  1.00000 

16.7500,  5.50000,  1.00000 

RR;  Make  the  remaining  symmetric  stairs.  Sixth  shift 
-0.250000,  0.000000,  0.500000 
0.000000,  0.000000,  90.0000,  0.000000 

PG:  Seventh  Stair  (two  vertical  plates) .  Horizontal 
4,  0 

17.0000,  4.50000,  0.500000 
17.0000,  5.50000,  0.500000 

16.7500,  5.50000,  0.500000 

16.7500,  4.50000,  0.500000 

PG:  Seventh  Stair  (two  vertical  plates) .  Vertical 
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4,  0 

16.7500,  5.50000,  0.500000 

16.7500,  4.50000,  0.500000 

16.7500,  4.50000,  1.00000 

16.7500,  5.50000,  1.00000 

RR:  Make  the  remaining  symmetric  stairs.  Seventh  shift 
-0.250000,  0.000000,  0.500000 
0.000000,  0.000000,  90.0000,  0.000000 

PG:  Eighth  Stair  (two  vertical  plates).  Horizontal 
4,  0 

17.0000,  4.50000,  0.500000 
17.0000,  5.50000,  0.500000 

16.7500,  5.50000,  0.500000 

16.7500,  4.50000,  0.500000 

PG:  Eighth  Stair  (two  vertical  plates) .  Vertical 
4,  0 

16.7500,  5.50000,  0.500000 

16.7500,  4.50000,  0.500000 

16.7500,  4.50000,  1.00000 

16.7500,  5.50000,  1.00000 

RR:  Make  the  remaining  symmetric  stairs.  Eighth  shift 
-0.250000,  0.000000,  0.500000 
0.000000,  0.000000,  90.0000,  0.000000 

PG:  Ninth  Stair  (two  vertical  plates) .  Horizontal 
4,  0 

17.0000,  4.50000,  0.500000 
17.0000,  5.50000,  0.500000 

16.7500,  5.50000,  0.500000 

16.7500,  4.50000,  0.500000 

PG:  Ninth  Stair,  Vertical  0.1cm.  higher  to  avoid  warning 
(p.l61) 

4,  0 

16.7500,  5.50000,  0.500000 

16.7500,  4.50000,  0.500000 

16.7500,  4.50000,  1.00100 

16.7500,  5.50000,  1.00100 

CE :  We  have  to  revert  to  global  coordinates .  AGAIN 


RR;  Back  to  global  coordinates 
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2.00000,  0.000000,  -4.00000 
0.000000,  0.000000,  90.0000,  0.000000 

CM:  CG:  Test  tube  at  new (0,0,0).  Should  be  placed  at 
origin.  CM:  OK! 

CM:  0.000000,  0.000000,  0.000000 

CM:  0.000000,  0.000000,  90.0000,  0.000000 

CM:  1.00000,  1.00000 

CE:  0.000000,  90.0000,  10.0000,  90.0000 

CM:  We  have  already  put  the  fo3rward  vertical  plate  that 
CM :  surrounds 

CM:  the  opening.  We  now  put  the  remaining  three,  trying  to 
CE:  eliminate  the  warnings  (thus  Zmax=5.001) 

PG:  REST  of  the  frame  that  surrotinds  the  opening:  LEFT 
4,  0 

14.7500,  4.50000,  5.00100 

15.7500,  4.50000,  5.00100 

15.7500,  4.50000,  4.70000 

14.7500,  4.50000,  4.70000 

PG:  REST  of  the  frame  that  surro'onds  the  opening:  RIGHT 
4,  0 

14.7500,  5.50000,  5.00100 

14.7500,  5.50000,  4.70000 

15.7500,  5.50000,  4.70000 

15.7500,  5.50000,  5-00100 

PG:  REST  of  the  frame  that  surroiinds  the  opening:  AFT 
4,  0 

15.7500,  4.50000,  5.00100 

15.7500,  5.50000,  5.00100 

15.7500,  5.50000,  4.70000 

15.7500,  4.50000,  4.70000 

CM:  We  should  now  "attach"  a  specific  material  with  certain 
CM:  characteristics  in  the  overhead  of  the  first  floor  in 
CM:  order  to  simulate  the  opening  where  the  stairs  end! 

CM:  Center 

CM:  of  the  opening:  15  5  5  From  x=14->15  y=4.5->5.5  z=5. 

CM:  THIS  UNORTHODOX  METHOD  DOES  NOT  WORK  PROPERLY,  SO  WE 
CM:  CHOOSE 

CM:  TO  WORK  ANALYTICALLY,  BY  CONSTRUCTING  EACH  OTHER  PART 
OF  CM:  THE 
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CE;  OVERHEAD  BUT  THE  OPENING  ITSELF  AND  THEN  "GLUE"  THEM 
TOGETHER! 


CM:  ZD:  Simulated  opening.  Permitivity  I.OICM:  1 
CM:  1 

CM:  0.0500000,  1.01000,  0.000000,  1.00000,  0.000000 
CM: 

CM:  P6:  This  is  an  opening  in  the  overhead 
CM:  4,  0 

CM:  15.7500,  4.50000,  5.00000 
CM:  14.7500,  4.50000,  5.00000 
CM:  14.7500,  5.50000,  5.00000 
CM:  15.7500,  5.50000,  5-00000 
CM: 

CM:  ZD:  Back  again  to  P.E.C. 

CE:  0 


CM:  We  now  t3ry  to  construct  the  HANDRAILS  as  cylinders  v/ith 
CM:  small  radius  and  adequate  angle.  USE  THE  Z  AXIS 
CE:  POSITION-  ROTATE 

CG:  LEFT  HANDRAIL 
17.0000,  4.50000,  0.250000 
335.000,  0.000000,  0.000000,  -999.000 
0.100000,  0.100000 

0.000000,  45.0000,  5.00000,  45.0000 

CG:  RIGHT  HANDRAIL 
17.0000,  5.50000,  0.250000 
335.000,  0.000000,  0.000000,  -999.000 
0.100000,  0.100000 

0.000000,  45.0000,  5.00000,  45.0000 


CM: 


CM: 


END  OF  FIRST  FLOOR' 


CE: 


PL: 

T 

XQ: 


£N:  End 
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CE:  !!3iis  is  ttie  input  file  when  Antenna  is  located  in  a 
CE:  random  place.  We  will  present  only  the  commands 
different  from  the  previous  file. 

CE:  Put  one  antenna  in  (14,5,1). 

SG:  Put  the  source  at  (14,5,1)  position 
14.0000,  5.00000,  1.00000 
0.000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0.000000 
1.00000,  0.000000 

VN:  VOLUMETRIC  PLOT  at  Z=1.5  m. 

0.000000,  0.000000,  0.000000 

0.000000,  0.000000,  90.0000,  0.000000 

T  ‘ 

0.000000,  0.000000,  1.50000 
0.0500000,  0.0500000,  0.000000 
1,  400,  200 


EN:  End 

CE:  This  is  the  input  file  when  two  Antennas  are  located  in 
CE:  thft  room.  We  will  present  only  the  commands  different 
from  the  previous  file 

CM:  Put  2  sources  in  Z=4m  level  and  observation 
CM:  area  at  Z=lm.  Rear  tank  plates  overlap  to  avoid 
CE:  diffraction  inside  the  tank 


SG:  Put  the  source  at  (14,5,4)  position 
14.0000,  5.00000,  4.00000 
0.000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0.000000 
1.00000,  0.000000 


SG:  Put  ANOTHER  source  at  (6,5,4)  position 
6.0000,  5.00000,  4.00000 
0.000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0.000000 
1.00000,  0.000000 
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VN:  VOLUMETRIC  PLOT  at  Z=1 . 0  m. 

0.000000,  0-000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
T 

0.000000,  0.000000,  i. 00000 
0.0500000,  0.0500000,  0.000000 
1,  400,  200 

EN:  End 

CE:  This  is  the  input  file  for  room  to  room  case.  This  is 
a  different  file.  We  will  present  only  the  commands 
different  from  the  previous  file. 

ZD;  P.S.C. 

0 

FR:  We  will  use  a  frequency  of  2.4  GHz  (ISM  band} 

2.40000 

US:  (Necessary)  in  Wavelengths 
0 

SG:  Put  the  source  at  (14,5,1)  position 
14.0000,  5.00000,  1.00000 
0.000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0.000000 
1.00000,  0.000000 

CM:  Antenna  is  located  in  the  first  floor,  but 
CM:  observervation 

CM:  points  are  in  the  second  floor  (1.5m  above  the  floor) . 
CM:  Recall  that  the  problem  with  the  plate  edges 
diffraction 
CE;  still  exists 

VN:  VOLUMETRIC  PLOT  at  Z=6.5  m 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0-000000 
T 

1.00000,  1.00000,  6.50000 
0.0500000,  0.0500000,  0.000000 


131 


1,  400,  200 


CE:  Construct  the  ROOM,  according  to  the  RIGHT-HAMD  RULE 
(SAME  AS  BEFORE) 

CM:  - 

CM;  - END  OF  FIRST  FLOOR - 

CE:  - 


CE:  SECOND  FLOOR  (MISSILES  AND  FRAME  ALREADY  PUT) 

CM:  SECOND  OVERHEAD:  This  is  an  Overhead  at  z=10  level.  We 
CM:  will 

CE:  use  the  first  construction,  shifted  by  Sm  along  Z  axis 

CM:  PG:  level:  (0,0,10)  - SECOND  OVERHEAD - 

CM:  Theoretically - 

CM;  4,  0 

CM:  0.000000,  0.000000,  10.00000 
CM:  0.000000,  10.0000,  10.00000 
CM:  20.0000,  10.0000,  10.00000 
CE:  20.0000,  0.000000,  10.00000 

CE:  Construct  all  the  plates,  but  the  opening  itself 


PG:  RIGHT 
4,  0 

0.000000,  5.50000,  10.00000 
0.000000,  10.0000,  10.00000 
20.0000,  10.0000,  10.00000 
20.0000,  5.50000,  10.00000 

PG:  MIDDLE 
4,  0 

15.7500,  5.51000,  10.00000 
20.0000,  5.51000,  10.00000 
20.0000,  4.49000,  10.00000 

15.7500,  4.49000,  10.00000 

PG:  MIDDLE 
4,  0 

0.000000,  5.51000,  10.00000 

14.7500,  5.51000,  10.00000 
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14.7500,  4.49000,  10.00000 
0.000000,  4.49000,  10.00000 

PG:  LEFT 
4,  0 

0.000000,  0.000000,  10.00000 
0.000000,  4.50000,  10.00000 
20.0000,  4.50000,  10.00000 
20.0000,  0.000000,  10.00000 


CM:  In  case  we  want  the  overhead  without  the 
OPENING ( solid) , 

CE:  we  xmcomment  the  following  command: 

CM:  PG:  This  is  an  OPENING  in  the  second  overhead 
CM:  4,  0 

CM:  15.7500,  4.50000,  10.00000 
CM:  14.7500,  4.50000,  10.00000 
CM:  14.7500,  5.50000,  10.00000 
CE:  15.7500,  5.50000,  10.00000 

CM:  The  SECOND  Offset  Multi-level (up  and  down)  Ladderway 
CM:  will  be  put  in 

CM:  the  aft  part  of  the  room.  First  we  put  a  vertical  plate 
CM:  at  x=17 

CE:  and  z=0.5.  This  does  not  count  as  a  stair 

RR;  Use  the  first  construction  for  the  SECOND  LADDERWAY 

0.000000,  0.000000,  5.00000 

0.000000,  0.000000,  90.0000,  0.000000 

CM:  PG:  Initial  vertical  plate,  before  first  stair  Warnings 
CM:  4,  0 

CM:  17.0000,  4.50000,  0.000000 
CM:  17.0000,  4.50000,  0.500000 
CM:  17.0000,  5.50000,  0.500000 
CE:  17.0000,  5.50000,  0.000000 

PG:  First  Stair  (two  plates  vertical  each  other) .Horizontal 
4,  0 

17.0000,  4.50000,  0.500000 
17.0000,  5.50000,  0.500000 

16.7500,  5.50000,  0.500000 

16.7500,  4.50000,  0.500000 
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PG:  First  Stair  {two  vertical  plates) .  Vertical 

4,  0 

16.7500,  5.50000,  0.500000 

16.7500,  4.50000,  0.500000 

16.7500,  4.50000,  1.00000 

16.7500,  5.50000,  1.00000 

RR:  Make  the  remaining  symmetric  stairs.  First  shift 
-0.250000,  0.000000,  0.500000 
0.000000,  0.000000,  90.0000,  0.000000 

PG:  Second  Stair  (two  vertical  plates) .  Horizontal 
4,  0 

17.0000,  4.50000,  0.500000 
17.0000,  5.50000,  0.500000 

16.7500,  5.50000,  0.500000 

16.7500,  4.50000,  0.500000 

PG:  Second  Stair  {two  vertical  plates) .  Vertical 
4,  0 

16.7500,  5.50000,  0.500000 

16.7500,  4.50000,  0.500000 

16.7500,  4.50000,  1.00000 

16.7500,  5.50000,  1.00000 

RR:  Make  the  remaining  symmetric  stairs.  Second  shift 
-0.250000,  0.000000,  0-500000 
0.000000,  0.000000,  90-0000,  0.000000 

PG:  Third  Stair  {two  vertical  plates) .  Horizontal 
4,  0 

17.0000,  4.50000,  0.500000 
17.0000,  5.50000,  0.500000 

16.7500,  5.50000,  0.500000 

16.7500,  4.50000,  0.500000 

PG:  Third  Stair  (two  vertical  plates) .  Vertical 
4,  0 

16.7500,  5.50000,  0.500000 

16.7500,  4.50000,  0.500000 

16.7500,  4.50000,  1.00000 

16.7500,  5.50000,  1.00000 

RR:  Make  the  remaining  symmetric  stairs.  Third  shift 
-0.250000,  0.000000,  0.500000 
0.000000,  0.000000,  90.0000,  0.000000 


PG;  Forth  Stair  (two  vertical  plates) .  Horizontal 
4,  0 

17.0000,  4.50000,  0.500000 
17.0000,  5.50000,  0.500000 

16.7500,  5.50000,  0.500000 

16.7500,  4.50000,  0.500000 

PG:  Forth  Stair  (two  vertical  plates) .  Vertical 
4,  0 

16.7500,  5.50000,  0.500000 

16.7500,  4.50000,  0.500000 

16.7500,  4.50000,  1.00000 

16.7500,  5.50000,  1.00000 

RR:  Make  the  remaining  symmetric  stairs.  Forth  shift 
-0.250000,  0.000000,  0.500000 
0.000000,  0.000000,  90.0000,  0.000000 


PG:  Fifth  Stair  (two  vertical  plates) .  Horizontal 
4,  0 

17.0000,  4.50000,  0.500000 
17.0000,  5.50000,  0.500000 

16.7500,  5.50000,  0.500000 

16.7500,  4.50000,  0.500000 

PG:  Fifth  Stair  (two  vertical  plates) .  Vertical 
4,  0 

16.7500,  5.50000,  0.500000 

16.7500,  4.50000,  0.500000 

16.7500,  4.50000,  1.00000 

16.7500,  5.50000,  1.00000 

RR:  Make  the  remaining  symmetric  stairs.  Fifth  shift 
-0.250000,  0.000000,  0.500000 
0.000000,  0.000000,  90.0000,  0.000000 

PG:  Sixth  Stair  (two  vertical  plates) .  Horizontal 
4,  0 

17.0000,  4.50000,  0.500000 
17.0000,  5.50000,  0.500000 

16.7500,  5.50000,  0.500000 

16.7500,  4.50000,  0.500000 


PG:  Sixth  Stair  (two  vertical  plates) .  Vertical 
4,  0 
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16.7500,  5.50000,  0.500000 

16.7500,  4.50000,  0.500000 

16.7500,  4.50000,  1.00000 

16.7500,  5.50000,  1.00000 

RR:  Make  the  remaining  symmetric  stairs.  Sixth  shift 
-0.250000,  0.000000,  0.500000 
0.000000,  0.000000,  90.0000,  0.000000 


PG :  Seventh  Stair 
4,  0 

(two  vertical 

plates) . 

17.0000,  4.50000, 

0.500000 

17.0000,  5.50000, 

0.500000 

16.7500,  5.50000, 

0.500000 

16.7500,  4.50000, 

0.500000 

PG:  Seventh  Stair 
4,  0 

(two  vertical 

plates) . 

16.7500,  5.50000, 

0.500000 

15.7500,  4.50000, 

0.500000 

16.7500,  4.50000, 

1.00000 

16.7500,  5.50000, 

1.00000 

Horizontal 


Vertical 


RR:  Make  the  remaining  symmetric  stairs.  Seventh  shift 
-0.250000,  0.000000,  0.500000 
0.000000,  0.000000,  90.0000,  0.000000 

PG:  Eighth  Stair  (two  vertical  plates) .  Horizontal 
4,  0 

17.0000,  4.50000,  0.500000 
17.0000,  5.50000,  0.500000 

16.7500,  5.50000,  0.500000 

16.7500,  4.50000,  0.500000 

PG:  Eighth  Stair  (two  vertical  plates) .  Vertical 
4,  0 

16.7500,  5.50000,  0.500000 

16.7500,  4.50000,  0.500000 

16.7500,  4.50000,  1.00000 

16.7500,  5.50000,  1.00000 

RR:  Make  the  remaining  symmetric  stairs.  Eighth  shift 
-0.250000,  0.000000,  0.500000 
0.000000,  0.000000,  90.0000,  0.000000 


PG:  Ninth  Stair  (two  vertical  plates) .  Horizontal 
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4,  0 

17.0000,  4.50000,  0.500000 
17.0000,  5.50000,  0.500000 

16.7500,  5.50000,  0.500000 

16.7500,  4.50000,  0.500000 

PG:  Ninth  Stair,  Vertical.  O.lcin  higher  to  avoid  warnings 
{p.l61) 

4,  0 

16.7500,  5.50000,  0.500000 

16.7500,  4.50000,  0.500000 

16.7500,  4.50000,  1.00100 

16.7500,  5-50000,  1.00100 

CE:  We  have  to  revert  to  global  coordinates.  AGAIN 


RR:  Back  to  global  coordinates.  Z  TRANSLATIONS  HAS  CHANGED 
2.00000,  0.000000,  -9.00000 
0.000000,  0.000000,  90.0000,  0.000000 


CM:  CG:  Test  tiobe  at  new (0,0,0).  Should  be  placed  at 
origin.  CM:  OK! 

CM:  0.000000,  0.000000,  0.000000 

CM:  0.000000,  0.000000,  90.0000,  0.000000 

CM:  1.00000,  1.00000 

CE:  0.000000,  90.0000,  10.0000,  90.0000 


CM:  We  have  already  put  the  forward  vertical  plate  that 
surroxmds 

CM:  the  opening.  We  now  put  the  reinaining  three,  trying  to 
CE:  eliminate  the  warnings  (thus  Zmax=5.001) 


PG:  REST  of  the  frame  that  surrounds  the  opening:  LEFT 
4,  0 

14.7500,  4.50000,  10.00100 

15.7500,  4.50000,  10.00100 

15.7500,  4.50000,  9.70000 

14.7500,  4.50000,  9.70000 


PG:  REST  of  the  frame  that 
4,  0 

14.7500,  5.50000,  10.00100 

14.7500,  5.50000,  9.70000 

15.7500,  5.50000,  9.70000 

15.7500,  5.50000,  10.00100 


surrounds  the  opening:  RIGHT 
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PG:  REST  of  the  frame  that  surrounds  the  opening:  AFT 
4,  0 

15.7500,  4.50000,  10.00100 

15.7500,  5.50000,  10.00100 

15.7500,  5.50000,  9.70000 

15.7500,  4.50000,  9.70000 

CM:  ZD:  Simulated  opening.  Penaitivity  1.01 
CM:  1 
CM:  1 

CM:  0.0500000,  1.01000,  0.000000,  1.00000,  0.000000 
CM: 

CM:  PG:  This  is  an  opening  in  the  overhead 
CM:  4,  0 

CM:  15.7500,  4.50000,  5.00000 
CM:  14.7500,  4.50000,  5.00000 
CM:  14.7500,  5.50000,  5.00000 
CM:  15.7500,  5.50000,  5.00000 
CM: 

CM:  ZD:  Back  again  to  P.E.C. 

CE:  0 


CM:  We  now  try  to  construct  the  HANDRAILS  as  cylinders  with 
CM:  small  radius  and  adequate  angle.  USE  THE  Z  AXIS 
CE:  POSITION-  ROTATE 

RR:  Use  the  first  construction  for  the  SECOND  LADDERWAY 
0.000000,  0.000000,  5.00000 
0.000000,  0.000000,  90.0000,  0.000000 

CG:  LEFT  HANDRAIL 
17.0000,  4.50000,  0.250000 
335.000,  0.000000,  0.000000,  -999.000 
0.100000,  0.100000 

0.000000,  45.0000,  5.00000,  45.0000 

CG:  RIGHT  HANDRAIL 
17.0000,  5.50000,  0.250000 
335.000,  0.000000,  0.000000,  -999.000 
0.100000,  0-100000 

0.000000,  45.0000,  5.00000,  45.0000 


CE :  We  have  to  revert  to  global  coordinates .  AGAIN 

RR:  Back  to  global  coordinates.  Z  GOES  5m  DOWN: 

0.00000,  0.000000,  -5.00000 
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0-000000,  0.000000,  90.0000,  0.000000 

CM:  CG:  Test  tube  at  (0,0,0).  Should  be  placed  at  origin. 
CM:  OK! 

CM:  0.000000,  0.000000,  0.000000 

CM:  0.000000,  0.000000,  90.0000,  0.000000 

CM:  1.00000,  1.00000 

CE:  0.000000,  90.0000,  10.0000,  90.0000 

LY: 

T 

CM:  LP:  See  data  in  the  out  file  the  same  with  PL  command 
CE:  T 

PL: 

T 

XQ: 

EN:  End 


The  following,  are  the  input  sub-files  for  the  truncation 
of  the  CE:  room.  Each  room,  as  well  as  the  free  space 
volume,  is  truncated  into  several  CELLS,  in  order  to  avoid 
observation  points  inside  the  numerous  obstacles.  For  the 
shake  of  sin5>licity,  (we  would  need  at  least  700  pages  for 
all  the  files),  we  will  only  present  the  different  commands 
for  antenna  and  observation  points  placements. 

ROOM  DECOMPOSITION  INTO  CELLS  (FIRST  LOCATION) 

CELL  1 

SG:  Put  the  source  at  (14,5,1)  position 
14.0000,  5.00000,  1.00000 
0.000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0.000000 
1.00000,  0.000000 

VK:  VOLUMETRIC  PLOT  1  at  Z=1 . 5  m. 

0.000000,  0.000000,  0.000000 

0.000000,  0.000000,  90.0000,  0.000000 

T 

0.100000,  1.20000,  1.50000 


139 


0.250000,  0.250000,  0.000000 
.1,  71,  3 

CELL  2 

SG:  Put  the  source  at  (14,5,1 
14.0000,  5.00000,  1.00000 
0.000000,  0.000000,  90.0000, 
-1,  0.500000,  0.000000 
1.00000,  0.000000 

VN:  VOLUMETRIC  PLOT  1  at  Z=l. 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000, 

T 

0.100000,  2.20000,  1.50000 
0.250000,  0.250000,  0.000000 
1,  19,  7 

CELL  3 

SG:  Put  the  source  at  (14,5,1 
14.0000,  5.00000,  1.00000 
0.000000,  0.000000,  90.0000, 
-1,  0.500000,  0.000000 
1.00000,  0.000000 

VN:  VOLUMETRIC  PLOT  1  at  Z=l. 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000, 

T 

7.10000,  2.20000,  1.50000 
0.250000,  0.250000,  0.000000 

I,  7,  7 

CELL  4 

SG:  Put  the  source  at  (14,5,1) 
14.0000,  5.00000,  1.00000 
0.000000,  0.000000,  90.0000, 
-1,  0.500000,  0.000000 
1.00000,  0.000000 

VN:  VOLUMETRIC  PLOT  1  at  Z=l. 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000, 

T 

II. 1000,  2.20000,  1.50000 
0.250000,  0.250000,  0.000000 


position 

.000000 

m. 

.000000 


position 

.000000 

m. 

.000000 


position 

.000000 

m. 

.000000 
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1,  7,  7 

CELL  5 

SG:  Put  the  source  at  (14,5,1)  position 
14.0000,  5.00000,  1.00000 
0.000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0.000000 
1.00000,  0.000000 

VK:  VOLUMETRIC  PLOT  1  at  2=1.5  m. 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
T 

15-1000,  2.20000,  1.50000 
0.250000,  0.250000,  0.000000 
1,  11,  7 

CELL  6 

SG:  Put  the  source  at  (14,5,1)  position 
14.0000,  5.00000,  1.00000 
0.000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0.000000 
1.00000,  0.000000 

VN:  VOLUMETRIC  PLOT  1  at  2=1.5  m. 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  o.tooooo 

T 

0.100000,  4.20000,  1.50000 
0.250000,  0.250000,  0.000000 
1,  71,  7 

CELL  7 

SG:  Put  the  source  at  (14,5,1)  position 
14.0000,  5.00000,  1.00000 
0.000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0.000000 
1.00000,  0.000000 


VN:  VOLUMETRIC  PLOT  1  at  2=1.5  m. 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
T 

0.100000,  6.20000,  1.50000 
0.250000,  0.250000,  0.000000 
1,  19,  7 


141 


CELL  8 

SG:  Put  the  source  at  {14,5,1)  position 
14.0000,  5.00000,  1.00000 
0.000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0.000000 
1.00000,  0.000000 


VN:  VOLUMETRIC  PLOT  1  at  Z=1.5  m. 

0.000000,  0.000000,  0.000000  ' 

0.000000,  0.000000,  90.0000,  0.000000 
T 

7.10000,  6.20000,  1.50000 
0.250000,  0.250000,  0.000000 

I,  7,  7 

CELL  9 

SG:  Put  the  source  at  (14,5,1)  position 
14.0000,  5.00000,  1.00000 
0.000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0-000000 
1.00000,  0.000000 

VM:  VOLUMETRIC  PLOT  1  at  2=1.5  m. 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
T 

II. 1000,  6.20000,  1-50000 
0.250000,  0.250000,  0.000000 
1,  7,  7 

CELL  10 

SG:  Put  the  source  at  (14,5,1)  position 
14.0000,  5.00000,  1.00000 
0.000000,  0.000000,  90-0000,  0.000000 
-1,  0.500000,  0.000000 
1.00000,  0.000000 

VN:  VOLUMETRIC  PLOT  1  at  2=1.5  m. 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
T 

15.1000,  6.20000,  1.50000 
0.250000,  0.250000,  0.000000 
1,  11,  7 
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CELL  11 

SG:  Put  the  source  at  (14,5,1)  position 
14.0000,  5.00000,  1.00000 
0.000000,  0-000000,  90.0000,  0.000000 
-1,  0.500000,  0.000000 
1.00000,  0.000000 

VN:  VOLUMETRIC  PLOT  1  at  Z=1.5  m. 

0.000000,  0.000000,  0.000000 

0.000000,  0.000000,  90.0000,  0.000000 

T 

0.100000,  8.20000,  1.50000 
0.250000,  0.250000,  0.000000 

1,  71,  3 

FREE  SPACE 
CELL  1 


SG:  Put  the  source  at  (14,5,1)  position 
14.0000,  5.00000,  1.00000 
0.000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0.000000 
1.00000,  0.000000 

VN:  VOLUMETRIC  PLOT  1  at  Z=1.5  m. 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
T 

0.100000,  1.20000,  1.50000 
0.250000,  0.250000,  0.000000 
1,  71,  3 


CELL  2 

SG:  Put  the  source  at  (14,5,1)  position 
14.0000,  5.00000,  1.00000 
0.000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0.000000 
1.00000,  0.000000 

VN:  VOLUMETRIC  PLOT  1  at  Z=1.5  m. 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
T 

0.100000,  2.20000,  1.50000 


143 


0.250000,  0.250000,  0.000000 
1,  19,  7 


LY: 


CELL  3 

SG:  Put  the  source  at  (14,5, 
14.0000,  5.00000,  1.00000 
0.000000,  0.000000,  90.0000, 
-1,  0.500000,  0.000000 
1.00000,  0.000000 


VN:  VOLUMETRIC  PLOT  1  at  Z=1 . 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000, 

T 

7.10000,  2.20000,  1.50000 
0.250000,  0.250000,  0.000000 
1,  7,  7 


CELL  4 


SG:  Put  the  source  at  (14,5, 
14.0000,  5.00000,  1.00000 
0.000000,  0.000000,  90.0000, 
-1,  0.500000,  0.000000 
1.00000,  0.000000 


VN:  VOLUMETRIC  PLOT  1  at  Z=l. 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000, 

T 

11.1000,  2.20000,  1.50000 
0.250000,  0.250000,  0.000000 
1,  7,  7 


CELL  5 

SG:  Put  the  source  at  (14,5,1) 
14.0000,  5.00000,  1.00000 
0.000000,  0.000000,  90.0000,  C 


position 

.000000 

m. 

.000000 


position 

.000000 

m. 

.000000 


position 

.000000 
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-1,  0.500000,  0.000000 

1.00000,  0.000000 


VN:  VOLUMETRIC  PLOT  1  at  Z=1.5  m. 

0.000000,  0.000000,  0.000000 

0.000000,  0.000000,  90.0000,  0.000000 

T 

15-1000,  2.20000,  1.50000 
0.250000,  0.250000,  0.000000 

1,  11,  7 


CELL  6 

SG:  Put  the  source  at  (14,5,1}  position 
14.0000,  5.00000,  1.00000 
0.000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0.000000 
1.00000,  0.000000 


VN:  VOLUMETRIC  PLOT  1  at  Z=1.5  m. 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
T 

0.100000,  4.20000,  1.50000 
0.250000,  0.250000,  0.000000 
1,  71,  7 


EN 

CELL  7 

SG:  Put  the  source  at  (14,5,1)  position 
14.0000,  5.00000,  1.00000 
0.000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0.000000 
1.00000,  0.000000 

VN:  VOLUMETRIC  PLOT  1  at  Z=1 . 5  m. 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
T 

0.100000,  6.20000,  1.50000 
0.250000,  0.250000,  0.000000 
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1,  19,  7 


CELL  8 


SG:  Put  the  source  at  (14,5,1)  position 
14.0000,  5.00000,  1.00000 
0-000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0.000000 
1.00000,  0.000000 

W:  VOLUMETRIC  PLOT  1  at  Z=1.5  m. 

0.000000,  0.000000,  0.000000 

0.000000,  0.000000,  90.0000,  0.000000 

T 

7.10000,  6.20000,  1.50000 
0.250000,  0.250000,  0.000000 

1,  7,  7 


CELL  9 

SG:  Put  the  source  at  (14,5,1)  position 
14.0000,  5.00000,  1.00000 
0.000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0.000000 
1.00000,  0.000000 

VN:  VOLUMETRIC  PLOT  1  at  Z=1.5  m. 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0,000000 
T 

11.1000,  6.20000,  1.50000 
0.250000,  0.250000,  0.000000 
1,  7,  7 


CELL  10 


SG:  Put  the  source  at  (14,5,1)  position 
14.0000,  5.00000,  1.00000 
0.000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0.000000 
1.00000,  0.000000 

VN:  VOLUMETRIC  PLOT  1  at  Z=1.5  m. 
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0.000000,  0.000000,  0.000000 
0-000000,  0.000000,  90.0000,  0.000000 
T 

15.1000,  6-20000,  1.50000 
0.250000,  0.250000,  0.000000 
1,  11,  7 


CELL  11 

SG:  Put  the  source  at  (14,5,1)  position 
14.0000,  5.00000,  1.00000 
0.000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0.000000 
1.00000,  0.000000 

VN:  VOLUMETRIC  PLOT  1  at  Z=1.5  m. 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
T 

0.100000,  8.20000,  1.50000 
0.250000,  0.250000,  0.000000 
1,  71,  3 


ROOM  DECOMPOSITION  INTO  CELLS  (SECOND  LOCATION) 
CELL  1 

SG:  Put  the  source  at  (2,2,2)  position 
2.00000,  2.00000,  2.00000 
0.000000,  0.000000,  90.0000,  0.000000 
-1, ,  0.500000,  0.000000 
1.00000,  0.000000 

VN:  VOLUMETRIC  PLOT  1  at  Z=2  m. 

0-000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
T 

0.100000,  1.20000,  2.00000 
0.250000,  0.250000,  0.000000 
1,  71,  3 


CELL  2 

SG:  Put  the  source  at  (2,2,2)  position 
2.00000,  2.00000,  2.00000 
0.000000,  0.000000,  90.0000,  0.000000 
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-1,  0.500000,  0.000000 

1.00000,  0.000000 


VN:  VOLUMETRIC  PLOT  1.  at  Z=2  m. 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
T 

0.100000,  2.20000,  2.00000 
0.250000,  0.250000,  0.000000 
1,  19,  7 

CELL  3 

SG:  Put  the  source  at  (2,2,2)  position 
2.00000,  2.00000,  2.00000 
0.000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0.000000 
1.00000,  0.000000 


VN:  VOLU5IETRIC  PLOT  1  at  Z=2m 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
T 

7.10000,  2.20000,  2.00000 
0.250000,  0.250000,  0.000000 
1,  7,  7 

CELL  4 

SG:  Put  the  source  at  (2,2,2)  position 
2.00000,  2.00000,  2.00000 
0.000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0.000000 
1.00000,  0.000000 


VN:  VOLUMETRIC  PLOT  1  at  Z=2m. 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
T 

11.1000,  2.20000,  2.00000 
0.250000,  0.250000,  0.000000 
1,  7,  7 

CELL  5 

SG:  Put  the  source  at  (2,2,2)  position 

2.00000,  2.00000,  2.00000 
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0.000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0.000000 
1.00000,  0.000000 

VN:  VOLUMETRIC  PLOT  1  at  Z=2  m. 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
T 

15.1000,  2.20000,  2.00000 
0.250000,  0.250000,  0.000000 
1,  11,  7 

CELL  6 

SG:  Put  the  source  at  (2,2,2)  position 
2.00000,  2.00000,  2.00000 
0.000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0.000000 
1.00000,  0.000000 

VN:  VOLUMETRIC  PLOT  1  at  Z=2m. 

0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
T 

0.100000,  4.20000,  2.00000 
0.250000,  0.250000,  0.000000 
1,  71,  7 

CELL  7 

SG:  Put  the  source  at  (2,2,2)  position 
2.00000,  2.00000,  2.00000 
0.000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0.000000 
1.00000,  0.000000 

VN:  VOLUMETRIC  PLOT  1  at  Z=2m. 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
T 

0.100000,  6.20000,  2.00000 
0.250000,  0.250000,  0.000000 
1,  19,  7 

CELL  8 

SG:  Put  the  source  at  (2,2,2)  position 

2.00000,  2.00000,  2.00000 

0.000000,  0.000000,  90.0000,  0.000000 
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-1,  0.500000,  0.000000 

1.00000,  0.000000 


VN:  VOLUMETRIC  PLOT  1  at  Z=2  in. 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
T 

7.10000,  6.20000,  2.00000 
0.250000,  0.250000,  0-000000 
1,  7,  7 

CELL  9 

SG:  Put  the  source  at  (2,2,2)  position 
2.00000,  2.00000,  2.00000 
0.000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0.000000 
1.00000,  0.000000 


VN:  VOLUMETRIC  PLOT  1  at  Z=2m. 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
T 

11.1000,  6.20000,  2.00000 
0.250000,  0.250000,  0.000000 
1,  7,  7 

CELL  10 

SG:  Put  the  source  at  (2,2,2)  position 
2.00000,  2.00000,  2.00000 
0.000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0.000000 
1.00000,  0.000000 

VN:  VOLUMETRIC  PLOT  1  at  Z=2  m. 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
T 

15.1000,  6.20000,  2.00000 
0.250000,  0.250000,  0.000000 
1,  11,  7 

CELL  11 

SG:  Put  the  source  at  (2,2,2)  position 
2.00000,  2.00000,  2.00000 
0.000000,  0.000000,  90.0000,  0.000000 
-1,  0-500000,  0.000000 
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1.00000,  0.000000 


VN:  VOLUMETRIC  PLOT  1  at  Z=2ra 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  I 
T 

0.100000,  8.20000,  2.00000 
0.250000,  0.250000,  0.000000 
1,  71,  3 

FREE  SPACE 
CELL  1 

SG:  Put  the  source  at  (2,2,2) 


2.00000, 

2.00000,  2 

.00000 

0.000000, 

0.000000, 

90.0000,  1 

-1,  0.500000,  0.000000 

1.00000, 

0.000000 

VN:  VOLUMETRIC  PLOT 

1  at  Z=2m 

0.000000, 

0.000000, 

0.000000 

0.000000, 

T 

0.000000, 

90.0000,  1 

0.100000, 

1.20000, 

2.00000 

0.250000, 
1,  71,  3 

0.250000, 

0.000000 

CELL  2 

SG:  Put  the  source 

at  (2,2,2) 

2.00000, 

2.00000,  2 

.00000 

0.000000, 

0.000000, 

90.0000, 

-1,  0.500000,  0.000000 

1.00000, 

0.000000 

VN:  VOLUMETRIC  PLOT 

1  at  Z=2in 

0.000000, 

0.000000, 

0.000000 

0.000000, 

T 

0.000000, 

90.0000, 

0.100000, 

2.20000, 

2.00000 

0.250000, 
1,  19,  7 

0.250000, 

0.000000 

CELL  3 


.000000 


position 

1.000000 

!. 000000 


position 

1.000000 

1.000000 
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SG:  Put  the  source  at  (2,2,2)  position 
2-00000,  2.00000,  2.00000 
0.000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0.000000 
1.00000,  0.000000 

VN:  VOLUMETRIC  PLOT  1  at  Z=2m. 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
T 

7.10000,  2.20000,  2.00000 
0.250000,  0.250000,  0.000000 
1,  7,  7 


CELL  4 


SG:  Put  the  source  at  (2,2,2)  position 
2.00000,  2.00000,  2.00000 
0.000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0.000000 
1.00000,  0.000000 


VN:  VOLUMETRIC  PLOT  1  at  Z=2  m. 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
T 

11.1000,  2.20000,  2.00000 
0.250000,  0.250000,  0.000000 
T,  7,  7 


CELL  5 


SG:  Put  the  source  at  (2,2,2)  position 
2.00000,  2.00000,  2.00000 
0.000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0.000000 
1.00000,  0.000000 

VH:  VOLUMETRIC  PLOT  1  at  Z=2m. 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
T 
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15.1000,  2.20000,  2.00000 
0.250000,  0.250000,  0.000000 
1,  11,  7 


CELL  6 

SG :  Put  the  source  at  (2,2,2) 
2.00000,  2.00000,  2.00000 
0.000000,  0.000000,  90.0000,  i 
-1,  0.500000,  0.000000 
1.00000,  0.000000 

VN:  VOLUMETRIC  PLOT  1  at  Z=2in 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  I 
T 

0.100000,  4.20000,  2.00000 
0.250000,  0.250000,  0.000000 
1,  71,  7 


CELL  7 

SG:  Put  the  source  at  (2,2,2) 
2.00000,  2.00000,  2.00000 
0.000000,  0.000000,  90.0000,  ( 
-1,  0.500000,  0.000000 
1.00000,  0.000000 

VN:  VOLUMETRIC  PLOT  1  at  Z=2in 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  i 
T 

0.100000,  6.20000,  2.00000 
0.250000,  0.250000,  0.000000 
1,  19,  7 


CELL  8 

SG:  Put  the  source  at  (2,2,2) 
2.00000,  2.00000,  2.00000 
0.000000,  0.000000,  90.0000, 
-1,  0.500000,  0.000000 


position 

1.000000 

1.000000 


position 

1.000000 

1.000000 


position 

1.000000 
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1.00000,  0.000000 


VN:  VOLUMETRIC  PLOT  1  at  Z=2m. 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
T 

7.10000,  6.20000,  2.00000 
0.250000,  0.250000,  0.000000 
1,  7,  7 


CELL  9 


SG:  Put  the  source  at  (2,2,2)  position 
2.00000,  2.00000,  2.00000 
0.000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0.000000 
1.00000,  0.000000 


VN:  VOLUMETRIC  PLOT  1  at  Z=2m. 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
T 

11.1000,  6.20000,  2.00000 
0.250000,  0.250000,  0.000000 
1,  7,  7 


CELL  10 


SG:  Put  the  source  at  (2,2,2)  position 
2.00000,  2.00000,  2.00000 
0.000000,  0.000000,  90.0000,  0-000000 
-1,  0.500000,  0.000000 
1.00000,  0.000000 

VN:  VOLUMETRIC  PLOT  1  at  Z=2m. 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
T 

15.1000,  6.20000,  2.00000 
0.250000,  0.250000,  0.000000 
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1,  11,  7 


CELL  11 

SG:  Put  the  source  at  (2,2,2)  position 
2.00000,  2.00000,  2.00000 
0.000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0.000000 
1.00000,  0.000000 

VE:  VOLUMETRIC  PLOT  1  at  Z=2  m. 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
T 

0.100000,  8.20000,  2.00000 
0.250000,  0.250000,  0.000000 
1,  71,  3 


ROOM  DECOMPOSITION  INTO  CELLS  (THIRD  LOCATION) 
CELL  1 

SG:  Put  the  source  at  {2, 2, 2. 2)  position 
2.00000,  2.00000,  2.20000 
0.000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0.000000 
1.00000,  0.000000 

VN:  VOLUMETRIC  PLOT  1  at  Z=2.2  rti. 

0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
T 

0.100000,  1.20000,  2.200000 
0.250000,  0.250000,  0.000000 
1,  71,  3 


CELL  2 

SG:  Put  the  source  at  (2, 2, 2. 2)  position 

2.00000,  2.00000,  2.20000 

0.000000,  0.000000,  90.0000,  0.000000 

-1,  0.500000,  0.000000 

1.00000,  0.000000 
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VN:  VOLUMETRIC  PLOT  1  at  Z=2.2  iti. 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
T 

0.100000,  2.20000,  2.200000 
0.250000,  0.250000,  0.000000 
1,  19,  7 

CELL  3 

SG:  Put  the  source  at  (2, 2, 2. 2)  position 
2.00000,  2.00000,  2.20000 
0.000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0.000000 
1.00000,  0.000000 

VN:  VOLUMETRIC  PLOT,  1  at  Z=2.2m 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
T 

7.10000,  2.20000,  2.200000 
0.250000,  0.250000,  0.000000 
1,  7,  7 

CELL  4 

SG:  Put  the  source  at  (2, 2, 2. 2)  position 
2.00000,  2.00000,  2.20000 
0.000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0.000000 
1.00000,  0.000000 


VN:  VOLUMETRIC  PLOT  1  at  Z=2.2m. 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
T 

11.1000,  2.20000,  2.200000 
0.250000,  0.250000,  0.000000 
1,  7,  7 

CELL  5 

SG:  Put  the  source  at  {2, 2, 2. 2)  position 
2.00000,  2.00000,  2.20000 
0.000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0.000000 
1.00000,  0.000000 
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VN:  VOLUMETRIC  PLOT  1  at  Z=2.2  m. 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
T 

15.1000,  2.20000,  2.200000 
0.250000,  0.250000,  0.000000 
1,  11,  7 

CELL  6 

SG:  Put  the  source  at  {2, 2, 2. 2)  position 
2.00000,  2.00000,  2.20000 
0.000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0.000000 
1.00000,  0.000000 


VN:  VOLUMETRIC  PLOT  1  at  Z=2.2m. 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
T 

0.100000,  4.20000,  2.200000 
0.250000,  0.250000,  0.000000 
1,  71,  7 

CELL  7 

SG:  Put  the  source  at  (2, 2, 2. 2)  position 
2.00000,  2.00000,  2.20000 
0.000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0.000000 
1.00000,  0.000000 


VN:  VOLUMETRIC  PLOT  1  at  Z=2.2m. 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
T 

0.100000,  6.20000,  2.200000 
0.250000,  0.250000,  0.000000 
1,  19,  7 

CELL  8 

SG:  Put  the  source  at  (2, 2, 2. 2)  position 
2.00000,  2.00000,  2.20000 
0.000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0.000000 


157 


1.00000,  0.000000 


W:  VOLUMETRIC  PLOT  1  at  Z=2.2  m. 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
T 

7.10000,  6.20000,  2.200000 
0.250000,  0.250000,  0.000000 

I,  7,  7 

CELL  9 

SG:  Put  the  source  at  {2, 2, 2. 2)  position 
2.00000,  2.00000,  2.20000 
0.000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0.000000 
1.00000,  0.000000 

VN:  VOLUMETRIC  PLOT  1  at  Z=2.2m. 

0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
T 

II. 1000,  6.20000,  2.20000 
0.250000,  0.250000,  0.000000 
1,  7,  7 

CELL  10 

SG:  Put  the  source  at  {2, 2, 2. 2)  position 
2.00000,  2.00000,  2.20000 
0.000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0.000000 
1.00000,  0.000000 

VN:  VOLUMETRIC  PLOT  1  at  Z=2.2  tn. 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
T 

15.1000,  6.20000,  2.200000 
0.250000,  0.250000,  0.000000 
1,  11,  7 

CELL  11 

SG:  Put  the  source  at  {2, 2, 2. 2)  position 
2.00000,  2.00000,  2.20000 
0.000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0.000000 
1.00000,  0.000000 
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VN:  VOLUMETRIC  PLOT  1  at  Z=2 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000, 
T 

0.100000,  8.20000,  2.200000 
0.250000,  0.250000,  0.000000 
1,  71,  3 

FREE  SPACE 
CELL  1 


SG:  Put  the  source  at  (2,2,2 
2.00000,  2.00000,  2.20000 
0.000000,  0.000000,  90.0000, 
-1,  0.500000,  0.000000 
1.00000,  0.000000 


VN:  VOLUMETRIC  PLOT  1  at  Z=2 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000, 
T 

0.100000,  1.20000,  2.200000 
0.250000,  0.250000,  0.000000 
1,  71,  3 


CELL  2 

SG:  Put  the  source  at  (2,2,2 
2.00000,  2.00000,  2.20000 
0.000000,  0.000000,  90.0000, 
-1,  0.500000,  0.000000 
1.00000,  0.000000 


VN:  VOLUMETRIC  PLOT  1  at  Z=2 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000, 
T 

0.100000,  2.20000,  2.200000 
0.250000,  0.250000,  0.000000 
1,  19,  7 


CELL  3 


2in. 

0.000000 


2)  position 
0.000000 

2m. 

0.000000 


2)  position 
0.000000 

2m. 

0.000000 
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SG:  Put  the  source  at  (2,2,2 
2.00000,  2.00000,  2.20000 
0.000000,  0.000000,  90.0000, 
-1,  0.500000,  0.000000 
1.00000,  0.000000 


VN:  VOLUMETRIC  PLOT  1  at  Z=2 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000, 
T 

7.10000,  2.20000,  2.200000 
0.250000,  0.250000,  0.000000 
1,  7,  7 


CELL  4 


SG:  Put  the  source  at  (2,2,2 
2.00000,  2.00000,  2.20000 
0.000000,  0.000000,  90.0000, 
-1,  0.500000,  0-000000 
1.00000,  0.000000 

VN:  VOLUMETRIC  PLOT  1  at  Z=2 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000, 
T 

11.1000,  2.20000,  2.200000 
0.250000,  0.250000,  0.000000 
1,  7,  7 


CELL  5 


SG:  Put  the  source  at  (2,2,2 
2.00000,  2.00000,  2.20000 
0.000000,  0.000000,  90.0000, 
-1,  0.500000,  0.000000 
1.00000,  0.000000 


VN:  VOLUMETRIC  PLOT  1  at  Z=2 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000, 
T 

15.1000,  2.20000,  2.200000 


.2)  position 
0.000000 

.2m. 

0.000000 


2)  position 
0.000000 

2  m. 

0.000000 


2)  position 
0.000000 

2m. 

0.000000 
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0.250000,  0.250000,  0.000000 
1,  11,  7 


CELL  6 

SG:  Put  the  source  at  (2,2,2 


2.00000, 

2.00000,  2 

.20000 

0.000000, 

0.000000, 

90.0000, 

-1,  0.500000,  0.000000 

1.00000, 

0.000000 

VN:  VOLUMETRIC  PLOT 

1  at  Z=2 

0.000000, 

0.000000, 

0.000000 

0.000000, 

T 

0.000000, 

90.0000, 

0.100000, 

4.20000, 

2.200000 

0.250000, 
1,  71,  7 

0.250000, 

0.000000 

CELL  7 

SG:  Put  the  source 

at  (2,2,2 

2.00000, 

2.00000,  2 

.20000 

0.000000, 

0.000000, 

90.0000, 

-1,  0.500000,  0.000000 

1.00000, 

0.000000 

VN:  VOLUMETRIC  PLOT 

1  at  Z=2 

0.000000, 

0.000000, 

0.000000 

0.000000, 

0.000000, 

90.0000, 

T 

0.100000, 

6.20000, 

2.00000 

0.250000, 

0.250000, 

0.000000 

1,  19,  7 

CELL  8 

SG:  Put  the  source  at  (2,2,2 
2.00000,  2.00000,  2.20000 
0.000000,  0.000000,  90.0000, 
-1,  0.500000,  0.000000 
1.00000,  0.000000 


2)  position 
0.000000 

2m. 

0.000000 


2)  position 
0.000000 

2m. 

0.000000 


.2)  position 
0.000000 
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VW:  VOLUMETRIC  PLOT  1  at  Z=2.2in. 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
T 

7.10000,  6.20000,  2.200000 
0.250000,  0.250000,  0.000000 
1,  7,  7 


CELL  9 


SG:  Put  the  source  at  (2, 2, 2. 2)  position 
2.00000,  2.00000,  2.20000 
0.000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0.000000 
1.00000,  0.000000 

VN:  VOLUMETRIC  PLOT  1  at  Z=2.2m. 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
T 

11.1000,  6.20000,  2.200000 
0.250000,  0.250000,  0.000000 
1,  7,  7 


CELL  10 


SG:  Put  the  source  at  (2, 2, 2. 2)  position 
2.00000,  2.00000,  2.20000 
0.000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0.000000 
1.00000,  0.000000 

VN:  VOLUMETRIC  PLOT  1  at  Z=2.2ra. 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
T 

15.1000,  6.20000,  2.200000 
0.250000,  0.250000,  0.000000 
1,  11,  7 
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CELL  11 


SG:  Put  the  source  at  (2, 2, 2.2)  position 
2.00000,  2.00000,  2.20000 
0.000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0.000000 
1-00000,  0-000000 


VN:  VOLUMETRIC  PLOT  1  at  Z=2.2  m. 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
T 

0.100000,  8.20000,  2.200000 
0.250000,  0.250000,  0.000000 
1,  71,  3 


EN:  END 

RCX>M  DECOMPOSITION  INTO  CELLS  (FORTH  LOCATION) 
CELL  1 

SG:  Put  the  source  at  (4,3,1)  position 
4.00000,  3.00000,  1.00000 
0.000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0.000000 
1.00000,  0.000000 

VN:  VOLUMETRIC  PLOT  1  at  Z=1  m. 

0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
T 

0.100000,  1.20000,  1.00000 
0.250000,  0.250000,  0.000000 
1,  71,  3 


CELL  2 

SG:  Put  the  source  at  (4,3,1)  position 
4.00000,  3.00000,  1.00000 
0.000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0.000000 
1.00000,  0.000000 


VN:  VOLUMETRIC  PLOT  1  at  Z=1  m. 
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0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
T 

0.100000,  2.20000,  1.00000 
0.250000,  0.250000,  0.000000 
1,  19,  7 

CEUj  3 

SG:  Put  the  source  at  (4,3,1)  position 
4.00000,  3.00000,  1.00000 
0.000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0.000000 
1.00000,  0.000000 
VN:  VOLUMETRIC  PLOT  1  at  Z=lni 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
T 

7.10000,  2.20000,  1.00000 
0.250000,  0.250000,  0.000000 

I,  7,  7 

CELL  4 

SG:  Put  the  source  at  (4,3,1)  position 
4.00000,  3.00000,  1.00000 
0.000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0.000000 
1.00000,  0.000000 

VN:  VOLUMETRIC  PLOT  1  at  Z=lm. 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
T 

II. 1000,  2.20000,  1.00000 
0.250000,  0.250000,  0.000000 
1,  7,  7 

CELL  5 

SG:  Put  the  source  at  (4,3,1)  position 
4.00000,  3.00000,  1.00000 
0.000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0.000000 
1.00000,  0.000000 

VN:  VOLUMETRIC  PLOT  1  at  Z=1  m. 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 


164 


15.1000,  2.20000,  1.00000 
0.250000,  0.250000,  0.000000 
1,  11,  7 

CELL  6 

SG:  Put  the  source  at  (4,3,1)  position 
4.00000,  3.00000,  1.00000 
0.000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0.000000 
1.00000,  0.000000 


VN:  VOLUMETRIC  PLOT  1  at  Z=lm. 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
T 

0.100000,  4.20000,  1.00000 
0.250000,  0.250000,  0.000000 
1,  71,  7 

CELL  7 

SG:  Put  the  source  at  (4,3,1)  position 
4.00000,  3.00000,  1.00000 
0.000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0,000000 
1.00000,  0.000000 


VN:  VOLUMETRIC  PLOT  1  at  Z=lin. 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
T 

0.100000,  6.20000,  1.00000 
0.250000,  0.250000,  0.000000 
1,  19,  7 

CELL  8 

SG:  Put  the  source  at  (4,3,1)  position 
4.00000,  3.00000,  1.00000 
0.000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0.000000 
1.00000,  0.000000 


VN:  VOLUMETRIC  PLOT  1  at  Z=lm. 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
T 
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7.10000,  6.20000,  1.00000 
0.250000,  0.250000,  0.000000 

I,  7,  7 

CELL  9 

SG:  Put  the  source  at  (4,3,1)  position 
4.00000,  3.00000,  1.00000 
0.000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0.000000 
1.00000,  0.000000 

VN:  VOLUMETRIC  PLOT  1  at  Z=lm. 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
T 

II. 1000,  6.20000,  1.0000 
0.250000,  0.250000,  0.000000 
1,  7,  7 

CELL  10 

SG:  Put  the  source  at  (4,3,1)  position 
4.00000,  3.00000,  1.00000 
0.000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0.000000 
1.00000,  0.000000 

VN:  VOLUMETRIC  PLOT  1  at  Z=1  m. 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
T 

15.1000,  6.20000,  1.00000 
0.250000,  0.250000,  0.000000 
1,  11,  7 

CELL  11 

SG:  Put  the  source  at  (4,3,1)  position 
4.00000,  3.00000,  1.00000 
0.000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0.000000 
1.00000,  0.000000 


VN:  VOLUMETRIC  PLOT  1  at  Z=lm. 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
T 

0.100000,  8.20000,  1.00000 
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0.250000,  0.250000,  0.000000 
1,  71,  3 

FREE  SPACE 
CELL  1 


SG:  Put  the  source  at  (4,3,1)  position 
4.00000,  3.00000,  1.00000 
0.000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0.000000 
1.00000,  0.000000 

VN:  VOLUMETRIC  PLOT  1  at  Z=lin. 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
T 

0.100000,  1.20000,  1.00000 
0.250000,  0.250000,  0.000000 
1,  71,  3 


CELL  2 


SG:  Put  the  source  at  (4,3,1)  position 
4.00000,  3.00000,  1.00000 
0.000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0.000000 
1.00000,  0.000000 


VN:  VOLUMETRIC  PLOT  1  at  Z=lm. 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
T 

0.100000,  2.20000,  1.00000 
0.250000,  0.250000,  0.000000 
1,  19,  7 


CELL  3 

SG:  Put  the  source  at  (4,3,1)  position 
4.00000,  3.00000,  1.00000 
0.000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0.000000 
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1.00000,  0.000000 
VN:  VOLUMETRIC  PLOT  1  at  Z=lm. 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
T 

7.10000,  2.20000,  1.00000 
0.250000,  0.250000,  0.000000 
1,  7,  7 


CELL  4 

SG:  Put  the  source  at  (4,3,1)  position 
4.00000,  3.00000,  1.00000 
0.000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0.000000 
1.00000,  0.000000 
VN:  VOLUMETRIC  PLOT  1  at  Z=1  m. 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
T 

11.1000,  2.20000,  1.00000 
0.250000,  0.250000,  0.000000 
1,  7,  7 


CELL  5 

SG:  Put  the  source  at  (4,3,1)  position 
4.00000,  3.00000,  1.00000 
0.000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0.000000 
1.00000,  0.000000 
VN:  VOLUMETRIC  PLOT  1  at  Z=lm. 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
T 

15.1000,  2.20000,  1.00000 
0.250000,  0.250000,  0.000000 
1,  11,  7 


CELL  6 

SG:  Put  the  source  at  (4,3,1)  position 
4.00000,  3.00000,  1.00000 
0.000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0.000000 
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1.00000,  0.000000 
VN:  VOLtMETRIC  PLOT  1  at  Z=lm 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000, 

T 

0.100000,  4.20000,  1.00000 
0.250000,  0.250000,  0.000000 
1,  71,  7 


CELL  7 

SG:  Put  the  source  at  (4,3,1) 
4.00000,  3.00000,  1.00000 
0.000000,  0.000000,  90-0000,  i 
-1,  0.500000,  0.000000 
1.00000,  0.000000 
VN:  VOLUMETRIC  PLOT  1  at  Z=lm 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  I 
T 

0.100000,  6.20000,  1.00000 
0.250000,  0.250000,  0.000000 
1,  19,  7 


CELL  8 

SG :  Put  the  source  at  (4,3,1) 
4.00000,  3.00000,  1.00000 
0.000000,  0.000000,  90.0000,  I 
-1,  0.500000,  0.000000 
1.00000,  0.000000 
VN:  VOLUMETRIC  PLOT  1  at  Z=lm 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  I 
T 

7.10000,  6.20000,  1.00000 
0.250000,  0.250000,  0.000000 
1,  7,  7 


CELL  9 

SG:  Put  the  source  at  (4,3,1) 
4.00000,  3.00000,  1.00000 


.000000 


position 

1.000000 

1-000000 


position 

1.000000 

1.000000 


position 
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0.000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0.000000 
1.00000,  0.000000 
VN:  VOLUMETRIC  PLOT  1  at  Z=lm. 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
T 

11.1000,  6.20000,  1.00000 
0.250000,  0.250000,  0.000000 
1,  7,  7 

CELL  10 

SG:  Put  the  source  at  {4,3,1)  position 
4.00000,  3.00000,  1.00000 
0.000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0.000000 
1.00000,  0.000000 
VN:  VOLUMETRIC  PLOT  1  at  Z=lm. 
0.000000,  0-000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
T 

15.1000,  6.20000,  1.00000 
0.250000,  0.250000,  0.000000 
1,  11,  7 

CELL  11 

SG:  Put  the  source  at  (4,3,1)  position 
4.00000,  3.00000,  1.00000 
0-000000,  0.000000,  90.0000,  0.000000 
-1,  0.500000,  0.000000 
1.00000,  0.000000 

VN;  VOLUMETRIC  PLOT  1  at  Z=1  la. 
0.000000,  0.000000,  0.000000 
0.000000,  0.000000,  90.0000,  0.000000 
T 

0.100000,  8.20000,  1.00000 
0.250000,  0.250000,  0.000000 
1,  71,  3 


EN:  END  OF  FORTH  LOCATION 
End  of  all  NEC -BSC  files 
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APPENDIX  C.  MATHCAD  FILES 

MathCAD  file  used  in  the  con^jarison  with  free  space 
{Con5>lea:  Metallic  Room) 


Plot  of  the  BSC  Results 
Read  the  Data 

NEC_BSC_data  :=  READPRN("C:/Temp/thesis/complex_six  pIate_box_test.|!? 
M  :=  rows(NEC_BSC_data) 

M  =  7.2x10^ 

N  :=  cols(NEC_BSC_data) 

N  =  9 

phi:=  NEC_BSC_data 


NEC_BSC_data 


n:=0..M-i 


14.857  I  20.731  I  15.877  1  18.91  I  18.512  I  18.312  13.591  10.528  26.676  31.365 


N.x  is  the  #  points  in  the  outer  loop 


rows(Etot)  =  7.2  x  10^ 
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Ny:=^ 


rowsi 


(^t) 


N. 


Ny  =  60 
q  :=  0 ..  Ny  -  1 


star^  :=  q-Nx 

stopq  :=  startq  +  Nx-  1 

test  :=  submatri3(Etot ,  staiiq ,  stopq ,  0 ,  o) 

rows(test)  =  120 

cols(test)  =  60 

- > 

Etot_dB  :=  20-lo^Etot) 

test_dB  :=  submatri5(]^ot_dB » startq ,  stopq ,  0 ,  o) 

niax(test_dB)  =  39.564 

nui<test_dB)  =  -195.229 

max_field_dR=  max(test_dB) 

dyn_range:=  23 

field_lin3it=  inax_j5eId_dB-  dyn_range 
p  :=  0..rows(test_dB)  -  1 

fiel4,q  :=  if(test_dBp,q  >  field_limi«est_d^,q,fieid_Iii4l 

test_lin  :=  submatri}(]^t  >  startq ,  stopq ,  0 ,  o) 
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END  of  MathCAD  file  used  in  the  comparison  with  free  space 
(CoxDplex  Metallic  Room) 

MathCAD  file  used  in  the  coaparison  with  free  space  (Room 
With  Two  Antennas) 


Plot  of  the  BSC  Results 
Read  the  Data 

NEC_BSC_data  :=  feADPI^"C:/Temp/thesis/fl_2ant_diff_zaxis 
M  :=  rows(NEC_BSC_data) 

M=  Sxio"^ 


N 

N 


cols(NEC_BSC_data) 


phi  :=  NEC_BSC_data 
phi  = 


0 

1 

2 

3 

4 

'5 

8 

9 

Q 

1 

~ 

1 

1 

1 

1 

1 

1 

1 

1 

^  NEC_BSC_dala 


Er:= 


10 


20 
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.,/b  : 

■:i'V  ^ 

■'  -3  -- 

g.'5^ 

4; 

0.015 

1  1.179 

2.438 

2.097 

1.952 

2.106 

Ee- 

- > 

2s^C_BSC_data 

10 

II 

v>iS' 

"0; 

-V' 

^  3 

8  ^ . '' 

■.5- 

:# 

0.024 

0.029 

0.299 

0.606 

0.813 

1.023 

. . . : _ ^ 

— — ^ 

NEC_BSC_data 

^0 

Eq  :=  10 

1} 

;;  3:^5 

0 

0.438 

0219 

0.406 

0.241 

0.327 

0201 

.  -V 

Eiot  := 


n:=0..M-l 


i';- 

5^0  -i'-; 

1 ' 

^"'2, 

J3.  ^ 

_  8 

0.439 

1.199 

I  2.489 

I  2.196 

2.14 

2.35 

1.115 

2,199 

3.363 

1.961 

rows^Efot)  =  8  X 

rows(Etot) 


Ny  =400 
q:=0..Ny-l 
starlq  :=  q-Nx 
stopq  :=  starlq  +  Nx-  1 

test  :=  submatri:(EK)t ,  starlq ,  stopq ,  0 , 0) 
rows(test)  =  200 
cols(test)  =  400 

Etot_dB  ~  20-lo^Etot) 

test_dB  :=  submatri:(Etot_dB  >  starlq ,  stopq ,  0 ,  o) 
max(test_dB)  =  38.777 
mii<test_dB)  =  -195.229 
max_field_dR=  max(test_dB) 
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When  dynamic  range  is  30  db,  we  get  some  errors,  thus  we 
set  the  following: 
dyn_range  :=  25 

field_limit=  max_field_dB-  dyn_rang€ 
p  :=  0..rows(test_dB)  -  1 

:=  if(test_d]^,q  >  field_limjtest_dBp,q,field_liit^t 

test_lin  :=  submatri)(  Etoi ,  starlq ,  stopq ,  0 ,  o) 
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END  of  yffafchgAP  file  used  in  tlie  coit5>arison  with  free  space 
(Room  With  Two  Antennas) 

MathCAD  file  used  in  the  comparison  with  free  space 
(Symmetric  placement  of  Antenna) 

Plot  of  the  BSC  Results 
Read  the  Data 

MlEC|BSGjiat|pl?iS33K?l^"C:/temp/Ihesis/New_lf^^ 

M  :=  rows(NEC_BSC_data) 

M  =  8xl0'^ 

N  :=  cols(NEC_BSC_<iata) 

N  =  9 

phi:=  NEC_BSC_ciata 
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Ny  =  400 
q:=0..Ny-l 
starlq  :=  q-Nx 
stopq  :=  startq  +  Nx-  1 

test  :=  submatri^Etot, starlq, stopq,0,0) 

rows(test)  =  200 
cols(test)  =  400 
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- > 

Etot_dB  :=  20-Iog{Etot) 

test_dB  :=  submatri)(Eto£_dB » startq ,  stopq ,  0 ,  o) 

ma;^test_dB)  =  65.84 

mii(test_dB)  =  -195.229 

max_field_dR=  max(test_dB) 

dyn_range:=  41 

field_limiJ=  max_field_dB-  dyn_range 
p  :=  0..  rows{test_dB)  -  1 

fiel<|,q  :=  if(test_d]^,q  >  field_limittest_d]^^q,field_lin:^t 

testjin  :=  submatri}(EK,t  s  starlq ,  stopq ,  0 ,  o) 


Ifield 


END  of  MatliCAD  file  used  in  tlie  ccaoparison  with  free  space 
(Symmetric  placement  of  Antenna) 
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MathCAD  file  used  in  the  con?)arison  with  free  space  (Random 
placement  of  Antenna) 


Plot  of  the  BSC  Results 
Read  the  Data 

K1EC_BSC_data:=READPRN("C:/remp/thesis/New_lf_r_op_an_pl_lant4i 
M  :=  rows(NEC_BSC_data) 

M  =  8x 

N  :=  cols(NEC_BSC_data) 

N  =  9 


phi:=NEC_BSC_data 


NEC_BSC_data 


:=  10 

^  0  0.39  0.247  0.252  0.257  0.263  0.268 


n:=0..M-l 


K.x  is  the  #  points  in  the  outer  loop 
rows(Etot)  =  8  X  10^ 


:':r7 

8 

0,281 

0.286 

0.291 

rows(Etot) 

iNy  - — 

Nx 

Ny  =  400 
q:=0..Ny-l 
starlq  :=  q-Nx 
stopq  :=  star^j  +  Nx-  1 

test  :=  submatri:;(Etoi ,  startq ,  stopq ,  0 ,  o) 

rows(test)  =  200 

cols(test)  =  400 

- ^ 

Etot.dB  :=  20-lo|(Etot) 

test_dB  :=  siibmatri5(Etot_dB » starlq ,  stopq ,  0 ,  o) 

inax(test_dB)  =  44.709 
min(test_dB)  =  -195.229 
max_field_dK=  inax(test_dB) 

when  dynamic  range  is  30  db,  we  get  some  errors,  thus  we 
set  the  following: 

dyn_range  :=  33 

field_iimit=  max_fieId_dB-  dynjrange 
p  :=  0..  rows(test_dB)  -  1 

fiel<|,q  :=  if(test_d]^^q  >  field_linijtest_d]^^q,fie!d_iin^t 

testjin  :=  submatri:!(Etot ,  startq ,  stopq ,  0 ,  o) 
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m 


END  of  HathCAD  file  used  in  the 
(Random  placement  of  Antenna) 


18; 


conparison  with  free  space 


file  used  in  the  comparison  with  free  space  (Room 
to  rocaa  scenaurio) 

Plot  of  the  BSC  Results 
Read  the  Data 

NEG_BSC_data  :=  READPRN  (X:/Teii3p/thesis/New_2F_r_op_aiial_pktes.txf 

M  :=  rows(NEC_BSC_data) 

M  =  8xl0'^ 

N  :=  coIs(NEC_BSC_data) 

N  =  9 

phi  :=  NEC_BSC_data 


6.5 1  6.5 


I«C  BSC  data 


MO  -‘JO  I  MO  -10  I  MO  -10  I  MO  -lOjMO  -10|MO  -10 


NEC  BSC  data 


MO  -10  I  MO  -10  I  MO  -10  I  MO  -10|M0  -10|M0  “10 


NEC  BSC  data 


1-10  -10 

ilk 

o 

t 

o 

1-10  -10 

1-10  -10 

1-10  -10 

O 

o 

Etot 

n:=0..M-l 


Six  is  the  #  points  in  the  outer  loop 
rows(E(ot)  =  8  X  10^ 


1.732-10  -10 

1.732-10  -10 

rows(Etot) 

n7~" 


Ny  =  400 
q:=0..Ny-l 


startq  :=  q-Nx 

stopq  :=  starlq  +  Nx-  1 

test  :=  submatri}(Etot ,  startq ,  stopq ,  0 ,  o) 
rows(test)  =  200 

cols(test)  =  400 

- > 

Etot.dB  :=  20-log(Etot) 

test_dB  :=  submatri}(E£ot_dB  .startq, stopq, 0,0) 

max(test_dB)  =  23.308 


mii<test_dB)  =  -195.229 
max_fieId_dB:=  iiiax(test_dB) 

When  dynamic  range  is  30  db,  we  get  some  errors,  thus  we 
set  the  following: 

4yn_range  :=  75 

field_limit=  max_field_dB-  dyn_range 


p  ;=  0..rows(test_dB)  -  1 

fielc^^q  :=  if(test_dBp^q  >  field_lin4test_d^,q,field_lin:^t 
testjin  :=  submatri}(Etot .  startq ,  stopq ,  0 ,  o) 
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END  of  MathCAD  file  used  in  the  comparison  with  free  space 
(Boom  to  rocaa  scenario)  . 


file  used  in  calculating  the  excess  gain/ loss 
random  variables  (First  location) 


NEC3SC_data:=READPRK"daiaroortt'' ) 

Xcoord-  NEC_BSC_data 
y^<5:=NECJBSC_daia^^^ 

Zcoord"  NEC_BSC_data 

antenna  ,  obser. points  NOT  in  the  same  level 

.  .  . > 

P  :=^(Wd  -  14.0)'  ^  (ycoord  “  5.0)^  (^i^ord  “  i-5)^ 


free_space_data  :=  READPRN"<htafiree" ) 


''  NEC_BSC_data  ^  '' 


10 


20 


NEC3SC_d3ta 


<S> 


Ey:=  10 


20 


?gCJSC_d2a 


E,:=10 


20 
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rows(Eto,)  =  1.539X  10^ 


/  «  \ 
Ixee^space^data  ' 


^_frec“*~ 


10 


20 


feee^space^data 


X5) 


^Jree‘~ 


20 


tree^gpace^daza 


<?) 


Ez_firee*“ 


20 


Since  the  HEC-BSC  dynamic  range  is  only  30dB,  we  will  have 
to  accoxmt  for  this  fact  by  "filtering”  out  all  values 
which  are  out  of  that  range.  This  is  done  in  the  following 
fashion : 


jj:=0..rDWs(Etot_<iB)-l 

lowest_allowed_Etot_value  :=  -  dyn_range 

lowest_allowed_Etot_value  =  12.948 


^otJtctjdB  (20‘IOg(Ejot_fre^) 

lowe$t_allowed_Etot_fr_vai  :=  ma:3{Eiocfree_dB)  <iyn_range 


lowest_aiIowed_Etot_fr_vaI  =  7.82 

Etotai  dB..  j  <  Iowest_aliowed_Etot_value  ,  lowest_aIlowed_Etot_value  .E^ot^dB  J 

Etotai  &ee  dB  iffeot  free  dB  <  lowest_.aIlowed3tot^fr...val ,  iowest^allowed JEtoufr^va!  ^E^otJ^^dB.. 

Converting  to  linear: 

- > 

Etotal-10 

- > 

E»t£;_Sxc_dB 

10  ^ 

Let  us  now  create  the  two  data  inatrices,  namely  the 
"datamatrix”  which  consists  of  the  room's  field  strength 
values  for  every  observed  point  and  the  "freespacedataOl” 
which  includes  the  corresponding  data  for  the  free  space 
case. 

datamatrix:=  augment  (p 
freespacedataOl  :=  augment(p , 
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At  this  point,  owing  to  the  nontmiform  layout  of  the 
distance  (from  the  antenna)  points  in  the  □  array,  we  will 
have  to  distinguish  ther,  into  several  different  groups 
using  the  following  approach: 
lower  ;=  floor{iniri(p  )) 
tipper  :=cea(maj(p)) 

wi-in-itnHm  distance  for  this  antenna  location: 

upper  =  15 
lower=ip 

Groups  separation  for  this  calculation  (in  cm) : 

^psamdon  ?=:'5Q 

"Step  factor"  for  this  calculation: 

./  100 

faaor:=ceil  - 

separation 

Ntmtber  of  distance  bins  for  this  calculation: 

bins  :=  upper -factor 

bins  =30 

,  upper  -  lower 

h  := - 

bills 

h  =  0.5 
n  :=  0..  bins 
kjj  :=  lower  -r  n-h 

It  should  be  noted  that  the  nxmber  of  bins  may  be  chosen  in 
a  fashion  which  best  suits  our  desirable  “distance  steps". 
For  instance,  should  we  choose  to  have  groups  points 
separated  by  50can,  where  the  maximum  distance  from  the 
antenna  is  49m,  our  intuitive  selection  should  be  98. 

We  next  introduce  a  new  ftinction  which  shall  be  utilized 
several  times  in  our  subsequent  calculations.  As  input,  it 
esroects  a  two  colxrmn  matrix  and  an  array.  It  looks  for 
those  values  of  the  first  column  (i.e.,  in  our  problem,  the 
distance  array  which  fall  within  our  made-up  groups  and 
then  creates  a  new  matrix  that  contains  all  those  second 
column's  values  whose  respective  first  column  value  satisfy 
our  requirement.  It  finally  returns  this  new  matrix  which, 
it  should  be  tanderlined,  normally  has  a  different  number  of 
elaoents  in  each  row. 
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build(M,binvector)  := 


A  <-  M 

B  <-  M 
k  4—  binvector 
for  l€  O..rows(k)  -  2 
D]  4 —  0 

for  i€  O..rows(k)  -  2 
for  j  €  0..  rows(M)  -  1 
Dj  4-  augmeBt{Bj,Di)  if  Iq  <  Aj  < 
for  i  €  0-  rows(D)  -  1 

Dj  4-  submatrix(Di,0,rows(Di)  -  l,0,coIs(Di)  -  2)  if  Dj  0 
D 

It  is  high  time  we  used  this  helping  function. 
datamatrix02:=  build(datainatxix,  k) 
freespacedata02 

datamatrixoJ  = 

The  following  new  fxxnction  is,  seemingly,  of  high 
importance  since  it  is  the  one  which  calculates  the 
"normalized  difference"  random  variables.  It  receives  two 
matrices  of  the  same  dimensions  and  returns  a  new  matrix 
which  in  each  row  contains  the  values  of  the  random 
variable  pertinent  to  that  row's  "distance  group.” 


:=  build(fi:eespacedata01  ,k) 


;  01  ■■ 

-  1 

-J'; 

;  5'-?, 

■6, 

0 

[1.13] 

[1,35] 

[1.31] 

[1.40] 

[1.67] 

[1,93] 

[1,127] 
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difference(X,Y)  := 


for  k  €  0-  rows(X)  - 1 

Ck-«-^0 

f  for  j.€  0..rows(X)-  l 
n  <r-  cois  (Xj  ) 

Al  4— :  Xj:  -.  ;  . 

B-e-Yj  ' -"K 
i  <-  0 

for  je  O-.n  -  f  K 


A 


r  A. 


diff;  4-  201od 


Aq. 


co!s(X^) 


li 

Cj  4-  Stack (Cj,<iiff) 

k  4-  TOWs(Cj) 

}4-co!s{Cj)  .  . 

Cj:  4-  subiuatrix(Cj,0,k  - 1,1,1-  l) 


oaie  "room  gain"  is  defined  here  in  dB 

We  now  get.  rid  of  the  redtmdant  first  row  zeros  and  then 
calculate  the  normalized  difference  variables. 
new_dataniatrixO2:=subinatiix(dat^atrixO21,rows(data0iatrixO2>  —  2,0,  cols  (datamatrixO?*  —  1) 
new_freespacedaia02  :=submatrix(ffeespacedata02  ,l,rows(freespace<iata02)  -  2,0,cois(fireespacedata02)  -  1) 


norm. 


di#  = 


"  1 " 

^5,;-'  4 

^0' 

[1.35I 

[1,35] 

[1,40] 

[1,67] 

[1,93]  I 

[1,127] 

|[1,lS] 

form  a  bug  for  it  returns  an  erroneous  matrix.  The  cause  of 
this  weird  occurrence  remains  unknown,  however  we  shall 
attenpt  to  effectively  overcome  this  mishap  by  applying  an 
easy  trick, 
i  :=  0..  rows(new_datamatrix02)  -  1 
rj  :=  co!s(new_datamatrixO^) 

corr_norm_diffj  :=  subinatrix(nonn_dlffi,0,rows(norm_diffi)  -  l,0,r;  -  l) 


corr_norm_diff 


1 

1;  '4 

0 

[1,35]  [1,31] 

[1,40] 

[1,67] 

[1,93] 

[1,127] 

[1,134] 

This  last  matrix  named  ”corr_norm_dif f "  is  the  one  we 
intended  to  find.  In  each  row  of  unegual  element  number,  it 
contains  a  different  random  variable  as  function  of  our 
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"devised"  groups.  Let  us  take  a  closer  look  in  one  of  them, 
say  the  5th. 


inean(test_l  q)  =  0.592 

THESE  ARE  %  DIFFERENCE 

Defining  the  number  of  the  bins: 

^ints':=^'2C)0 

^  0..  “  i 

j  :=  0..  rows(corr_nonn_diff)  -  1 

upper  j  :=  ma:{corr_norm_diff  j) 

iowerj  “  miii(corr_norm_diff  j) 

upper_max  :=  ina;(upper) 

lower_mm  ;=  min(lower) 

upper_inax  -  lower_min 
step  := - - - 

upper_inax  =  13.584 
lower jnin  =  -17.064 


step  =0.153 


genjntervals  „  :=  lower_imn  +  step  -ii 


row’s(gen_intervals )  =  200 


buildUL)  := 


•lo 


for  q€  l..rows(L)  -  1 
A  Lq 

4e^,<-stack(4ew.A) 


l^:=buil<mL) 
rowsfl^  =  28 
cols(lp^)  =  199 


189 


qq  I 

avg_pdf  qq  :=  0ieaii(Lpdf  ) 


^avg_pdf  =  0.998 


avg  pdf 


r  " 

■§11 

3.84-10  -4 

0 

0 

0 

3.436-10  -4| 

0 
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is^ypi  of  file  used  in  calculating  tiie  excess 

gain/ loss  random  variables  {First  location) 


WafchCAD  file  used  in  calculating  tbe  excess  gain/loss 
random  variables  {Second  location) 

Calculating  the  normalized  difference  random  variables 
Random  position  STACK 

In  this  program  we  are  aiming  at  calculating  all  the  random 
variables  of  the  normalized  difference  between  the  room's 
and  the  free  space  field  strength  (Etot  and  Etot_free 
respectively)  .We  start  by  reading  the  associated  data  files 
generated  by  NEC_BSC . 

!<iEC_BSC_datal  :=  READPRK'C^^temp^  ) 

Ml  :=rows(NECJBSC_datal) 

Ml  =213 

NI  :=cols(NEC_BSC_datal) 

N1  =  9 

NECJBSCjJata2  :=  lffiADPRK'*C:/ten^p^^ 

M2  :=  rows(NEC_BSC_data3 
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M2  =  133 

N2  :=  cols(NEC_BSC_data3 
N2  =  9 

NEC_BSC_data3  :=  READPRK"C:/temp/Thesis/x3a.txt" ) 
M3  :=rows(NEC_BSC_data3) 

M3  =49 

N3  :=  cols(NEC_BSC_data3) 

N3  =  9 

NEC_BSC_data4;=  READPRK”C:/tempmiesis/x4a.txt" ) 
M4  :=  rows(NEC_BSC_data4) 

M4  =  49 

N4  :=  coIs(NEC_BSC_dara4> 

N4  =  9 

NEC_BSC_data5  :=  READPRI^”C:/temp/Thesis/x5a.txt" ) 
M5  :=  rows(I<iEC_BSC_data5) 

M5  =  77 

N5  :=  cols(NEC_BSC_data5) 

N5  =  9 

NEC_BSC_data6:=  READPKK"C:/temp/Thesis/x6a.txt'' ) 
M6  ;=  rows(NEC_BSC_data^ 

M6  =  497 

N6  ;=  cols  ( NEC_BSC_daia^ 

N6  =  9 

OTC_BSC_data7  :=  READPRN:”C:/temp/Thesis/x7a.txt" ) 
M7  :=  rows(NEC_BSC_data7) 

M7  =  133 

N7  :=  cols(  NEC_BSC_data7) 

N7  =  9 

NEC_BSC_data8  :=  READPRN:"C:/temp/Thesis/x8a.txr  ) 
M8  :=rows(NEC_BSC_data8) 

M8  =  49 

N8  :=  coIs(NEC_BSC_data8) 

N8  =  9 

NEC_BSC_data9  :=  READPRl<^’'C:/temp/rhesis/x9a.txt”  ) 
M9  :=rows(NECJBSC_data9) 

M9  =  49 

N9  :=  coIs(NEC_BSC_data9) 

N9  =  9 

NEC 3SC_datalO;=  READPRN:''C:/temp/Thesis/xlOa.txt”  ) 
MIO  :=  rows(NEC_BSC_datal() 

MiO  =  77 

N10:=  coIs(NEC_BSC_datalQ 
N10  =  9 

NEC_BSC_datal  1  :=  READPRN;"Cy'tempmiesis/xl  la-m" ) 
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Mil  :=  rows(NEC_BSC_dataI  1) 

MU  =  213 

Nil  :=  cols(NEC_BSC_datalI) 

Nil  =  9 

W<i~  NEC_BSC_data 
y«^:=NEC3$C_data^^^ 

2^,d:=NEC_BSC_data^2) 

.  . -  .  . . > 

P  ■—■Ji^xxxrd  ~  -)”  (y coord  ■“  2)“  —  (zcoord  ~  ^)~ 
firee_space_datal  ) 

MMl  :=  rows{free_space_datal ) 

MMl  =  213 

NNl  ;=  cols(free_space_datal ) 

NNl  =  9 

fi'ee_space_;data2  :=READPRl^"Ci/tenip/Thesisifs2a.cct" ) 
MM2  :=  rows(free_space_data2 ) 

MM2  =  133 

NN2  :=  cols(ft'ee_space_data2 ) 

NN2  =  9 

&ee_space_data3  :=  READ‘PRN"C:yienjp/Tliesis/fs3a.txt'^ ) 
MM3  :=  rows(free_space_data3 ) 

MM3  =49 

NN3:=  cols(firee_space_da£a3 ) 

NN3  =  9 

free_space_data4  ?=  READPR]SCX:/teinp/rh^/fs4a;^^^^^ :) 
MM4  :=  rows(free_space_data4  ) 

MM4  =  49 

NN4:=  coIs(firee_space_data4  ) 

NN4=9 

free.spacejdataS :  ^=iaEi!aDI9?N"G:/teiBp/1^ 

MM5  :=  rows(free_space_data5  ) 

MM5  =  T7 

NN5  :=  co!s(free_space_data5 ) 

NN5  =  9 

free_space_data6  :=REM)PRl^''C:/temp/Tliesis/fs6a.txt"  ) 
MM6  :=  rows(free_space_data6 ) 

MM6  =  497 

NN6:=  cols(free_space_data6 ) 

NN6=9 

&ee_space_data7  —  ) 

MM7  :=rows(free_space_data7  ) 

MM7  =  133 


NN7  ?=  cols(free_space_data7 ) 

NN7  =  9 

:free_spaee_data8  :=  B£ADPI^’'C:/temp/Thesis/fs8a.txt"  ) 
MM8  :=  rows(fTee_space_data8  ) 

MM8  =  49 

NN8  :=  cols(free_space_data8 ) 

NN8  =  9 

free_space_data9  :=:RE^PRl<("G;/Emp/lTiesis/fs9  ) 
MM9  :=  rows(free_space_data9  ) 

MM9  =  49 

NN9  :=  cols  (free_space_data9 ) 

NN9  =  9 

free_space_datalO  :=READPRI^"C:/temp/Thesis/fslOa.txt" ) 
MMIO  :=  rows(free_space_datalO ) 

MMIO  =  77 

NN10:=  cols(free_space_datal0  ) 

NN10=9 

free_space_datal  1  :=READPRH”C:/temp/Illiesis/fslla.txt" } 
MMH  :=rows(free_space_datall ) 

MMll  =  213 

NNl  1  ;=  coIs(free_space_datal  1 ) 

NN11  =  9 


NECJBSC.data 


<3)  ^ 


^,:=1  10 


20 


NEC3$C_d3a 


Ej,  :=  10 


20 


NEC_BSC_dat2 


<7) 


Ej:=  10 


20 


row$(Ejot)  ==  1.539X  l(f 


10 


20 


fece^space^daia 


<5) 


20 


frec^space^data 


<7> 


E,w:=10 


20 
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since  the  NEC-BSC  dynamic  range  is  only  30dB,  we  will  have 
to  accoxmt  for  this  fact  by  "filtering"  out  all  values 
which  are  out  of  that  range.  This  is  done  in  the  following 
fashion : 


:=(2(^Iog(Etot)) 


ij:=0..rows(:^ot_dB)  -  1 

Iowest_allowedJEtot_vaIue  :=  ma:i(Eiot_dB)  “  <lyn_range 


lowest_allowedJEtot_value  =  22.107 

- ^ 

EE«_&ee_dB  (20’log(Etjn_fre^) 

lowest_alIowed_Etot_fr_val  :=  nxa3(Eioj_ftee_{jE^  -  dyn_range 


lowest_cdlowedJEtot_fir_vaI  =  22.47 

Etotai.cS.,  <  Iowest_aIlowed_Et>t_vaiue  ,Iowest_a!lowed_Etot_va!iie 

Etotai.freejdB..  “  <  lowest_ailowed_Etot_fr_val ,  lowest_allowed_Etot_fr_vaI  ,Etot_feee_dB.] 

Converting  to  linear: 

- > 

- >  . 

Let  us  now  create  the  two  datamatrices ,  namely  the 
"datamatrix"  which  consists  of  the  room's  field  strength 
values  for  every  observed  point  and  the  "freespacedataOl" 
which  includes  the  corresponding  data  for  the  free  space 
case. 

datainamx:=  augine]3t(p  .Eat^^i) 
freespacedataOl  :=  augment(p ,  ^caai_fe^ 

At  this  point,  owing  to  the  nonunifoima  layout  of  the 
distance  (from  the  antenna)  points  in  the  □  array,  we  will 
have  to  distinguish  them  into  several  different  groups 
using  the  following  approach: 
lower  :=  flooi{mii5(p  )) 
upper  :=  ceil(ma;{p)) 

Maximum  and  minimum  distance  for  this  antenna  locations 


upper  =  17 
lower  =  0 

Groups  separation  for  this  calculation  (In  cm) : 

separatioB  ?=  50 

"Step  factor"  for  this  calculations 
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./  100 

factor  :=  ceil  - 

^  separation 

factor  =  2 

Number  of  distance  bins  for  tbis  calculation: 

bins  :=  upper 'factor 

bins  =  34  : 

.  upper  -  lower 

h  := - 

bins 

h  =  0.5 
n  :=  0..  bins 
kj.  :=  lower  +  n*h 

It  should  be  noted  that  the  number  of  bins  may  be  chosen  in 
a  fashion  which  best  suits  our  desirable  “distance  steps". 
For  instance,  should  we  choose  to  have  groups  points 
separated  by  BOcm,  where  the  maximxom  distance  from  the 
antenna  is  49m,  our  intuitive  selection  should  be  98. 

We  next  introduce  a  new  function  which  shall  be  utilized 
several  times  in  our  subsequent  calculations.  As  input,  it 
expects  a  two  column  matrix  and  an  array.  It  looks  for 
those  values  of  the  first  colximn  (i.e.,  in  our  problem,  the 
distance  array  □..which  fail  within  our  made-up  groups  and 
then  creates  a  new  matrix  that  contains  all  those  second 
coltimn's  values  whose  respective  first  coluinn  value  satisfy 
our  requirement.  It  finally  returns  this  new  matrix  which, 
it  should  be  underlined,  normally  has  a  different  ntmijer  of 
elements  in  each  row. 

buiid(M,binvector)  :=  A  M 
B  M 
k  <-  binvector 
for  le  O..rows(k) -2 
Dj<-0 

for  i  e  0-  rows(k)  -  2 
for  j  €  0..  rows(M)  -  1 
Dj  <-  augment if  Iq  <  A j  < 
for  i  6  0..  rows(D)  -  1 

Dj  <—  submatrix(Di,0,rows(Di)  -  l,0,coIs(Di)  -  2)  if  Dj  0 
D 

It  is  high  time  we  used  this  helping  function. 
datamamx02:=  buiidCdatamatrbc,  k) 
ffeespacedata02 

datamatrixoi^  = 
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:=  build(ffeespacedataOi  ,k) 


The  following  new  fimction  is,  seemingly,  of  high 
importance  since  it  is  the  one  which  calculates  the 
"normalized  difference"  random  variables.  It  receives  two 
matrices  of  the  same  dimensions  and  returns  a  new  matrix 
which  in  each  row  contains  the  values  of  the  random 
variable  pertinent  to  that  row's  “distance  group." 


diff6rence(X,Y)  := 


for  ke  0..  r6ws(X)  -  I 

^for  je  0..rows(X)  -  1 
n  co!s(Xj)  : 

A  4“  Xj  -  ■ 

B-e-Yj 
0 

for  i€  0..n  —  1 
diff; 2(>}62l 


V 


B, 


'0,i 


n  4-  cols(Xj) 

Cj  4-  stack(Cj,diff) 

k  4-  rows(Cj) 

1 4-  cols(Cj) 


Cj  4-  submatrtx^CjjOjk  -1, 1,1  -  l) 


The  "roCTi  gain"  is  defined  here  in  dB 

We  now  get  rid  of  the  redundant  first  row  zeros  and  then 
calculate  the  normalized  difference  variables. 

new_dataniatnx02:=  subniatrix(<iataniatrix02  l,tows{dataniatrixO^  —  2,0,  cols  {datainatrix02i  —  1) 
new_freespacedata02  :=  sabnialiix(freespacedata02 ,  l,rows(freespacedata02)  -  2,0, cols (freespacedata02)  -  1) 
nonn_drff  :=  differeiice(new_dataiiiatrix02,  new_freespacedata02 ) 


diff: 


ymi 

"■',15 

M 

[1,32] 

[1,46]  [1,51] 

[1,51] 

[1,51]| 

[1,51] 

[1,51] 

Surprisingly  and  xinexpectedly  the  Mathcad  appears  to  suffer 
form  a  bug  for  it  returns  an  erroneous  matrix.  The  cause  of 
this  weird  occurrence  remains  unknown,  however  we  shall 
attempt  to  effectively  overcome  this  mishap  by  applying  an 
easy  trick, 
i  :=  0..  rows(new_datamatrix0:5  -  1 
r;  ;=  cols(new_datamatrixO^) 

corr_norm_diffi  :=  submatrix(norau<iiffi,0,rows(nonn_diffi)  -  LO.tj  -  l) 
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corr__norm_diff  =  — - - - - -  . 

0  I  [1 ,32]  I  [1 ,46]  I  [1 ,51]  I  [1 ,47]  |  [1 ,50]  |  [1 ,48]  |  [1 ,48] 

This  last  matrix  named  "corr_norm_diff "  xs  the  one  we 
intended  to  find.  In  each  row  of  xxnecjual  element  number,  it 
contains  a  different  random  variable  as  function  of  our 
"devised"  groups.  Let  us  take  a  closer  look  in  one  of  them, 
say  the  5th. 


1 

2 

,  &■  ■■  '4.'.' 

5... 

6 

7 

8 

6  2.395 

-0.754 

2.494  - 

1.278  -0.751 

-0.465 

0.245 

-2.726 

-0.447 

test_0  :=  vcorr_nonn_diff 

I  0 

teSt_Oo  =  ^  ' . 

^2.395 1  - 

mean(test_0o)  =  0.338 


1 

2 

3 

4 

5 

6 

7, 

8 

1.31 

-0.076 

2.766 

3.854 

-6.99 

2.108 

2.363 

5.748 

test_l  :=  \coiT_norm_diff  } _ _ _ 

p^5  rq  I  2  r^i 

test_lo=  :x_x_____x . .  . . 

2.38 1  1.31  I -0.076 1  2.766 1  3.! 

mean(test_l  o)  =  1-121 

THESE  ARE  %  DIFFERENCE 

Defining  the  nundser  of  the  bins? 

l^points 

ii 0..  Np^jjj*5  “  1 

j  :=  0..  rows(coiT_nomi_<liff)  -  1 
upper  j  :=  ina;(corr_iionn_diff  j) 
lower:  ;=  niin(corT_nonn_diff  j) 
upper_max  ;=  nia;<  upper) 
lower_inin  :=  inin(  lower) 

upper  max  -  lower_min 

^ep:= - - 

■^^pomts 

upper_max  =  13.48 
lower  min  =  -11.543 


step  =  0.125 

gen  Jntervals  :=  lower^min  ^  step  -ii 

rows(gen_intervals  )  =  200 
buildL(L):=  4cw<-l7 

for  q  e  1 ..  rows(L)  -  1 

A«-l/ 

1-sew  ^  Stack  aJ 

Lpef-buildLCD 
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qq;=0..  Npdfs-1 

avg_pdf  :=  mean(Lpdf  ) 


y^avg_pdf  =  I 
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ZVZJpdl  =  ■  \  . '  ■  -  - -  '  "■■■ill', - 

I  0  j5.894-10  -^l  0.474 1  5.^ 

niean_rv:=  ^  gen_intervalsqq-avg_pdf  qq 

qq  =0 

step  =0.125 
avg^pdf  S9  :=  0.005 

^avg_pdf  =0.531 
avg_pdf 

avg_pdf  := - - 

^avg_pdf 

^avg_pdf  =  1 

inean_sq_rv'  :=  ^  (gen_intervals  qq)^-avg_pdf  < 

qq=0 

mean_sq_rv  =  25.604 

var_rv  :=  mean_sq_rv  -  mean^rv'' 

var„rv  =  23,545 

mean_rv^  =  1.435 

Defining  Gaussian  pdf: 

-  (sen_iijxer\’3i%“  meaxt.fv)^ 

1  2'var  TV 

■gauss  ii  :=— . . — -e  ■  -  ■  ^ 


•>^2-7i:-var_rv 

room^ain^  :=  genjntervals 


91  ^  I  f2,  93 

5.734-10  -3  5.304-10  -3  7.015-10  '3 
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Antenna  at  location  2 


END  of  MathCAD  file  used  in  calculating  the  excess 
gain/ loss  random  variables  (Second  location) 


MathCAD  file  used  in  calculating  the  excess  gain/loss 
random  variables  (Third  location) 


NEC_BSC_data:=READPRH"<iataroom_b”  ) 
Xcoord-  NEC_BSC_data 
yeoord-  NEC_BSC_data 
2coot<j:=  NEC_BSC_data 

antenna  ,  obser. points  in  same 

.  .  ;  '  "ryv;  . 

P  -Ji^oord  ~  t  (ycoord  “  (^icooid  “ 
free_space_data  :=READPRK”datafree_b”  ) 


^  NEC_BSC_dasa  ^ 


Ex:= 


10 


20 


NEC_BSC_d2J2 


<s> 


Ev:=  10 


20 


NEC_BSC_data 


<7> 


E^:=  10 


20 


level 
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rows(Etot)  =  1.539X  l(f 


^ee_$pace_dat2  ^ 


^Jree'’~ 


10 


20 


feee_sp2C^_data 


<5> 


20 


free^space^dan 


.<?) 


E..&ee“10 


20 


Since  the  NEC-BSC  dynamic  range  is  only  30dB,  we  will  have 
to  account  for  this  fact  by  "filtering”  out  ail  values 
which  are  out  of  that  range.  This  is  done  in  the  following 
fashion: 

^ot_dB  “  (20'Iog{Etot]j 


jj;=0..rows(EK«_dB)- 1 

lowest_aDowed_Etot_value  :=  max(Eiot_<i0)  -  d>-n_range 

lowest_alIowed_Etot_valBe  =22.115 

-> 


Etot_firee_dB (20-log(Ei(5{_g^j 

lowest_aUowed_Etot_fr_val  :=  maj^]^ot_free_<iB)  ~  <iyn_range 


Iowest_alIowed_Etot_fr_val  =  22.47 

EtotaLdB  ..  dB.,|  <  lowest_allowed_Etot_value  ,  lowest_allowed_Etot_value  ,E^_dB,,l 

Ef^.,1  ftee  ^  --iffeot  free  dB  <  lowest_ailowed_EK>t_fr_vaI ,  lowest_alIowed_Etot_fr_val  ,Etot_&ee  dB  ) 

Converting  to  linear: 

- > 

- > 

BtotaJLfrec.dB 

Wfree-  10  ^ 

Let  us  now  create  the  two  datamatrices,  namely  the 
"datamatrix"  which  consists  of  the  room's  field  strength 
values  for  every  observed  point  and  the  "freespacedataOl" 
which  includes  the  corresponding  data  for  the  free  space 
case. 

datamatrix  —  augment(p  ,Et<^) 
freespacedataOl  :=  augmentfp  ,EAr<.;_faJ 
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At  this  point,  owing  to  the  nonxaniform  layout  of  the 
distance  (from  the  antenna)  points  in  the  □  array,  we  will 
have  to  distinguish  them  into  several  different  groups 
using  the  following  approach: 
lower  :=  flooi{niin(p )) 
upper  :=ceil(maj(p)) 

Maviimim  and  Tnin-irnmn  distance  for  this  antenna  location; 


upper  =  17 
lower  =  0 


Groups  separation  for  this  calculation 

separation  :=  50 

"Step  factor"  for  this  calculation; 

100  ^ 


factor  :=  cei’ 


separation 


(in  cm) : 


factor  =  2 

Humber  of  distance  bins  for  this  calculation; 

bins  :=  upper -factor 


bins  =  34 


,  upper  -  lower 

h  := - 

bins 

h  =  0.5 
n  :=  0..  bins 


ku  :=  lower  -f  n-h 

It  should  be  noted  that  the  number  of  bins  may  be  chosen  in 
a  fashion  which  best  suits  our  desirable  “distance  steps”. 
For  instance,  should  we  choose  to  have  groups  points 
separated  by  50cm,  where  the  maximum  distance  from  the 
antenna  is  49m,  our  intuitive  selection  should  be  98. 

We  next  introduce  a  new  function  which  shall  be  utilized 
several  times  in  our  subsequent  calculations.  As  input,  it 
expects  a  two  coliimn  matrix  and  an  array.  It  looks  for 
those  values  of  the  first  coliimn  (i.e.,  in  our  problem,  the 
distance  array  □. .which  fall  within  our  made-up  groups  and 
then  creates  a  new  m^atrix  that  contains  all  those  second 
column's  values  whose  respective  first  column  value  satisfy 
our  requirement-  It  finally  returns  this  new  matrix  which, 
it  should  be  underlined,  normally  has  a  different.  nunOjer  of 
elements  in  each  row. 
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buiId(M,binvector)  := 


A  <-  M 

B  4-  M 
k  4—  binvector 
for  l€  O..rows(k) -2 
Di4-0 

for  i  €  0..  rows(k)  -  2 
for  j  €  0..  rows(M)  -  I 
Dj  4-  augment (Bj,Dj)  if  k;  <  Aj  < 
for  i  e  0..  rows{D)  -  1 

Di  4-  subiiiatrix(Di,0,rows(Di)  -  l,0,cols(Dj)  -  2)  if  D;  0 
D 

It  is  high  time  we  used  this  helping  ftinction. 
datamatrix02:=  buiid(datainatrix,  k) 
ffeespacedata02 

datamatrixof^  = 

The  following  new  function  is,  seemingly,  of  high 
importance  since  it  is  the  one  which  calculates  the 
"normalized  difference"  random  variables.  It  receives  two 
matrices  of  the  same  dimensions  and  returns  a  new  matrix 
which  in  each  row  contains  the  values  of  the  random 
variable  pertinent  to  that  row's  "distance  group." 


;=  buUd{freespacedata01  ,k) 


4 

,  Qk 

[1.10] 

[1,32] 

i[1,46] 

[1,51] 

[1,47] 

[1,50] 

[1,48] 
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difference(  X,  Y)  := 


for  k  e  0..r6ws(X)  - 1 

Ci;<-0 

f  for  je  0..rows(X)  -  1 
n  cols(Xj) 


B<-  Yj  ■; 

14-0 

for  i  e  0..  n  -  1 


diffj  4-  20'log 


13  ^ —  COiS^Xj^ 

Cj  <—  stack (Cj.diff) 

Cj4-{Cj)^ 

k<-rows(Cj) 

I  4”“  COiS  (Cj} 

Cj  4-  submatrix{Cj,0,k 


The  "room  gain"  is  defined  here  in  dB 

We  now  get  rid  of  the  redimdant  first  row  zeros  and  then 
calculate  the  nontialized  difference  variables- 

new_datamatrix02:=submatrix(datamatrix02i,rows((ktamam  2,0,cois(datamatrix05  -  1) 

new_free$pacedata02  :=submatrix(fTeespacedata02  ,hrows(freespacedata02)  -  2,0, cols (freespacedata02)  -  i) 


nortn. 


difi^  = 


"OM:: 

1 

2 

-■_4'  " 

c§. . 

6 

0 

[1,32] 

f . % 

CD 

[1,51] 

[1,51] 

[1,51] 

[1,51] 

[1,51] 

OUJ- - - - 

form  a  bug  for  it  returns  an  erroneous  matrix.  The  cause  of 
this  weird  occurrence  remains  tonknown,  however  we  shall 
attempt  to  effectively  overcome  this  mishap  by  applying  an 


easy  trick, 
i  :=  0„  rows(new_datamatrix025  -  1 
Ti  :=  cols{new_daiamatrixO!^) 
corr_norm_diffi  := 

T 

corT_norm_diff  = 

This  last  matrix  named  ”corr_norm_dif f ”  is  the  one  we 
intended  to  find.  In  each  row  of  xmegual  element  number,  it 
contains  a  different  random,  variable  as  function  of  our 


submatnx(norm„diffj,0,rows(norm_diff;)  -  1,0,^  l) 


"0 

'1 

2 

4/ 

5 

Q--i 

0 

[1.32] 

[1,46] 

[1,51] 

[1,47] 

[1,50] 

[1,48] 

[1,48] 
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"devised*'  groups.  Let  us  take  a  closer  look  in  one  of  them, 
say  the  5th. 

/  j\  <o> 

test_0  :^\coTtjaonnjijfi  } 

lCStJ>Q  = 

mean(test_Oo)  =  0.197 
te$t_l  := 
test_Io  = 

mean(test_lo)  =  L387 

THESE  ARE  %  DIFFERENCE 

Defining  the  nniaber  of  the  bins: 

^points 

ii  :=  0..  —  i 

j  0.,  rows(coiT_Bonn_diff)  -  I 

upper  j  :=  ma)(corr_norm_diff j) 

loweij  :=  imi3(corr_norm_diiff  j) 

upper_max  max(upper) 

Iower_miii?=  miiKlower) 

upper_max lower^min 

step  :=  - 

N  • 

'^'points 

upper_max  =  14.215 
lower_min  =  -10.339 
step  =0.123 

gen^intervals  ^  :=  lower_min  4  step  ii 


rows(gen_intervals  )  =  200 


bu2ld^L):= 

for  q€  1..  rows(L)  -  1 
A 

^■aew'  ^  Stack 


Lp<if:=buiIdUL) 


cols(Lpdf)  =  199 
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qq:=  O..Npdfs-  1 

avg_pdf  qq  :=  inean(Lp(if  ) 


^avg_pdf  =  0.998 


gen  intervals  = 


'  0,,. 

1 

2 

3 

4 

5 

6 

7 

.0.1 

-10.339 

-10.216 

-10.094 

-9.971 

-9.848 

-9.725 

-9.603 

-9.48 
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avg_pdf 


y>^4'  / 

mrnm 

■86  ’ 

0.463 

15.015-10  -3 

6.468-10  -3 

4.614-10  -3  j 

3.483-10  ”3 

7.664-10  -3 

5.727*10  >3 

avg  pdf  81  :=  0.005 
^avg_pcif  =0.54 

avg_pdf 


avg_pdf  := 


^avg_pdf 
^avg_pdf  =1 


mean_rv  :=  ^  genjntervals  qo_-avg_pdf 
qq=0 

step  =  0.123 

m 


qq 


Npcifs-l 


mean 


.sq_rv  ;=  ^  (gen_intervals  Qq)^-avg_pdf 


qq=0 

mean_sq_rv  =  24.219 

'y 

var_rv  :=  mean_sq_rv  -  mean^rv*" 
var_rv‘  =  15.658 
mean^rv  =  2.926 

Defining  Gaussian  pdf : 


I  -'2-var_rv 

. ™  • . 


>/2'^-yar_rv 
room_gain^  :=  gen  ^intervals 
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•Rwn  of  MathCAD  file  used  in  calculating  the  excess 
gain/loss  random  variables  (Third  location) 

MathCAD  file  used  in  calculating  the  excess  gain/ loss 
random  variables  (Forth  location) 

Calculating  the  normalized  difference  random  variables 
random  position  STACK 

In  this  program  we  are  aiming  at  calculating  all  the  random 
variables  of  the  normalized  difference  between  the  room's 
and  the  free  space  field  strength  {Etot  and  Etot_free 
respectively) -We  start  by  reading  the  associated  data  files 
generated  by  NEC_3SC. 

NEC_BSC_aaral  ;=  READFR!^”C:/££mp/TSeSis/xlc.tx^^ 

Ml  :=  rows(  ]SIEC_BSC_datal) 

Ml  =  213 

Nl  :=  coIs(NEC_BSC_d^l) 

NI  =  9 

NECJBSC_daia2:=READPRI^X:/temp/aesis/x2c.txt") 

M2  :=  rows(NEC_BSC_daia2) 

M2  =  133 

N2  :=  cols(NEC_BSC_data3 
N2  =  9 

1^3SG_.<hta3  ;=  l®i®PE2^"C:^taap/Thesis/^ 

212 


M3  :=rows(NEC_BSC_data3) 

M3  =49 

N3  :=  cols(NEC_BSC_daia3) 

N3  =  9 

NEC3SC_dau4:=  READPRI^”C:/temp/Thesis/x4c.txt" ) 
M4  :=rows(NEC_BSC_data'!^ 

M4  =  49 

N4  :=  cols(NEC_BSC_data4> 

N4  =  9 

NEC_BSC_data5  :=  READPRN;"C:/temp/rhesjs/x5c.txt" ) 
M5  :=rows(NEC_BSC_data5) 

M5  =  77 

N5  :=  coIs(-NEC_BSC_data5) 

N5  =  9 

I^C_BSC_data6  :=  READPiyS[''G;/temp/Thesis/x6c.txr  ) 
M6  :=  rows(NEC_BSC_data@ 

M6  =  497 

N6  :=  cols(NEC_BSC_daia@ 

N6  =  9 

:NEGiJ(SG_jdat^:=:REApPRK''C:/^^ 

M7  :=rows(NEC_BSC_data7) 

M7  =  133 

N7  :=  cols(NEC_BSC_data7) 

N7  =  9 

NEC_BSC_data8  :=  READPRiX”C:/temp/Thesis/x8c.txt”  ) 
M8  :=  rows(NEC_BSC_dataS[) 

M8  =  49 

N8  :=  cols(NEC_BSC_data8) 

N8  =  9 

lsOEC_BSC_data9;=  READPR]S[''C;/temp/Thesis/x9c.txt”  ) 
M9  :=rows(>3EC_BSC_data9) 

M9  =  49 

N9  :=  cols(NEC_BSC_dau9) 

N9  =  9 

NEC_BSC_daialO:=  READPR]Si:"C;/temp/ThesisA;10c.txt" ) 
MIO  :=  rows(NEC_BSC_datalQ 
MIO  =  77 

N10:=  cols(NEC_BSC_datalQ 
N10=9 

NEC_BSC_datal  1  :=  READPRK'’C:/temp/Thesis/xI  Ic.txt" ) 
Ml  1  :=  rows(NEC_BSC_datal  0 
Ml  1  =  213 

N1 1  :=  cols(NEC_BSC_datai  I) 

Nil  =  9 
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5^^:=  NEC_BSC_data 
y«^,4:=NEC_BSC_<kta^^^ 
Zcoord  :=  NEC_BSC_<lata 


P  •sfi^^coofi  ^)”  ~  (y coord”  {?coord“  ^)* 

free^space_<iatal  :=  READPRM] "C:/terap/I1iesis/fslc.txt”  ) 
MMl  :=  rows{free_space_datal ) 

MMl  =213 

NNl  :=  cols(free_space_datal ) 

NNl  =  9 

free_space_<iata2  :=  READPR!S("C:/t^p/T&^is/fs2p-txr  ) 
MM2  :=rows(free_space_data2 ) 

MM2  =  133 

NN2  :=  cols(free_space_data2 ) 

NN2  =  9 

free_space_data3  :=  READPR|^‘!C-j/tdmp/Tliea:^S3c.W^ 
MM3  ;=  rows(free_space_data3 ) 

MM3  =49 

NN3  :=  coIs(firee_space_data3 ) 

NN3  =  9 

free_space_dau4  :=READPRi^”C:/temp/Thesis/fs4c.txt''  ) 
MM4  :=  rows(free_space_data4 ) 

MM4  =  49 

NN4  :=  cols(free_space_data4 ) 

NN4=9 

ftee_space_data5  :=  READPRI^”C:/t^|i/Ilesis/Mc.txr  ), 
MM5  :=  rows(free_space_data5  ) 

MM5  =  77 

NN5  :=  cols(free_space_data5 ) 

NN5  =  9 

&ee_space_data6  ;=  KEADPRN[”C:/temp/lliesis/fs6c.ixt” ) 
MM6  :=  rows{free_space_data6 ) 

MM6  =  497 

NN6;=  cols(free_space_data6 ) 

NN6=9 

free_space_data7  ;=  REAI>PKK^C'ytemp/Thesis/fs7c.m'' ) 
MM7  :=rows(firee_space_da£a7 ) 

MM7  =  133 

NN7  :=  cols(ftee_space_data7  ) 

NN7=9 

free_space_<iata8  ;=  REAI>PPIS^"Gi:/!em|/Iibeas/fe8c.m  ) 
MM8  ;=  rows(fi:ee_space_daia8 ) 
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MM8  =  49 

NN8  ;=  coIs(fTee_space_data8 ) 

NN8  =  9 

free_space_dau9  :=READPRK’'C:/temp/Thesis/ft9c.txt"  ) 
MM9  :=rows(free_$pace_data9 ) 

MM9  =  49 

NN9  :=  cols(free_space_data9 ) 

NN9  =  9 

free_space_datalO  :=  READPR]Si;"C:/temp/Thesis/fslOc.txt" ) 
MMIO  :=  rows(free_space_datalO ) 

MMIO  =  77 

NN10:=  cols(fTee_space_datalO ) 

NN10=  9 

fTee_Space_dataI  1  :=  READPRK"C:/temp/Thesis/fsllc.txt'’  ) 
MM  11  ;=  rows(free_space_datal  1  ) 

MM11  =  213 

NNl  1  :=  cols(free_space_datal  I ) 

NN11  =  9 


^  NEC_BSC_data 


Ex- 


10 


20 


NEC_BSC_d2t3 


<5> 


Ey  :=  10 


20 


NEC_BSC_daia 


<7) 


£2;=  10 


20 


Eot:=^/(Ex)--(Ev)'  +  (E2f 

rows(Etot)  =  1.539X  10^ 


/  <3}  \ 

'  free^space^data  ' 


10 


20 


to^space^data 


/5> 


E>-_free:=  10 


20 


frecjspace^data 


b) 


;_free'' 


10 


20 


Etot_frec’~  (Ey_fre^” 

Since  the  NEC-BSC  d3fnainic  range  is  only  30d3,  we  will  have 
to  account  for  this  fact  by  "filtering"  out  all  values 
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which  are  out  of  that  range, 
fashion: 


[■his  is  done  in  the  following 


EtoudB  •-  (20-Iog(^ot)) 


ij:=0..rows(Eiot_dB)  "  1 

lowest_aUowed_Etot_value  ;=  ina3(^oj_4B)  -  dyn_raiige 
iowest_aIIowed_Etot_value  =  28.39 

^ot_&ee_dB'~  (29'iog(E.(35_§^) 

lowest_allowed_Etot_lT_vaI  :=  -  dyn_range 

lowest_allowed_Etot_fr_vaI  =  28.51 

EtotaLdB.,  :=!:^|^Eiot_dB.J  <^ow®st_ailowed_Etot_value  ,Iowest_allowed_Etot_vaiue  ,Etot_dB. 

Et<«ai  free  dB  :=  free.dB  <  lowest_allowed_Etot_fr_val  ,lowest_alIowed_Etot_fr_val  ,Etot  free  dB,. 
-  “  ii  \  ~  ji  s 

Converting  to  linear: 


^^totaLdS 


20 


EiotaI_£roe_<!B 


20 


Etotal_&ee’“ 

Let  us  now  create  the  two  da tama trices,  nsimely  the 
"datamatrix"  which  consists  of  the  room's  field  strength 
values  for  eveiry  obse2rved  point  and  the  "freespacedataOl" 
which  includes  the  corresponding  data  for  the  free  space 
case. 

datamatrix"  augment(p 
■TeespacedataOi  :==  augmeDt(p  .EA»i_faJ 

At  this  point,  owing  to  the  nomoniform  layout  of  the 
distance  {from  the  antenna)  points  in  the  □  array,  we  will 
have  to  distinguish  them  into  several  different  groups 
using  the  following  approach: 
lower  ;=  floor{min(p )) 
upper  :=  ceil(ma3(p)) 

Ma v-i tmim  and  iwi -n-itmitw  distance  for  this  antenna  location: 
upper  =  15 
lower  =  0 

Groups  separation  for  this  calculation  (in  cm) : 

separation  :=  50 

"Step  factor"  for  this  calculation: 


factor  :=  ceil 


./  100 


separation 

factor  =  2 

Ntnnber  of  distance  bins  for  this  calculation: 
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bins  :=  upper-factor 

bins  =  30 

upper  -  lower 

n  ;= - 

bins 

h  =  0.5 


n  ;=  0..  bins 


kn  :=  lower  T- n-h 

It  should  be  noted  that  the  number  of  bins  may  be  chosen  in 
a  fashion  which  best  suits  our  desirable  "distance  steps". 
For  instance,  should  we  choose  to  have  groups  points 
separated  by  50cm,  where  the  maximtun  distance  from  the 
antenna  is  49m,  our  intuitive  selection  should  be  98. 

We  next  introduce  a  new  function  which  shall  be  utilized 
several  times  in  our  subsequent  calculations.  As  input,  it 
expects  a  two  colximn  matrix  and  an  array.  It  looks  for 
those  values  of  the  first  column  (i.e.,  in  our  problem,  the 
distance  array  □. .which  fall  within  our  made-up  groups  and 
then  creates  a  new  matrix  that  contains  all  those  second 
column's  values  whose  respective  first  coliomn  value  satisfy 
our  requir^ent.  It  finally  returns  this  new  matrix  which, 
it  should  be  underlined,  normally  has  a  different  number  of 
elements  in  each  row. 


build(M ,  binvector  ) 


A  «-  M 
B  M 


k  binvector 


for  le  O..rows(k)  -  2 
D]  4—  0 

for  ie  O..rows(k)  -  2 
for  j  €  0..  rows(M)  -  I 
E^j  <—  augment (BjjDj)  if  Iq  <  Aj  < 
for  i  €  0..  rows(D)  -  1 

Dj  4-  submatri)^Di,0,rows(l>,)  -  l,0,cols(Di)  -  2)  if  Dj  9^  0 
D 

It  is  high  time  we  used  this  helping  function. 
datamatrix02:=  build(datamatrix,  k) 


freespacedata02  :=build(freespacedata01  ,k) 


datamatrixO^ 


- 

:  ,1'  •, 

3 

4 

,  --5  ■ 

6 

0 

[1,14] 

[1.30] 

[1.33] 

[1.63] 

[1,71]l 

[1,76] 

[1,94] 

The  following  new  ftinction  is,  seemingly,  of  high 
importance  since  it  is  the  one  which  calculates  the 
"normalized  difference"  random  variables.  It  receives  two 
matrices  of  the  same  dimensions  and  returns  a  new  matrix 
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which  in  each  row  contains  the  values  of  the  random 
variable  pertinent  to  that  row's  "distance  group." 
difference(X,10:  := 


for  ke  0..rows(X)  -  1 

for  j  e  0-  rowsfX)  ~  1 
n  «-  cols(Xj) 

A  4— .Xj:  : 

B  4-  Yj 
i4-  0 

for  i€  O-h  -  1 
diffj  4-  20-iog| 


Bo. 


n  4- 


Cj  4*^  StfiCl 


k4-  rows(Cj) 
l4-coIs(Cj) 

Cj  4—  subinatrix(Cj,0,k  -  1, 1, 1  -  l) 


The  "room  gain"  is  defined  here  in  dB 

We  now  get  rid  of  the  redundant  first  row  zeros  and  then 
calculate  the  normalized  difference  variables. 

new_dataiiutrix02:=subniatrix(datarnatrix02,l,rows(datamatrix02l  -  2,0,cols(datamatxix0^  -  1) 
new_freespacedata02  :=submatrix(freespacedata02,hrows(ft'eespace<lata02)  -2,0,cols(freespacedaia02)  -  1) 


norm_diff  :=  difference(Hew_<iatamairix02,newjfreespace(iata02 ) 
nonn_difiF  = 


Surprisingly  and  unexpectedly  the  Mathcad  appears  to  suffer 
form  a  bug  for  it  returns  an  erroneous  matrix.  The  cause  of 
this  weird  occurrence  remains  'unknown,  however  we  shall 
atterr^t  to  effectively  overcome  this  mishap  by  applying  an 
easy  trick. 


.  '-'4  ■'* 

WfWS 

[1,30] 

[1,33] 

[1,63] 

[1,71] 

[1,76] 

[1,94] 

[1,99] 

i  :=  0..  row${new_datamatrix02!)  -  1 
n  cois(iiew_datamatrixO^) 
corT_nonn_diffj 

T 

corr_norm_diff"  = 

This  last  matrix  named  "corr_norm_diff "  is  the  one  we 
intended  to  find.  In  each  row  of  unequal  element  nxamber,  it 


subma£rix(norm_diffj,0,rows(nonn_ciififi)  -  l,0,rj  -  l) 


%v'3, . 

?ii># 

[1,30] 

[1,33] 

[1,63] 

1[1,71] 

[1,76] 

[1,94] 

[1.99] 
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contains  a  different  random  variable  as  fnnction  of  our 
"devised”  groups*  Let  us  take  a  closer  look  in  one  of  them, 
say  the  5th . 

te$t_0  := 

testae  0  = 


mean(test_Oo)  =  0.529 

( 

test_I  :=  \conunorm_diff  / 


Q 

:  .tlvi 

■7: 

0 

0.171 

-4.62 

0.389 

-0.842 

-2.145 

3.178 

2.166 

-3.219 

1.823 

1 _ 

mean(test_l  q)  =  '-0-346 

THESE  ARE  %  DIFFERENCE 

Defining  the  nninber  of  the  bins: 

^points 

ii  :=  0..  Npoinis ^ 

j  :=  0..  row$(coiT_norm_diff)  -  1 

upper  j  :=  ma>(corr_norm_diffj) 

iowerj  ;=min(corr„norm_diffj) 

upper_max  :=  ma)<upper) 

lower_miti  :=  min(lower) 

upper  max  -  iower_min 
step  := . . . . 

isj  . 

^^points 


/  (o) 

\corr_norm_diff  / 


1 

2  ’ 

3  4 

6 

7  ■ 

0 

0.078 

-2.592 

0.235 

-0.191 

3.503 

3.51 1 

-10.335 

-0.102 

upper_max  =  12.092 
iower_min  =  -10.335 
step  =0.112 

genjnterv^als  „  :=  lower^min  +  step  -ii 


?  ii  U  '  \  ‘ 


rows(genJntervals  )  =  200 


buildKD  := 


Lneu'^  lo 


for  q€  l.,rows{L)  -  1 

Lnew<-Stack(Lj,e^„A) 


Lp<if:=buildUL) 


cols{Lpdf)  =  199 
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Antenna  ai  location  4 


room^gaxndB.. 


u 


END  of  MathCAD  file  used  in  calculating  the  excess 
gain/ loss  random  variables  (Forth  location) 


End  of  all  MathCAD  files 
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